load clown.mat;
colormap('gray');
figure(1); image(X);
for k=[200 50 30 10]
  [U,S,V]=svd(X);
  [m,n] = size(S);
  colormap('gray');
  figure(2); image(U(:,1:k)*S(1:k,1:k)*V(:,1:k)');
  colormap('gray');
  title(['fattore di compressione = 1 : ' num2str((m*n)/(m*k+k+k*n))]);
  disp(['fattore di compressione = 1 : ' num2str((m*n)/(m*k+k+k*n))]);
  pause
end;
