// writing number as a sum of two squares

b := 10^20;
for n in [b .. b + 100] do
    s, x, y := NormEquation(1,n);
    if s then
        printf "%o = %o^2 + %o^2\n", n, x, y;
    end if;
end for;