X=imread('saturn.tif');
imagesc(X);
colormap(gray)
X = double(X);
print -deps saturn
figure
[U,S,V] = svd(X);
for i = 1:4
subplot(2,2,i),plot(U(:,i))
end
print -deps SaturnEvecs