2 курс / Вычислительная физика / Задание 3 / Половинное деление

var x,a,b,ga,al,be,fun,funal,funga,e:real;
begin
writeln('введи нижний предел a');
readln(a);
writeln('введи верхний предел b');
readln(b);
writeln('введи точность e');
readln(e);
al:=a;
be:=b;
fun:=x*x*x*x+2*x*x-x-1;
while be-al>e do
begin;
ga:=(al+be)/2;
funga:=ga*ga*ga*ga+2*ga*ga*ga-ga-1;
funal:=al*al*al*al+2*al*al*al-al-1;
if sign(funga)=sign(funal)then al:=ga else be:=ga;
end;
writeln('ответ равен = ', ga);
end.