a = imread('PSWebImageBitMap0001.bmp',  'bmp');
b = a(:,:,2);

c = histeq(b);

figure(4)
clf
imagesc(c)
colormap(gray)
title('Histogram equalization')
