Вот собственно функция.

Я частично код сделал, но он не стабильно работает, помогите пожалуйста.

PROGRAM Tabl;
VAR
x,z, y,y1,y2, xmin, xmax, dx : Real;
BEGIN
Writeln('введите min и max');
Writeln('введите min');
Readln(xmin);
Writeln('введите max');
Readln(xmax);
y1:=xmin+(xmin/0.001);
y2:=(xmax-xmin)/2.83;
if y1<y2 then dx := 0.1
else dx:=0.25;
x :=xmin;
z:=0.564;
WHILE x<=xmax DO begin
y := ((xmin-xmax)*(xmin-xmax)+ln(x))/(1+sqr(sin(z+x)))+ln((xmin-xmax)/x+sqrt(x));
Writeln(x:6:3, y:6:3);
x := x+dx;
end ;
Readln
END.

1) Что это y1:=xmin+(xmin/0.001) ?
2) Почему dx := 0.1 ?
3) В ваших переменных sqrt(x) = sqrt(z)
4) Зачем переименовали переменные?? Проверять тяжело.