The function GIdivisor computes the set of divisors which are in the first quadrant of its Gaussian integer argument z.

$\mathrm{with}\left(\mathrm{GaussInt}\right)\:$

$\mathrm{GIdivisor}\left(2+3\mathrm{I}\right)$

$\left\{1\,2+3\mathrm{I}\right\}$

$\mathrm{GIdivisor}\left(3+7\mathrm{I}\right)$

$\left\{1\,1+\mathrm{I}\,3+7\mathrm{I}\,5+2\mathrm{I}\right\}$

$\mathrm{GIdivisor}\left(161+240\mathrm{I}\right)$

$\left\{1\,4+\mathrm{I}\,15+8\mathrm{I}\,52+47\mathrm{I}\,161+240\mathrm{I}\right\}$