lut = uint8(255*(C2/max(C2)));
[Hlut, xlut] = imhist(intlut(hist1, lut),256);
tmpC = cumsum(Hlut);
k = max(tmpC) / max(Hlut);
Clut = tmpC / k;

figure(4);
subplot(4, 1, 1);
imshow(intlut(hist1, lut));
title('LUT converted');

subplot(4,1,2);
histeq(hist1, 256);
title('histeq treated');

subplot(4,1,3);
plot(xlut, Hlut);
title('histogram');

subplot(4,1,4);
plot(xlut, Clut);
title('cumulated histogram');