function test_pca(r)

if nargin == 0; r = 0.5; end

M = imread('cameraman2.jpg');
M = double(M(:,:,1));

[A,B] = my_pca(M,r);

ratio = (prod(size(A))+prod(size(B)))/prod(size(M));
fprintf('Storage used: %4.2f%%\n',ratio*100);

subplot(121); imshow(M,  []); title('Original')
subplot(122); imshow(A*B,[]); title('Condensed')