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.

Summary

Licensing

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.

Mathematica Program

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];