此math图片可使用矢量图形重新创建为SVG文件。这具有很多好处；更多信息请参见共享资源:待清理媒体。如果存在此图片的SVG格式，请将其上传，然后将此模板替换为{{vector version available|新图片名称}}。 建议您将SVG命名为“Gaussian primes.svg”，这样在使用Vector version available（或Vva）模板时就不需要填写新图片名称参数。

我，本作品著作权人，释出本作品至公有领域。这适用于全世界。 在一些国家这可能不合法；如果是这样的话，那么： 我无条件地授予任何人以任何目的使用本作品的权利，除非这些条件是法律规定所必需的。

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