This math image could be re-created using vector graphics as an SVG file. This has several advantages; see Commons:Media for cleanup for more information. If an SVG form of this image is available, please upload it and afterwards replace this template with {{vector version available|new image name}}. It is recommended to name the SVG file “Gaussian primes.svg”—then the template Vector version available (or Vva) does not need the new image name parameter.

I, the copyright holder of this work, release this work into the public domain. This applies worldwide. In some countries this may not be legally possible; if so: I grant anyone the right to use this work for any purpose, without any conditions, unless such conditions are required by law.

poly[x_, y_] := 
    Polygon[{{x + 1/2, y + 1/2}, {x - 1/2, y + 1/2},
             {x - 1/2, y - 1/2}, {x + 1/2, y - 1/2}}];

primelist = 
    Module[{l, b}, l = {}; 
      Do[b = 0; 
        While[a^2 + b^2 < 500, If[PrimeQ[a + b I, GaussianIntegers -> True],
            l = Join[l, {poly[a, b], poly[a, -b], poly[-a, b], poly[-a, -b]}]]; 
          b = b + 1], {a, 0, 22}]; l];

gauss = Apply[Show, Join[Map[Graphics, primelist], {AspectRatio -> Automatic}]];

Export["gaussian prime.png", gauss, ImageSize -> 1200];