-
Gaussian primes en.wikipedia.org/wiki/Gaussian_integer#Gaussian_primes p[x_, y_]:= Polygon[{{x+.5, y+.5}, {x-.5, y+.5}, {x-.5, y-.5}, {x+.5, y-.5}}]; Apply[Show, Join[Map[Graphics, Module[{l, b}, l = {}; Reap[ Do[If[PrimeQ[a + b I], Sow[p[a, b]]], {a, -50, 50}, {b, -50, 50}]]]]]]
-
Original source code: commons.wikimedia.org/wiki/File:Gaussian_primes.png
davidjamesweir’s Twitter Archive—№ 988
