Построить астроиду

uses graphABC;
var x,y,x0,y0,b:integer;
t:real;
begin
x0:=220;
y0:=220;
line(0,220,640,220);
line(220,0,220,480);
t:=0;
b:=50;
while t<=6.28 do
begin
x:=x0+round(b*cos(t)*cos(t)*cos(t));
y:=y0+round(b*sin(t)*sin(t)*sin(t));
putpixel(x,y,25);
t:=t+0.001;
sleep(1);
end;
end.

uses crt,Graph;
var
grDriver: Integer; grMode: Integer;
ErrCode: Integer;
 x,y,x0,y0,b:integer;
 t:real;
 
begin
grDriver := detect;
InitGraph(grDriver, grMode,'C:\Program Files\TP\BGI');
ErrCode := GraphResult;
if ErrCode = grOk then
begin
x0:=220;
y0:=220;
line(0,220,640,220);
line(220,0,220,480);
t:=0;
b:=50;
while t<=6.28 do
begin
x:=x0+round(b*cos(t)*cos(t)*cos(t));
y:=y0+round(b*sin(t)*sin(t)*sin(t));
putpixel(x,y,25);
t:=t+0.001;
end;
readkey;
end;
CloseGraph;
end.

uses crt,graph;
var i,xc,yc,x1,y1:integer;
    x,y,t,m:real;
    s:string;
    r:longint;
begin
clrscr;
repeat
write('Vvedite radius [50..200] r=');
readln(r);
until r in [50..200];
xc:=0;
initgraph(xc,yc,'');
xc:=getmaxX div 2;
yc:=getmaxY div 2;
m:=(yc-2)/2.4;
Line(0,yc,2*xc,yc);
OuttextXY(2*xc-20,yc-10,'X');
Line(xc,0,xc,2*yc);
OuttextXY(xc+10,10,'Y');
line(0,yc,2*xc,yc);
line(xc,0,xc,2*yc);
t:=0;
while t<=2*pi do
 begin
  x:=r*(sqr(cos(t))*cos(t));
  y:=r*(sqr(sin(t))*sin(t));
  x1:=xc+round(x);
  y1:=yc-round(y);
  putpixel(x1,y1,12);
  t:=t+0.001;
 end;
setfillstyle(6,12);
floodfill(xc,yc,12);
outtextXY(50,40,'Astroida');
outtextXY(50,60,'x=r*cos^3(t)  y=r*sin^3(t)');
str(r,s);
outtextXY(50,80,'pri R='+s);
For i:=1 to 3 do
 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));
  Str(i,s);
  OuttextXY(xc+round(i*m),yc+10,s);
  OuttextXY(xc-round(i*m),yc+10,'-'+s);
  OuttextXY(xc-25,yc+round(i*m),'-'+s);
  OuttextXY(xc-15,yc-round(i*m),s);
 End;
 readln;
readln
end.