Спирограф

uses crt,graph;
type kard=record{тип запись с полями}
          x,y,{точки кардиоиды}
          x1,y1:integer;{центр движущейся окружности}
          end;
mass=array[1..5000] of kard;{массив точек чуть с запасом}
{заполнение массива}
function Nod(m,n:integer):integer;
begin
if m mod n=0 then Nod:=n
else Nod:=Nod(n mod m,m);
end;
procedure Point(cx,cy,a,b,d:integer;var k:mass;var p:integer);
var t:real;
    n:real;
begin
p:=0;{количество точек спирали}
t:=0;{начало рисования}
n:=b/Nod(b,a);
while t<=2*n*pi do{делаем нужное число оборотов}
 begin
  p:=p+1;
  {точки спирали по параметрическому уравнению, сх, су - центр}
  k[p].x:=cx-round((a-b)*cos(t)+d*cos(a*t/b));
  k[p].y:=cy+round((a-b)*sin(t)-d*sin(a*t/b));
  {центр круга вычисляем}
  k[p].x1:=cx-round((a-b)*cos(t));
  k[p].y1:=cy+round((a-b)*sin(t));
  t:=t+0.01;
 end;
end;
{рисование}
procedure Spiral(var k:mass;n,cv,cx,cy,a,b:integer);
var i,j:integer;
begin
i:=1;
while i<n-5 do
 begin
  {рисуем цветом фона - стираем}
  setcolor(0);
  moveto(k[1].x,k[1].y);
  setlinestyle(0,0,1);
  for j:=1 to i do
  lineto(k[j].x,k[j].y);
  setlinestyle(0,0,1);
  Circle(k[j].x1,k[j].y1,b);
  setlinestyle(0,0,3);
  Circle(k[j].x,k[j].y,1);{отверстие}
  i:=i+5;{шаг вперед, можешь поменять}
  {рисуем в цвете}
  setcolor(12);
  setlinestyle(0,0,1);
  moveto(k[1].x,k[1].y);
  for j:=1 to i do
  lineto(k[j].x,k[j].y);{спираль}
  setcolor(8);
  setlinestyle(0,0,1);
  Circle(cx-a,cy,a);{неподвижная окружность}
  Circle(k[j].x1,k[j].y1,b);{подвижная}
  setcolor(12);
  setlinestyle(0,0,3);
  Circle(k[j].x,k[j].y,1);{отверстие}
  delay(50);{задержку тоже можно подобрать самому}
 end;
end;
var x0,y0,a,b,d,cv,n,i:integer;
    t:real;
    k:mass;
begin
clrscr;
{назначем радиус колеса}
repeat
write('Radius kolesa [20..200] a=');
readln(a);
until a in [20..200];
{назначем радиус диска}
repeat
write('Radius diska >=10 <',a,' b=');
readln(b);
until b in [10..a-1];
{расстояние отверстия от центра диска}
repeat
write('Rasstoyanie ot centra diska >=2 <',b,' d=');
readln(d);
until d in [2..b-1];
{переход в графический режим}
x0:=0;
initgraph(x0,y0,'');
{белый фон}
setbkcolor(15);
{центр экрана}
x0:=getmaxX div 2;
y0:=getmaxY div 2;
Point(x0,y0,a,b,d,k,n);
Spiral(k,n,12,x0+a,y0,a,b);
readln
end.