% Routine 4.2  importing photoimages 
X = imread('canyon.jpg'); 
X = double(X);
Y = (X(:,:,1)+X(:,:,2)+X(:,:,3))/3;
X = Y/3;
X = round(X);
[m,n] = size(X);
image(X)
axis off
colormap gray