Напечатать все пифагоровы числа не привышающие 30 (алгоритм)

Перебираешь все взаимно простые m > n разной четности, так чтобы m*m было меньше 30. И формируешь тройки по формуле:
k*(m*m-n*n, 2*m*n, m*m + n*n)

Например так:
function gcd(a,b: integer): integer;
begin
  if b = 0 then gcd := a
  else if (a > b) and (b > 0) then gcd := gcd(b, a mod b)
  else gcd := gcd(b, a);
end;

procedure Print(m,n:integer);
var i, a, b, c: integer;
begin
  repeat
    i := i + 1;
    a := (m*m-n*n)*i;
    b := 2*m*n*i;
    c := (m*m+n*n)*i;
    if c < 30 then
      writeln(a, ' ', b, ' ', c)
    else break;
  until false;
end;

var m, n, s: integer;
begin
  s := Trunc(Sqrt(30));
  for m := 1 to s do
    for n := m + 1 to s do 
      if (gcd(n, m) = 1) and ((m mod 2)<>(n mod 2)) then
        Print(n, m)
end.