DescriptionEisenstein primes.png	Eisenstein primes with norm less than 500
Date	6 December 2007
Source	own work, with Mathematica 6.0
Author	Hack
Permission (Reusing this file)	Public Domain
Other versions	File:Eisenstein primes.svg is a vector version of this file. It should be used in place of this PNG file when not inferior. File:Eisenstein primes.png → File:Eisenstein primes.svg For more information, see Help:SVG. In other languages Alemannisch ∙ Bahasa Indonesia ∙ Bahasa Melayu ∙ British English ∙ català ∙ čeština ∙ dansk ∙ Deutsch ∙ eesti ∙ English ∙ español ∙ Esperanto ∙ euskara ∙ français ∙ Frysk ∙ galego ∙ hrvatski ∙ Ido ∙ italiano ∙ lietuvių ∙ magyar ∙ Nederlands ∙ norsk bokmål ∙ norsk nynorsk ∙ occitan ∙ Plattdüütsch ∙ polski ∙ português ∙ português do Brasil ∙ română ∙ Scots ∙ sicilianu ∙ slovenčina ∙ slovenščina ∙ suomi ∙ svenska ∙ Tiếng Việt ∙ Türkçe ∙ vèneto ∙ Ελληνικά ∙ беларуская (тарашкевіца) ∙ български ∙ македонски ∙ нохчийн ∙ русский ∙ српски / srpski ∙ татарча/tatarça ∙ українська ∙ ქართული ∙ հայերեն ∙ বাংলা ∙ தமிழ் ∙ മലയാളം ∙ ไทย ∙ 한국어 ∙ 日本語 ∙ 简体中文 ∙ 繁體中文 ∙ עברית ∙ العربية ∙ فارسی ∙ +/−

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[ Table[{x - y/2, Sqrt[3] y/2}, {i, 1, 6}] + {{1/2, 1/(2 Sqrt[3])},
            {0, 1/Sqrt[3]}, {-1/2, 1/(2 Sqrt[3])}, {-1/2, -1/(2 Sqrt[3])},
            {0, -1/Sqrt[3]}, {1/2, -1/(2 Sqrt[3])}} ];

eisenprm[x_, y_] :=
    Module[{n, p}, n = x^2 - x*y + y^2;
      If[PrimeQ[n], Return[True]];
      p = Sqrt[n]; If[Mod[p, 3] == 2 && PrimeQ[p], True, False]]

primelist = 
    Module[{l}, l = {}; 
      Do[b = 0; While[a^2 - a*b + b^2 < 500, If[eisenprm[a, b],
            l = Join[l, {poly[a, b], poly[-a, -b], poly[-b, a - b], 
                  poly[b, -a + b], poly[-a + b, -a], poly[a - b, a]}]]; 
          b = b + 1], {a, 0, 22}]; l];

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

Export["eisenstein prime.png", eisen, ImageSize -> 600, 
    ImageResolution -> 144];