Mathematicaで算額奉納事始め
Mathematicaで算額奉納事始めはじめに 年賀状としてMathematicaで算額を記述した。年賀状なので、問題は平易な物にし、また画題は簡易な物にした。正十二面体の頂点座標は黄金比ϕを使って描画する。Mathematicaでそれら頂点座標と面リストを入力して正多面体を描画した。ネットでは正十二面体の黄金比の頂点座標を±記号を使って、数個の頂点座標を示してこれが正十二面体の頂点座標ですと云った例が見掛けるが、これは正しくない。正十二面体の頂点座標は20個有り、それに対応する面リストが必須で、これが無いと正十二面体は正確に描画できない。Mathematicaでは黄金比をϕと記す。正多面体は五個しかない、正十二面体と正二十面体には黄金比ϕ が潜んでいる。正多面体は、数学の宝の採石場で、玉を掘り当てるのは数学マニアの醍醐味である。古来より正多面体は紙とペンで研究されてきた。一昔前はHemmi計算尺、Tiger計算機から電卓、コンピュータが一般的になると、ソフトはBasic ,Fortan、Algol、Cが使われた。今は豊富な正多面体関数を包含しているMathematicaが簡易だ。シリーズ物 Mathematicaで黄金比の事始め Mathematicaで正多面体の黄金比の三角関数の表 Mathematicaで正多面体の三角関数を解くϕ= (1 + Sqrt[5])/2N[1/ϕ] = 0.6180339887498948N[ϕ] = 1.618033988749895N[ϕ2] = 2.6180339887498951.問題:二つの立体の間の体積を求める。2.略解2.1.正十二面体の体積は (15+7Sqrt[5] )a3/42.2.正六面体の体積はb33.解:一辺の長さをa= 2 にすると、正十二面体の体積は (15+7Sqrt[5] )a3/4=7.6631a3になる。一辺の長さをb= 2ϕ=(1+Sqrt[5])=3.2361 にすると、正六面体(Cube)の体積はb3になる。よって二つの立体の間の体積は7.6631a3-b3=7.6631*23-3.23613=61.3048-33.8895=27.4153になる。4.正多面体の頂点座標4.1.正十二面体正十二面体の頂点座標は黄金比ϕを使う。黄金比ϕ2,ϕ,1,0を使う。正十二面体の頂点座標は{{0, -1, -ϕ2}, {0, +1, -ϕ2}, {0, -1, +ϕ2}, {0, +1, +ϕ2}, {-1, -ϕ2, 0}, {+1, -ϕ2,0}, {-1, +ϕ2, 0}, {+1, +ϕ2,0}, {-ϕ2, 0, -1}, {-ϕ2, 0, +1},{+ϕ2, 0, -1}, {+ϕ2, 0, +1}, {-ϕ, -ϕ, -ϕ}, {-ϕ, -ϕ, +ϕ}, {-ϕ, +ϕ, -ϕ}, {-ϕ, +ϕ, +ϕ]}, {+ϕ, -ϕ, -ϕ}, {+ϕ, -ϕ, +ϕ}, {+ϕ, +ϕ, -ϕ}, {+ϕ, +ϕ, +ϕ}}になる。この場合辺の長さは2になる。RegularDodecahedronVertexData ={{0, -1, -ϕ2}, {0, +1, -ϕ2}, {0, -1, +ϕ2}, {0, +1, +ϕ2}, {-1, -ϕ2, 0}, {+1, -ϕ2,0}, {-1, +ϕ2, 0}, {+1, +ϕ2,0}, {-ϕ2, 0, -1}, {-ϕ2, 0, +1},{+ϕ2, 0, -1}, {+ϕ2, 0, +1}, {-ϕ, -ϕ, -ϕ}, {-ϕ, -ϕ, +ϕ}, {-ϕ, +ϕ, -ϕ}, {-ϕ, +ϕ, +ϕ]}, {+ϕ, -ϕ, -ϕ}, {+ϕ, -ϕ, +ϕ}, {+ϕ, +ϕ, -ϕ}, {+ϕ, +ϕ, +ϕ}}In[]:RegularDodecahedronVertexIndex ={{1, 2, 19, 11, 17}, {2, 1, 13, 9, 15}, {3, 4, 16, 10, 14}, {4, 3,18, 12, 20}, {5, 6, 18, 3, 14}, {6, 5, 13, 1, 17}, {7, 8, 19, 2,15}, {8, 7, 16, 4, 20}, {9, 10, 16, 7, 15}, {10, 9, 13, 5,14}, {11, 12, 18, 6, 17}, {12, 11, 19, 8, 20}}dp2 = RegularDodecahedron =Polyhedron[RegularDodecahedronVertexData,RegularDodecahedronVertexIndex]Graphics3D[{EdgeForm[{Thin, Blue}], FaceForm[{Pink, Opacity[0.7]}],dp2}, Boxed -> False]出力図は後記。4.2.正六面体ϕ= (1 + Sqrt[5])/22.1.黄金比ϕを使う。正六面体の頂点座標は{{-ϕ, -ϕ, -ϕ}, {-ϕ,-ϕ, ϕ}, {-ϕ, +ϕ, -ϕ}, {-ϕ, +ϕ, +ϕ},{+ϕ, -ϕ, -ϕ}, {+ϕ, -ϕ, ϕ}, {+ϕ, +ϕ, -ϕ}, {+ϕ, +ϕ, +ϕ}}になる。この場合辺の長さは2ϕになる。In[]:RegularCubeVertexData ={{-ϕ, -ϕ, -ϕ}, {-ϕ,-ϕ, ϕ}, {-ϕ, +ϕ, -ϕ}, {-ϕ, +ϕ, +ϕ},{+ϕ, -ϕ, -ϕ}, {+ϕ, -ϕ, ϕ}, {+ϕ, +ϕ, -ϕ}, {+ϕ, +ϕ, +ϕ}}RegularCubeVertexIndex = {{1, 2, 4, 3}, {1, 3, 7, 5}, {1, 5, 6, 2}, {8, 6, 5, 7}, {8, 4, 2,6}, {8, 7, 3, 4}}cp3 = RegularCube =Polyhedron[RegularCubeVertexData, RegularCubeVertexIndex]Graphics3D[{EdgeForm[{Thin, Blue}], FaceForm[{Pink, Opacity[0.7]}],cp3}, Boxed -> False]出力図は後記。これら各正多面体の頂点座標を黄金比ϕで記述した描画の副産物として、正十二面体の中に正六面体が収まる。正十二面体黄金比頂点座標ϕ2,ϕ,1,0→正六面体黄金頂点座標比ϕ。In[]:Graphics3D[{EdgeForm[{Thin, Blue}], FaceForm[{Orange, Opacity[0.3]}],cp3, FaceForm[{Pink, Opacity[0.2]}], dp2}, Boxed -> False,ViewPoint -> {10000, 10000, 10000}]正十二面体の中に正六面体が丁度収まる。頂点が飛び出たり、立体が内部でころころしない。Out[]:番外編デジタル正月飾り 植物を切り取ったり、痛めつける事がない、究極の自然に優しい正月飾りである。年末年始に咲く福寿草で目出度い黄金、木瓜と果実が付く南天を使って紅白を表し、正月飾りを仕立てた。添付品として3Dモデルの扇子と青竹を使用した。青葉、青竹は新しい成長を表す。古代エジプトでナイルの氾濫後、新しい草が芽吹く。そこにエジプト文明が開化した。キリスト教はこのプロセスを自己流に利用し、全てが流された後、そこから新しい芽を吹くという事実を、死と復活の教義に結び付けた。完