ジアン角単位系を用いて,円周率πの近似値1億6千万桁の計算することに成功した.計算式履いたってシンプルなもので,
π ~ Im( (-1)^ε )) / ε for extreamly small ε.(ε はジアン角単位系の極めて小さい角度)
を単純に算術演算しているだけだ.Julia には BigFloat という任意長精度の浮動小数点型あるので,実装は容易だ.計算の所要時間は2180秒,36分で.これが速いのか遅いのかはよく分からないが,EXEならもっと速くなることは確かなので,やってみることにしよう.
