16 % This should be exact.
17 actual = rank_k_approximation(A, 3);
19 unit_test_equals("Full rank approximation of a matrix is itself", ...
26 actual = rank_k_approximation(A, 13);
28 unit_test_equals("Full rank approximation of a matrix is itself", ...
36 actual = rank_k_approximation(A, 2);
37 unit_test_equals("Rank 2 approximation of a matrix is correct", ...
39 (norm(actual - expected) <= 1e-10));
45 actual = rank_k_approximation(A, 1);
46 unit_test_equals("Rank 1 approximation of a matrix is correct", ...
48 (norm(actual - expected) <= 1e-10));