Численные методы. Решить нелинейное уравнение методом половинного деления

uses graphABC;
const e=0.0001;
function F(x:real):real;
begin
F:=exp(ln(2)*x)-3*x-2;
end;
var k:array[1..20] of real;
    xn,xk,a,b,c,x,max,min,mx,my,dx,y1:real;
    x0,y0,n,i:integer;
    s:string;
begin
xn:=-5;xk:=5;{интервал по Х}
n:=0;
a:=xn-0.1;
b:=a+0.5;
repeat
//выделяем интервалы корней через 0.5
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;
  n:=n+1;
  k[n]:=(a+b)/2;
 end ;
a:=b;
b:=a+0.5;
until a>xk;
{найдем максимум и минимум для масштаба по У}
x:=xn;
max:=F(xn);
min:=F(xn);
while x<=xk+0.1 do
 begin
  if F(x)>max then max:=F(x);
  if F(x)<min then min:=F(x);
  x:=x+0.1;
 end;
x0:=windowwidth div 2;
mx:=(x0-30)/xk;{масштаб по Х}
y0:=windowheight div 2;
if max>abs(min) then my:=(y0-20)/max{по У}
else my:=(y0-20)/-min;
line(0,y0,windowwidth,y0);{оси}
line(x0,0,X0,windowheight);
for i:=1 to round(max) do{максимальное количество засечек в одну сторону}
 begin
  line(x0-3,y0-round(i*my),x0+3,y0-round(i*my));{засечки на оси У}
  line(x0-3,y0+round(i*my),x0+3,y0+round(i*my));
  {подпись оси У}
  str(i,s);
  if i mod 2=0 then
   begin
    textout(x0-35,y0-round(i*my),s);{соответственно засечкам}
    textout(x0-35,y0+round(i*my),'-'+s);
   end;
  line(x0+round(i*mx),y0-3,x0+round(i*mx),Y0+3); {засечки на оси Х}
  line(x0-round(i*mx),y0-3,x0-round(i*mx),Y0+3);
  {подпись оси Х}
   textout(x0+round(i*mx),y0+10,s);
   textout(x0-round(i*mx),y0+10,'-'+s);
 end;
{центр}
textout(x0+5,y0+10,'0');
{подписи концов осей}
textout(windowwidth-10,y0-10,'X');
textout(x0+5,10, 'Y');
{график}
x:=xn;
dx:=0.001;
while x<=xk do
 begin
  setpixel(x0+round(x*mx),y0-round(F(x)*my),clBlue);
  x:=x+dx;
 end;
setpencolor(14);
setfontsize(11);
setfontcolor(clRed);
for i:=1 to n do
 begin
  setbrushcolor(clRed);
  circle(x0+round(k[i]*mx),y0-round(F(k[i])*my),4);
  str(k[i]:0:4,s);
  setbrushstyle(bsClear);
  textout(x0+round(k[i]*mx)-20,y0+10,s);
 end;
{название}
textout(50,20,'График функции y=2^x-3x-2');
textout(80,40,'на интервале [-5;5]');
end.