Найти все натуральные числа a, b и c из интервала от 1 до 20, для которых выполняется равенство: a*a + b*b = c*c — Pascal(Паскаль)

VAR
  a, c: Integer;
  b: Real;
BEGIN
{Самым простым способом - перебором, хотя есть методы и получше}
  for c := 1 to 20 do
    for a := 1 to c-1 do
      begin
        b := SQRT(c*c - a*a);
        if int(b) = b  then
          WriteLn('a = ',a:3,', b = ',b:3:0,', c = ',c:3,
            ': ',
            a*a:6,' + ',b*b:6:0,' = ',c*c:6);
      end;
END.