Построение графика функций в полярной системе координат p^2=a^2*tg(φ)

//p^2=a^2*tg(f) p=a*sqrt(sin(f)/cos(f))
uses graphABC;
var x0,y0,d,r,i,x1,y1:integer;
    f,p,a:real;
    s:string;
begin
x0:=windowwidth div 2; //öåíòð
y0:=windowheight div 2;
a:=5;
//ÊÎÎÐÄÈÍÀÒÍÀß  ÑÅÒÊÀ
d:=round((y0-30)/10);//øàã ïî îñè ðàäèóñà â ïèêñåëÿõ
r:=10*d;//ìàêñèìàëüíûé ðàäèóñ
setbrushstyle(bsClear);
for i:=0 to 10 do
if i mod 2=0 then
 begin
  circle(x0,y0,i*d);//ðèñóåì îêðóæíîñòè
  //ïèøåì øêàëó
  str(i,s);
  textout(x0+i*d-15,y0+5,s);
 end;
f:=pi/6;//øàã ïî êðóãó=30 ãðàä
for i:=1 to 12 do//äåëèì íà 12 ÷àñòåé
 begin
  //ðèñóåì ðàäèóñû ïóíêòèðîì
  setpenstyle(psDash);
  x1:=x0+round(r*cos(i*f));
  y1:=y0-round(r*sin(i*f));
  line(x0,y0,x1,y1);
  //ïèøåì ïîäïèñè
  str(i*30,s);
  case i of
  1..2:textout(x1+5,y1-10,s+#176);
  3:textout(x1-5,y1-15,s+#176);
  4,5:textout(x1-20,y1-10,s+#176);
  6:textout(x1-25,y1-5,s+#176);
  7,8:textout(x1-20,y1,s+#176);
  9:textout(x1-10,y1+5,s+#176);
  10,11:textout(x1+5,y1,s+#176);
  end;
 end;
//ðèñóåì îñåâûå ëèíèè ñïëîøíîé ëèíèåé
setpenstyle(psSolid);
line(x0-r,y0,x0+r,y0);
line(x0,30,x0,windowheight-30);
//ÃÐÀÔÈÊ
f:=0;
setpenwidth(2);
setpencolor(clBlue);
while f<2*pi do
 begin
  if (sin(f)/cos(f))>0 then
   begin
    p:=sqrt(a*a*sqrt(sin(f)/cos(f)));
    setpixel(x0+round(d*p*cos(f)),y0-round(d*p*sin(f)),clBlue);
   end;
  f:=f+0.001;
 end;
setfontcolor(clBlue);
setfontsize(12);
setfontstyle(fsBold);
textout(x0+150,30,'p^2=a^2*tg(f)');
textout(x0+180,50,'a=5')
end.