Any prime number congruent to 1 modulo 4 can be written as a sum of two squares. The pa2b2 command finds such a decomposition.
Example. Input:
Output:
indeed 17=42+12