Построить график функции в декартовой системе координат: y=dctg(kx)

uses graphABC;
function F(x,d,k:real):real;
begin
F:=d*cos(k*x)/sin(k*x);
end;
var xn,xk,max,mx,my,dx,y1:real;
    d,k,x:real;
    x0,y0,i:integer;
    s:string;
begin
d:=0.5;
k:=2;
xn:=-2*pi;xk:=2*pi;{интервал по Х}
x0:=windowwidth div 2;
mx:=(x0-30)/xk;{масштаб по Х}
y0:=windowheight div 2;
my:=(y0-20)/50;{по У}
line(0,y0,windowwidth,y0);{оси}
line(x0,0,x0,windowheight);
for i:=1 to round(xk)+1 do{максимальное количество засечек в одну сторону}
  begin
   line(x0-3,y0-round(i*my*10),x0+3,y0-round(i*my*10));{засечки на оси У}
   line(x0-3,y0+round(i*my*10),x0+3,y0+round(i*my*10));
   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);
   {подпись оси У}
    textout(x0-35,y0-round(i*my*10),s);{соответственно засечкам}
    textout(x0-35,y0+round(i*my*10),'-'+s);
    {подпись оси Х}
    str(i,s);
    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.0001;
while x<=xk do
 begin
  if abs(sin(x))>0.01 then
  setpixel(x0+round(x*mx),y0-round(F(x,d,k)*my),clGreen);
  x:=x+dx;
 end;
{название}
textout(x0+50,20,'График функции y=dctg(kx)');
textout(x0+50,40,'d=0.5 k=2.0');
textout(x0+50,60,'на интервале [-2pi;2pi]');
end.