intror:= function(p)
return ForAny(Elements(GF(p)), a -> ForAny(Elements(GF(p)),
b-> IsZero(a^2 + b^2) and not (IsZero(a) and IsZero(b))));
end;