Построить кривую лемнискаты в полярных координатах

uses crt,graph;
var x0,y0,d,gd,gm,i,x1,y1,fn,fk:integer;
    r,a,f,x,xr,yr:real;
    s:string;
begin
clrscr;
write('a=');readln(a);
fn:=0;
fk:=360;
gd:=0;
initgraph(gd,gm,'');
Setbkcolor(1);
x0:=getmaxX div 2;
y0:=getmaxY div 2;
{КООРДИНАТНАЯ  СЕТКА}
d:=round((y0-20)/(a*sqrt(2)));{шаг по оси радиуса в пикселях}
for i:=1 to round(a*sqrt(2)) do
 begin
  {рисуем окружности}
   circle(x0,y0,i*d);
  {пишем шкалу}
  str(i,s);
  outtextXY(x0+i*d+3,y0+10,s);
 end;
{рисуем касательные пунктиром}
  Setlinestyle(1,0,1);
  line(x0-y0,y0+y0,x0+y0,y0-y0);
  line(x0+y0,y0+y0,x0-y0,y0-y0);
{рисуем осевые линии сплошной линией}
Setlinestyle(0,0,1);
line(10,y0,getmaxX-10,y0);
line(x0,0,x0,getmaxY);
outtextXY(x0+5,y0+10,'0');
outtextXY(getmaxX-15,y0+10,'L');
{ГРАФИК}
x:=fn;
while x<=fk do
 begin
  if cos(2*x)>=0 then
   begin
    r:=a*sqrt(2*cos(2*x));
    xr:=r*cos(x);
    yr:=r*sin(x);
    x1:=x0+round(xr*d);
    y1:=y0-round(yr*d);
    putpixel(x1,y1,14);
   end;
  x:=x+0.01;
 end;
readln
end.

uses crt,graph;
var x0,y0,d,gd,gm,i,x1,y1,fn,fk:integer;
    r,a,f,x,xr,yr:real;
        s:string;
        begin
        clrscr;
        write('a=');readln(a);
        fn:=0;
        fk:=360;
        gd:=0;
        initgraph(gd,gm,'');
        Setbkcolor(1);
        x0:=getmaxX div 2;
        y0:=getmaxY div 2;
 
        d:=round((y0-20)/(a*sqrt(2)));
        for i:=1 to round(a*sqrt(2)) do
         begin
 
              circle(x0,y0,i*d);
 
                  str(i,s);
                    outtextXY(x0+i*d+3,y0+10,s);
                     end;
 
                       Setlinestyle(1,0,1);
                         line(x0-y0,y0+y0,x0+y0,y0-y0);
                           line(x0+y0,y0+y0,x0-y0,y0-y0);
 
                           Setlinestyle(0,0,1);
                           line(10,y0,getmaxX-10,y0);
                           line(x0,0,x0,getmaxY);
                           outtextXY(x0+5,y0+10,'0');
                           outtextXY(getmaxX-15,y0+10,'L');
 
                           x:=fn;
                           while x<=fk do
                            begin
                              if cos(2*x)>=0 then
                                 begin
                                     r:=a*sqrt(2*cos(2*x));
                                         xr:=r*cos(x);
                                             yr:=r*sin(x);
                                                 x1:=x0+round(xr*d);
                                                     y1:=y0-round(yr*d);
                                                         putpixel(x1,y1,14);
                                                            end;
                                                              x:=x+0.01;
                                                               end;
                                                               readln
                                                               end.