2 курс / Численные методы / Задание 5 / Золотое сечение

program zolotoesechenie;
var z,a,b,eps,d,w,q,r: real;
function f(x:real):real;
begin
f:=(x-3)*(x-3)-1+sin(x);
end;
begin
writeln('Функция f=(x-3)^2-1+sin(x)');
r:=(sqrt(5)-1)/2;
writeln('Нужно ввести a,b,eps');
readln(a,b,eps);
writeln('f(a)=',f(a):0:4,' f(b)=',f(b):0:4);
d:=eps/2;
q:=a+(b-a)*r;
w:=b-(b-a)*r;
if f(a)>(f(a+d)) then z:=1 else z:=-1;
while abs(b-a)>eps do
begin
q:=a+(b-a)*r;
w:=b-(b-a)*r;
if (z*f(w))<(z*f(q)) then b:=q else a:=w;
end;
if z=1 then writeln('min=',w:0:4);
if z=-1 then writeln('max=',q:0:4);
readln
end.