Графика: Трехлепестковая роза ρ=sin3φ

uses graph;
var f,r,m:real;
    rz,xc,yc,x,y,i:integer;
    s:string;
begin
initgraph(xc,yc,'');
xc:=getmaxX div 2;{начало координат}
yc:=getmaxY div 2;
rz:=200;{размер картинки}
m:=rz/10;{масштаб для координатной сетки}
f:=0; {угол}
{координатная сетка}
line(xc,0,xc,2*yc);
outtextXY(xc+5,30,'Y');
line(0,yc,2*xc,yc);
outtextXY(2*xc-20,yc-15,'X');
for i:=1 to 10 do
 begin
  str(i/10:0:1,s);
  if i mod 2=0 then
   begin
    line(xc+round(i*m),yc-5,xc+round(i*m),yc+5);
    outtextXY(xc+round(i*m),yc+15,s);
    line(xc-round(i*m),yc-5,xc-round(i*m),yc+5);
    outtextXY(xc-round(i*m),yc+15,'-'+s);
    line(xc+5,yc-round(i*m),xc-5,yc-round(i*m));
    outtextXY(xc-30,yc-round(i*m),s);
    line(xc+5,yc+round(i*m),xc-5,yc+round(i*m));
    outtextXY(xc-35,yc+round(i*m),'-'+s);
   end
  else
   begin
    line(xc+round(i*m),yc-3,xc+round(i*m),yc+3);
    line(xc-round(i*m),yc-3,xc-round(i*m),yc+3);
    line(xc+3,yc-round(i*m),xc-3,yc-round(i*m));
    line(xc+3,yc+round(i*m),xc-3,yc+round(i*m));
   end
 end;
{график}
while f<=2*pi do
 begin
  r:=sin(3*f);
  x:=xc+round(rz*r*cos(f));
  y:=yc-round(rz*r*sin(f));
  putpixel(x,y,12);
  f:=f+0.001
 end;
readln
end.