clear
close all;  %This closes all open figures

load mandrill
image(X)
colormap(map)
[U,S,V]=svd(X);
NumVecs=[3,6,10,50,150,450];
for k=1:6
    j=NumVecs(k);
    H=U(:,1:j)*S(1:j,1:j)*V(:,1:j)';
    figure(1)
    subplot(2,3,k)
    image(H)
    figure(2)
    subplot(2,3,k)
    imagesc(H-X);
end