Построить график функции

uses graphABC;
function F(x:real):real;
begin
F:=x*x*x-x*x-x+1;
end;
var xn,xk,x,max,mx,my,dx,x1:real;
    x0,y0,i,k:integer;
    s:string;
begin
x0:=windowwidth div 2;{центр экрана}
y0:=windowheight div 2;
line(0,y0,windowwidth,y0);{оси}
line(x0,0,x0,windowheight);
xn:=-3;xk:=3;{интервал по Х}
mx:=(x0-30)/xk;{масштаб по Х}
if abs(F(xn))>abs(xk)then max:=abs(F(xn))
else max:=(F(xk));
my:=(y0-30)/max;
for i:=1 to 3 do{максимальное количество засечек в одну сторону}
 begin
  line(x0+round(i*mx),y0-3,x0+round(i*mx),y0+3); {засечки на оси Х}
  textout(x0+round(i*mx),y0+10,inttostr(i));
  line(x0-round(i*mx),y0-3,x0-round(i*mx),y0+3);
  textout(x0-round(i*mx),y0+10,inttostr(-i));
  line(x0+3,y0-round(i*my*10),x0-3,y0-round(i*my*10)); {засечки на оси Х}
  textout(x0-30,y0-round(i*my*10)-5,inttostr(i*10));
  line(x0+3,y0+round(i*my*10),x0-3,y0+round(i*my*10));
  textout(x0-30,y0+round(i*my*10)-5,inttostr(-i*10));
  str(i,s);
 end;
{центр}
textout(x0+5,y0+10,'0');
{подписи концов осей}
textout(windowwidth-10,y0-10,'X');
textout(x0-10,10, 'Y');
{график}
x:=xn;
dx:=0.001;
k:=0;
setbrushstyle(bsClear);
setpencolor(clBlue);
textout(x0+80,y0+50,'Ноли функции');
max:=0;
while x<=xk do
 begin
  x:=x+dx; {наращиваем х}
  setpixel(x0+round(x*mx),y0-round(F(x)*my),clRed);
  if abs(F(x))<dx/1000 then
   begin
    inc(k);
    x1:=x;
    setpenwidth(2);
    circle(x0+round(x1*mx),y0,4);
    str(x1:4:1,s);
    textout(x0+100,y0+50+k*20,'x'+inttostr(k)+'='+s);
   end;
 end;
textout(60,20,'График функции F(x)=x^3-x^1-x+1');
end.