正65537角形

1の65537乗根を根号だけで求めてみました。いや、式は一ヶ月以上前に計算で求めていたんですが、検証に時間がかかりました。

というのも通常のIntelのCPUのdoubleでは浮動小数点の精度が悪すぎて、計算できないんですよ。そして前回検証に利用したcalcは変数が多すぎるとコケル・・というわけで、試行錯誤の末128bitの精度をもって計算したら無事計算できました。

計算式（2105行）→結果

答え: 0.9999999954042475651948709986627380845520+0.00009587233620022679702185665995171896365861i
→検証

作図は約2105回もの根号を求めないといけないわけで（もうちょっと減るとは思うけど）、実際にやるのは無理っぽい。128bit級の作図精度も必要ですし・・
Posted by issei

カテゴリ: 数学

関連

コメント一覧

MPはいかが? GNU MP. あるいはlispとかで筆算させるとか。
kaz
Kazさん> GNU MPは算術関数系に不満があって、ちょっと使えず、MPFRって言うのを使いました。こちらは一通り揃って、どう丸めるかなどまで指定できて、かなり良い感じ。整数演算ならGNU MPでも（・∀・）ｲｲ!!と思うんですが。
Isami
（￣∧￣（＿ ＿）ﾌﾑﾌﾑ・・・・・なるほど・・・・ ・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・。 ・・・・・・ごめんなさい。嘘つきました。（_　_。）笑
ゴンゾー