function imhist(im)
hist(reshape(double(im),1,prod(size(im))),0:255)