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 ∙ العربية ∙ беларуская (тарашкевіца) ∙ български ∙ বাংলা ∙ català ∙ нохчийн ∙ čeština ∙ dansk ∙ Deutsch ∙ Ελληνικά ∙ English ∙ British English ∙ Esperanto ∙ español ∙ eesti ∙ euskara ∙ فارسی ∙ suomi ∙ français ∙ Frysk ∙ galego ∙ Alemannisch ∙ עברית ∙ हिन्दी ∙ hrvatski ∙ magyar ∙ հայերեն ∙ Bahasa Indonesia ∙ Ido ∙ italiano ∙ 日本語 ∙ ქართული ∙ 한국어 ∙ lietuvių ∙ македонски ∙ മലയാളം ∙ Bahasa Melayu ∙ မြန်မာဘာသာ ∙ norsk bokmål ∙ Plattdüütsch ∙ Nederlands ∙ norsk nynorsk ∙ norsk ∙ occitan ∙ polski ∙ prūsiskan ∙ português ∙ português do Brasil ∙ română ∙ русский ∙ sicilianu ∙ Scots ∙ slovenčina ∙ slovenščina ∙ српски / srpski ∙ svenska ∙ தமிழ் ∙ ไทย ∙ Türkçe ∙ татарча / tatarça ∙ українська ∙ vèneto ∙ Tiếng Việt ∙ 中文 ∙ 中文（中国大陆） ∙ 中文（简体） ∙ 中文（繁體） ∙ 中文（马来西亚） ∙ 中文（新加坡） ∙ 中文（臺灣） ∙ +/−

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