Исследование алгебраического уравнения функции и построение его графика

uses graphABC;
const a=-2;//интервал
      b=2;
      e=0.0001;//точность
function f(x:real):real;
begin
f:=-x*x*x*x+2*x*x+3;
end;
function f1(x:real):real;
begin
f1:=-4*x*x*x+4*x;
end;
function Bisec(xn,xk:real):real;
var x,c:real;
begin
repeat
c:=(xn+xk)/2;
if F(xn)*F(c)<0 then xk:=c
else xn:=c;
until xk-xn<e;
Bisec:=(xn+xk)/2;
end;
var x,dx,mx,my,x1,x2,x3,y:real;
    x0,y0,n,i:integer;
    s1,s2:string;
begin
setwindowsize(700,500);
centerwindow;
//начало координат
x0:=windowwidth div 2;
y0:=windowheight div 2;
//оси
line(0,y0,windowwidth,y0);
textout(windowwidth-20,y0-20,'X');
line(x0,0,x0,2*y0);
textout(x0+5,0,'Y');
textout(x0+5,y0+10,'0');
//масштабы по осям
mx:=(x0-20)/2;
my:=(y0-20)/abs(f(a));
//подписи по осям
for i:=1 to trunc(abs(f(a)))+1 do
 begin
  line(x0+round(i*mx/2),y0-3,x0+round(i*mx/2),y0+3);
  line(x0-round(i*mx/2),y0-3,x0-round(i*mx/2),y0+3);
  textout(x0+round(i*mx/2),y0+10,floattostr(i/2));
  textout(x0-round(i*mx/2),y0+10,floattostr(-i/2));
  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));
  textout(x0-30,y0-round(i*my)-5,inttostr(i));
  textout(x0-35,y0+round(i*my)-5,inttostr(-i));
 end;
 //график функции
setpencolor(clBlue);
x:=a;
dx:=0.01;
while x<=b do
 begin
  if x=a then moveto(x0+round(x*mx),y0-round(f(x)*my))
  else lineto(x0+round(x*mx),y0-round(f(x)*my));
  x:=x+dx;
 end;
//график производной
setpencolor(clRed);
x:=a;
dx:=0.01;
while x<=b do
 begin
  if x=a then moveto(x0+round(x*mx),y0-round(f1(x)*my))
  else lineto(x0+round(x*mx),y0-round(f1(x)*my));
  x:=x+dx;
 end;
//нахождение точек экстремумов через корни производной
dx:=0.1;
x:=a;
i:=0;
setfontcolor(clGreen);
setpencolor(clGreen);
if abs(f1(x))<e then
 begin
  i:=i+1;
  str(x:7:3,s1);
  str(f(x):7:3,s2);
  setbrushstyle(bsClear);
  textout(20,i*20,'x'+inttostr(i)+'='+s1+' f(x)='+s2);
  setbrushstyle(bsSolid);
  setbrushcolor(clGreen);
  circle(x0+round(x*mx),y0-round(f1(x)*my),3);
 end;
while x<b do
 begin
  x1:=x;
  x2:=x+dx;
  if f1(x1)*f1(x2)<=0 then
   begin
    x3:=Bisec(x1,x2);
    y:=f1(x3);
    str(x3:7:3,s1);
    str(f(x):7:3,s2);
    i:=i+1;
    setbrushstyle(bsClear);
    textout(20,i*20,'x'+s1+'='+s2);
    setbrushstyle(bsSolid);
    setbrushcolor(clGreen);
    circle(x0+round(x3*mx),y0-round(y*my),3);
   end;
  x:=x+dx;
 end;
end.