[3,3,3,4007500000] dual polyhedron
“辺寸が 10mm で 地球赤道周と同じ 正多角形での 反角柱 antiprism の
双対多面体の長さはいくらか？” という問題について お伝えします。
私は 15桁を越す数値計算には 多倍長電卓LM というフリーウエアーを用いています。
( 参照 URL=http://www.vector.co.jp/soft/win95/personal/se242555.html )
以下がそれで作ったプログラムです。C言語的なソフトです。
//--------------コピー開始------------------ a1=40075*1000*100; // n角形を指定 約40,075 km b1= pi*2/a1; // 360/角数 c1=0.5/sin(b1/2); // 外接円柱半径 d1=0.5*tan(b1/4); // n角の辺・心寸 e1=sqrt(3)/2; // 3角かど・辺寸 f1=sqrt(e1^2-d1^2); // 3角かど・辺寸 軸面投影 g1=sqrt((f1/2)^2+c1^2); // 外接球半径 h1=0.5/g1; // 角錐高 i1=asin(h1); // 稜仰角 j1=g1*cos(i1); // 稜芯寸 k1=pi*2/a1/2; // 双n接合角/2 l1=0.5/tan(k1); // 双n辺心寸 m1=asin(l1/j1); // 双n仰角 n1=j1/cos(m1); // n頂芯寸 print ""; print "[3,3,3,4007500000] dual polyhedron"; print ""; print "n頂芯寸"; print n1; print ""; print "光年"; // 9 460 730 472 580 800 m print n1*2/(9460730472580800*100); //------------コピー終わり------------------以下が 20桁指定での計算結果です
[3,3,3,4007500000] dual polyhedron n頂芯寸 = 939478161669236009.35 光年 = 1.9860584008645899692
2013年7月28日
