Вычислить методом половинного деления значение корня уравнения на интервале [-1;0] — Pascal(Паскаль)

uses crt;
function F(x:real):real;
begin
F:=exp(x)+sqrt(1+exp(2*x))-2;
end;
 
var a,b,c,x,e:real;
begin
clrscr;
writeln('F(x)=exp(x)+sqrt(1+exp(2*x))-2');
a:=-1;
b:=0;
repeat
write('Точность: ');readln(e);
until (e>0)and(e<1);
if F(a)*F(b)<0 then
  begin
   repeat
   c:=(a+b)/2;
   if F(a)*F(c)<=0 then b:=c
   else a:=c;
   until b-a<e;
   x:=(a+b)/2;
  end;
writeln('x=',x:0:5);
readln
end.