Спирали Котеса на Pascal ABC

uses graphABC;
var x0,y0,d,i,x1,y1,fn,fk:integer;
    r,a,f,x,xr,yr,m:real;
    s:string;
begin
write('a=');readln(a);
x0:=windowwidth div 2;
y0:=windowheight div 2;
{КООРДИНАТНАЯ  СЕТКА}
d:=round((y0-40)/(round(a)));{шаг по оси радиуса в пикселях}
m:=2*(y0-40);{масштаб для графика}
setpencolor(clBlue);
setbrushstyle(bsClear);
for i:=1 to round(a) do
 begin
  {рисуем окружности}
  circle(x0,y0,i*d);
  {пишем шкалу}
  str(i/a:0:1,s);
  textout(x0-i*d+5,y0+10,s);
 end;
f:=30;{шаг по кругу=30 град}
for i:=1 to 12 do{делим на 12 частей}
 begin
  {рисуем радиусы пунктиром}
  Setpenstyle(psDash);
  x1:=x0+round((y0-40)*cos(i*f*pi/180));
  y1:=y0-round((y0-40)*sin(i*f*pi/180));
  line(x0,y0,x1,y1);
  {пишем подписи}
  str(i*30,s);
  if (x1>=x0)and(y1<y0) then textout(x1+10,y1-15,s)
  else if (x1>=x0)and(y1>=y0) then textout(x1+10,y1+10,s)
  else if (x1<x0)and(y1<y0) then textout(x1-40,y1-15,s)
  else if(x1<x0)and(y1>=y0)then textout(x1-40,y1+10,s)
 end;
{рисуем осевые линии сплошной линией}
Setpenstyle(psSolid);
line(10,y0,windowwidth-10,y0);
line(x0,0,x0,windowheight);
textout(x0+5,y0+10,'0');
textout(25,y0+10,'A');
{ГРАФИК}
x:=0;{начальный угол}
while x<=360 do{полный оборот}
 begin
    r:=cos(a*x);{уравнение в полярных координатах}
    xr:=r*cos(x);{приведение к Декартовым}
    yr:=r*sin(x);
    x1:=x0+round(m*xr/2);{приведение к экранным}
    y1:=y0-round(m*yr/2);
    setpixel(x1,y1,clRed);
    x:=x+0.01;
 end;
end.