Эпициклоида

uses graph,crt;
var gd,gm,a,b,p,q,x,y:integer;
    t,k:real;
{функция наибольший общий делитель для сокращения чисел a,b}
function nod(p,q:integer):integer;
var  z:integer;
begin
repeat
if p>q then p:=p mod q
else q:=q mod p;
until(p=0)or(q=0);
nod:=p+q;
end;
 
begin
clrscr;
write('Vvedite a='); readln(a);
write('Vvedite b='); readLn(b);
{выберем большее для назначения 2*q*pi}
if a>b then
 begin
  x:=a;
  a:=b;
  b:=x;
 end;
{сокращаем}
p:=a div nod(a,b);
q:=b div nod(a,b);
gd:=detect;
initgraph(gd,gm,'');
k:=(getmaxY div 2)/(2*a+b);{масштабный коэффициент}
t:=0;
while t<=2*q*pi do
 begin
  x:=round(k*((b+a)*cos(t)-a*cos(t*(b+a)/a)));
  y:=round(k*((b+a)*sin(t)-a*sin(t*(b+a)/a)));
  putpixel(getmaxX div 2+x,getmaxY div 2-y,14);
  t:=t+0.001;
 end;
readLn
end.