Интервалы поиска корня

Program task_1;
uses crt;
function F(a,b,c,x:real):real;
begin
F:=exp(-a*x)-b;
end;
procedure Dihot(a,b,c,xn,xk,e:real;var x:real);
var xs:real;
begin
repeat
 xs:=(xn+xk)/2;
 if F(a,b,c,xn)*F(a,b,c,xs)<=0 then xk:=xs
 else xn:=xs;
until xk-xn<e;
x:=(xn+xk)/2;
writeln('method dilennya navpil:');
writeln('x=',x:0:1);
end;
procedure Iter(a,b,c,xn,xk,e:real;var x:real);
var x1,x2,f1,f2:real;
    n,i:longint;
begin
{method iteracij}
n:=round((xk-xn)/e);
x1:=xn; f1:=F(a,b,c,x1);
for i:=1 to n do
  begin
   x2:=x1+(xk-xn)/n;
   f2:=F(a,b,c,x2);
   if f1*f2 <= 0 then
    begin
     writeln('method iteracij:');
     writeln('x=',(x2 + x1)/2:0:5);
    end;
   x1:=x2; f1:=f2;
  end;
end;
 
var a,b,c,xn,xk,x,e:real;
begin
clrscr;
writeln ('vvesti znachennya a= , b= ');
readln (a,b);
writeln ('vvesti meghi xn= xk= ');
readln (xn, xk);
if F(a,b,c,xn)=0 then
  begin
     writeln('x=xn=',a:9:6);
     readln;
     exit;
  end;
if F(a,b,c,xk)=0 then
   begin
     writeln('x=xk=',b:9:6);
     readln;
     exit;
   end;
if F(a,b,c,xn)*F(a,b,c,xk)>0 then
  begin
   writeLn('meghi intervaly zadani nevirno');
   readln;
   exit;
  end;
writeln('F(x)=exp(-a*x)-b');
writeln('a=  b=');
writeln;
readln;
readln;
writeln('tochnist= ');
readln(e);
Dihot(a,b,c,xn,xk,e,x);
Iter(a,b,c,xn,xk,e,x);
writeln;
writeln('tochnist= ');
readln (e);
Dihot(a,b,c,xn,xk,e,x);
Iter(a,b,c,xn,xk,e,x);
writeln;
writeln('tochnist= ');
readln (e);
Dihot(a,b,c,xn,xk,e,x);
Iter(a,b,c,xn,xk,e,x);
readln
end.