Программирование с использованием методов одномерной оптимизации - Pascal

program lab1; 
const dx=0.0001;
function f(x:real):real; 
begin
if x=0 then f:=-1 else 
if (x=1)or(x=-1) then f:=-1 else if (x<1)and(x>-1) then f:=5*(exp(ln(e)*(-x))+4*x-((exp(x)*3))/3)
else f:=5*(exp(ln(e)*(-x))+4*x-((exp(x)*3))/3)
end;
function f1(x:real):real; 
begin 
f1:=(f(x+dx)-f(x))/dx; 
end;
function f2(x:real):real; 
begin 
f2:=(f1(x+dx)-f1(x))/dx; 
end; 
var 
h,x,x1,a,b:real; 
begin 
writeln('введите х'); 
read(x); 
writeln('введите шаг поиска отрезка'); 
read(h); 
if (f(x-h)>=f(x))and(f(x)<=f(x+h)) then begin a:=x+h; b:=x-h;end else
begin
if (f(x-h)>=f(x))and(f(x)>=f(x+h)) then x1:=x+h else 
if (f(x-h)<=f(x))then begin x1:=x-h;h:=(-1)*h;end; 
 
while f(x1)<=f(x) do 
begin 
x:=x1; 
x1:=x+h; 
writeln('x=',x,' x1=',x1); 
end; 
writeln('h=',h); 
if h>0 then begin a:=x+h;b:=x-h;end 
else begin a:=x-h;b:=x+h;end; 
end;
writeln('a=',a:0:3,' b=',b:0:3); 
x:=a; 
while f1(x)>0.00001 do 
begin 
x:=x-(f1(x)/f2(x)); 
writeln(x,'  ',f1(x));
end; 
writeln('точка min=',x:0:4); 
end.