Программа рисования гипоциклоиды

uses crt,graphabc;
var x2,y2,y3,r1,r2,a,ws, hs,centerY,centerX,x1,y1:integer;
i,n,b,x3,j:real;
s,t,h:string;
begin
readln(r1,r2);
y3:=r1-r2;
SetWindowSize(1980,1080);
ws:=WindowWidth;
hs:=WindowHeight;
centerX:=round(ws/2);
centerY:=round(hs/2);
i:=r1/r2;
Circle(centerX,centerY,r1);
for a:=2 to 5000 do
begin
x1:=round(r2*(i-1)*(cos(a)+cos((i-1)*a)/(i-1)))+centerx;
y1:=round(r2*(i-1)*(sin(a)-sin((i-1)*a)/(i-1)))+centery;
Circle(round(x1),round(y1),1);
end;
end.

uses graphABC;
type point=record
           x,y:integer
           end;
var xc,yc,x,y,x1,y1,r,r1,n,i:integer;
    u:real;
    a:array[1..2000]of point;
begin
xc:=windowwidth div 2;
yc:=windowheight div 2;
r:=yc-30; //радиус внешней окружности
r1:=r div 3;//радиус внутренней окружности в 3 раза меньше
//составляем массив точек кривой
u:=0;
i:=0;
while u<=6*pi do //по обороту на каждую сторону
 begin
  i:=i+1;
  //параметрическое уравнение гипоциклоиды при r/r1=3
  a[i].x:=xc+round(2*r/3*cos(u/3)+r/3*cos(2*u/3));
  a[i].y:=yc-round(2*r/3*sin(u/3)-r/3*sin(2*u/3));
  u:=u+0.01
 end;
//рисуем
u:=0;
lockdrawing;
n:=0;
while u<=6*pi do
 begin
  clearwindow;
  setpencolor(clBlack);
  circle(xc,yc,r); //внешняя окружность
  x:=xc+round((r-r1)*cos(u/3));//координаты внутренней окружности
  y:=yc-round((r-r1)*sin(u/3));
  circle(x,y,r1);//рисуем ее
  n:=n+1;//добавляем точку
  line(a[n].x,a[n].y,x,y);//рисуем прямую линию
  setpencolor(clRed);
  moveto(a[1].x,a[1].y);//рисуем кривую
  for i:=1 to n do lineto(a[i].x,a[i].y);
  sleep(10);//задержка, можно уменьшить
  u:=u+0.01;//наращиваем угол
  redraw //перерисовываем
 end
end.

uses graphABC;
var xc,yc,x,y,r,r1:integer;
    u:real;
begin
xc:=windowwidth div 2;
yc:=windowheight div 2;
r:=yc-30;
r1:=r div 3;
circle(xc,yc,r);
u:=0;
while u<=6*pi do
 begin
  x:=xc+round(2*r/3*cos(u/3)+r/3*cos(2*u/3));
  y:=yc-round(2*r/3*sin(u/3)-r/3*sin(2*u/3));
  setpixel(x,y,clRed);
  u:=u+0.001
 end;
end.

uses graphABC,crt;
type point=record
           x,y:integer
           end;
function NOD(a,b:integer):integer;
begin
if a=b then NOD:=a
else if a>b then NOD:=NOD(a-b,b)
else NOD:=NOD(b-a,a);
end;
var xc,yc,x,y,x1,y1,r,r1,n,nd,k,i:integer;
    m,u:real;
    a:array[1..10000]of point;
begin
//вводим радиусы целые числа, сократимые до взаимно простых, типа r1/r=1/3;1/4;
//2/3;3/4; 3/5 и т.д.
write('Радиус внешней окружности r=');
read(r);
write('Радиус внутренней окружности r1=');
read(r1);
hidecursor;
xc:=windowwidth div 2;
yc:=windowheight div 2;
m:=r1/r;//коэффициент для уравнений
nd:=NOD(r,r1); //нод
k:=r div nd;//знаменатель соеращеной дроби
//составляем массив точек кривой
u:=0;
i:=0;
while u<=2*k*pi do //по обороту на каждую сторону
 begin
  i:=i+1;
  //параметрическое уравнение гипоциклоиды при r/r1=3
  a[i].x:=xc+round((r-m*r)*cos(m*u)+m*r*cos(u-m*u));
  a[i].y:=yc-round((r-m*r)*sin(m*u)-m*r*sin(u-m*u));
  u:=u+0.01
 end;
//рисуем
u:=0;
lockdrawing;
n:=0;
while u<=2*k*pi do
 begin
  clearwindow;
  setpencolor(clBlack);
  circle(xc,yc,r); //внешняя окружность
  x:=xc+round((r-r1)*cos(m*u));//координаты внутренней окружности
  y:=yc-round((r-r1)*sin(m*u));
  circle(x,y,r1);//рисуем ее
  n:=n+1;//добавляем точку
  line(a[n].x,a[n].y,x,y);//рисуем прямую линию
  setpencolor(clRed);
  moveto(a[1].x,a[1].y);//рисуем кривую
  for i:=1 to n do lineto(a[i].x,a[i].y);
  sleep(10);//задержка, можно уменьшить
  u:=u+0.01;//наращиваем угол
  redraw //перерисовываем
 end
end.

a:array[1..100000]of point;

 u:=u+0.1