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

uses graph;
function F(x:real):real;
begin
F:=x*x*x+2*x*x+x;
end;
var xn,xk,max,min,x,y,mx,my,dx:real;
    x0,y0,i,px,py:integer;
    s:string;
begin
x0:=0;
initgraph(x0,y0,'');
x0:=getmaxX div 2;{центр экрана}
y0:=getmaxY div 2;
xn:=-5;xk:=5;{интервал по Х}
mx:=620/(xk-xn);{масштаб по Х}
{найдем максимум и минимум для масштаба по У}
max:=F(xn);min:=max;
x:=xn;
while x<=xk do
 begin
  x:=x+0.01;
  if F(x)>max then max:=F(x)
  else if F(x)<min then min:=F(x);
 end;
if max>abs(min) then my:=(y0-20)/max{по У}
else my:=(y0-20)/abs(min);
setbkcolor(15);{цвет фона}
setcolor(8);{цвет линий}
line(0,y0,getmaxX,y0);{оси}
line(x0,0,X0,getmaxY);
for i:=1 to 20 do{максимальное количество засечек в одну сторону}
  begin
    line(x0-3,y0+round(10*i*my),x0+3,y0+round(10*i*my));{засечки на оси У}
    line(x0-3,y0-round(10*i*my),x0+3,y0-round(10*i*my));
    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);
    str(i*10,s);
    {подпись оси У}
    outtextXY(x0-35,y0+round(10*i*my),'-'+s);{соответственно засечкам}
    outtextXY(x0-35,y0-round(10*i*my),s);
    {подпись оси Х}
    str(i,s);
    outtextXY(x0+round(i*mx),y0+10,s);
    outtextXY(x0-round(i*mx),y0+10,'-'+s);
  end;
{центр}
outtextXY(x0+5,y0+10,'0');
{подписи концов осей}
outtextXY(getmaxX-10,y0-10,'X');
outtextXY(x0-10,10, 'Y');
setcolor(8);
{график}
x:=xn;
dx:=0.001;
while x<=xk do
   begin
     x:=x+dx; {наращиваем х}
     px:=x0+round(x*mx); {изменение по оси Х}
     y:=F(x); {график в реальных координатах}
     py:=y0-round(y*my);{в экранных}
     putpixel(px,py,8);
   end;
readln
end.

initgraph(x0,y0,'');