Задача на метод деления пополам и метод итераций

uses crt;
function F(a,b,c,x:real):real;
begin
F:=b*(1-exp(-a*x))-c;
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('Метод деления отрезка пополам:');
writeln('x=',x:9:6);
end;
procedure Iter(a,b,c,xn,xk,e:real;var x:real);
var x1,x2,f1,f2:real;
    n,i:longint;
begin
{Метод итераций}
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('Метод итераций:');
     writeln('x=',(x2 + x1)/2:9:6);
    end;
   x1:=x2; f1:=f2;
  end;
end;
 
var a,b,c,xn,xk,x,e:real;
begin
clrscr;
a:=0.2;b:=2;c:=1;
xn:=-1;xk:=5;
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('Границы интервала заданы некорректно!');
   readln;
   exit;
  end;
writeln('F(x)=b*(1-exp(-a*x))-c');
writeln('a=0.2; b=2; c=1');
writeln;
writeln('Точность = 0.1: ');
e:=0.1;
Dihot(a,b,c,xn,xk,e,x);
Iter(a,b,c,xn,xk,e,x);
writeln;
writeln('Точность = 0.05: ');
e:=0.05;
Dihot(a,b,c,xn,xk,e,x);
Iter(a,b,c,xn,xk,e,x);
writeln;
writeln('Точность = 0.001: ');
e:=0.001;
Dihot(a,b,c,xn,xk,e,x);
Iter(a,b,c,xn,xk,e,x);
readln
end.