Графика Паскаль. Параметрические уравнения

uses graph;
var gd,gm:integer;
 
function step(a:real;b:integer):real;
var i:byte;
s:real;
begin
step:=1;
for i:=1 to b do
s:=s*a;
Step:=s;
end;
 
procedure Cissoida(a:real);
var
   x,y:real;
   t:real;
begin
t:=0;
 
 x:=(a*step(t,2))/(1+sqr(t));
 y:=(a*step(t,3))/(1+sqr(t));
 moveto(round(x)+300,round(y)+300);
     while t<=2*pi do
         begin
         x:=(a*step(t,2))/(1+step(t,2));
         y:=(a*step(t,3))/(1+step(t,2));
         putpixel(round(x)+300,round(y)+300,4);
         t:=t+0.001;
         end;
 
end;
begin
gd:=detect;
initgraph(gd,gm,' ');
setbkcolor(15);
setcolor(8);
line(300,0,300,1000);
line(0,300,1000,300);
Cissoida(1.6);
readln;
closegraph;
end.

uses graph;
procedure Cissoida(xc,yc:integer;a:real);
var  x,y,t:real;
begin
t:=-2*pi;
while t<=2*pi do
 begin
  x:=a*t*t/(1+t*t);
  y:=a*t*t*t/(1+t*t);
  putpixel(xc+round(x),yc-round(y),4);
  t:=t+0.001;
 end;
 
end;
var gd,gm,xc,yc:integer;
begin
gd:=detect;
initgraph(gd,gm,' ');
setbkcolor(15);
setcolor(8);
xc:=getmaxX div 2;
yc:=getmaxY div 2;
line(0,yc,2*xc,yc);
line(xc,0,xc,2*yc);
Cissoida(xc,yc,300);
readln;
end.