DescripcióEisenstein primes.png	Eisenstein primes with norm less than 500
Data	6 de desembre de 2007
Font	own work, with Mathematica 6.0
Autor	Hack
Permís (Com reutilitzar aquest fitxer)	Public Domain
Altres versions	File:Eisenstein primes.svg és una versió vectorial (SVG) d'aquest fitxer. En cas de ser millor, hauria de ser emprada en lloc d'aquesta imatge tramada. File:Eisenstein primes.png → File:Eisenstein primes.svg Per a més informació pel que fa als gràfics vectorials, llegiu la transició a SVG en Commons. També hi ha informació quant a la compatibilitat del MediaWiki amb les imatges SVG. En altres idiomes 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 ∙ українська ∙ ქართული ∙ հայերեն ∙ বাংলা ∙ தமிழ் ∙ മലയാളം ∙ ไทย ∙ 한국어 ∙ 日本語 ∙ 简体中文 ∙ 繁體中文 ∙ עברית ∙ العربية ∙ فارسی ∙ +/−

Jo, el titular del copyright d'aquesta obra, l'allibero al domini públic. Això s'aplica a tot el món. En alguns països això pot no ser legalment possible, en tal cas: Jo faig concessió a tothom del dret d'usar aquesta obra per a qualsevol propòsit, sense cap condició llevat d'aquelles requerides per la llei.

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