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

uses graphABC;
function F(x:real):real;
begin
F:=cos(x);
end;
var xn,xk,x,y,mx,dx:real;
    x0,y0,my,i,px,py:integer;
    s:string;
begin
x0:=windowwidth div 2;{центр экрана}
y0:=windowheight div 2;
xn:=-3;xk:=3;{интервал по Х}
mx:=(windowwidth-60)/(xk-xn);{масштаб по Х}
my:=y0-40;{по У}
setpencolor(clBlack);{цвет линий}
line(0,y0,windowwidth,y0);{оси}
line(x0,0,X0,windowheight );
line(x0-3,y0+my,x0+3,y0+my);
textout(x0-15,y0+my,'-1');
line(x0-3,y0-my,x0+3,y0-my);
textout(x0-15,y0-my,'1');
for i:=1 to round(4*pi)+2 do{максимальное количество засечек в одну сторону}
  begin
    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,s);
    textout(x0+round(i*mx)-15,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-10,10, 'Y');
//setpencolor(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);{в экранных}
     setpixel(px,py,clRed);
   end;
textout(60,30,'График функции y=cosx');
end.