描述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 ∙ العربية ∙ беларуская (тарашкевіца) ∙ български ∙ বাংলা ∙ 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 ∙ 中文 ∙ 中文（中国大陆） ∙ 中文（简体） ∙ 中文（繁體） ∙ 中文（马来西亚） ∙ 中文（新加坡） ∙ 中文（臺灣） ∙ 中文（臺灣） ∙ +/−

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

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