[Pascal] Циклические вычислительные процессы

[1Sol]

Всем доброго времени суток.
Прошу помощи в решении задачи:
Ввести целое N , если N>0, то найти все натуральные решения
уравнения k*k=x*x+y*y, если x,y,k меняются в диапазоне натуральных
чисел от 1 до N . Например, для N=10:
3, 4, 5
4, 3, 5
6, 8, 10
8, 6, 10
Если быть точным, прошу объяснить, как сделать так, чтобы программа находила 3 неизвестных из уравнения и диапазона.

[Mad_Cat]

for  k:=1 to n do
for x:=1 to n do
for y:=1 to n do
if sqr(k)=sqr(x)+sqr(y) then writeln(x:5,y:5,k:5);

[1Sol]

премного благодарен

[Denis0077]

Сорри за беспокойствие,но есть аналогичный вопрос.
Ввести целое k , если k > 0, то найти все пары натуральных чисел,
сумма квадратов которых равна числу k . Например, k = 25:
3 4
4 3
Как-то так...

[1Sol]

да почти также:
for x:=1 to k do
for y:=1 to k do
if k=sqr(x)+sqr(y) then
writeln(x:5,y:5);