Чтото тут с формулой не то, потому что в отладчике после
for i:=1 to MAX do inc(Hist[1+trunc((Data[i]-beta+lambda/2)/step)]);
массиd Hist всё равно содержал нулевые элементы
и по этому после for i:=1 to N do if Max_Num<Hist[i] then Max_Num:=Hist[i];
Max_Num было = 0;
и в конечном счете в строке
rectangle(
round(i*GetMaxX/N),
round(2*GetMaxY/3),
round((i+1)*GetMaxX/N),
round(2*GetMaxY/3*(1-Hist[i+1]/Max_Num)));
вылетала ошибка при делении на 0. но почему-то код ошибки был левый.