描述Eisenstein primes.png	Eisenstein primes with norm less than 500
日期	2007年12月6日
來源	own work, with Mathematica 6.0
作者	Hack
授權許可 (重用此檔案)	Public Domain
其他版本	File:Eisenstein primes.svg是本檔案的向量版本。 如果品質不低，就應該優先使用該檔案，而非PNG檔案。 File:Eisenstein primes.png → File:Eisenstein primes.svg 更多資訊請參閱Help:SVG/zh。 其他語言 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 ∙ українська ∙ ქართული ∙ հայերեն ∙ বাংলা ∙ தமிழ் ∙ മലയാളം ∙ ไทย ∙ 한국어 ∙ 日本語 ∙ 简体中文 ∙ 繁體中文 ∙ עברית ∙ العربية ∙ فارسی ∙ +/−

我，此作品的版權所有人，釋出此作品至公共領域。此授權條款在全世界均適用。 這可能在某些國家不合法，如果是的話： 我授予任何人有權利使用此作品於任何用途，除受法律約束外，不受任何限制。

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