clf;
im=readp5('lena.pgm');
subplot(2,2,1);
image(im); colormap(gray(256)); axis('image'); axis('off');
title('lena.pgm');

[h ind] = comp_hist(im);
subplot(2,2,2);
plot(ind,h);
