この計算ツールでできること
このツールは、急速に収束する3つの有名な無限級数のいずれかを項ごとに足し合わせて、円周率πを計算します。選べるのは、ラマヌジャンが1914年に発表した第1公式、同じく1914年の第2公式、そしてチュドノフスキー兄弟による1987年の公式です。これらの級数が驚異的なのは、1項加えるごとに正しい小数桁が一気に何桁も増える点で、わずか数項だけで倍精度の限界までπを再現できます。これは純粋数学であり、地域や国による規則の違いは一切なく、どこでも同じように使えます。
使い方
プルダウンから公式を選び、足し合わせる項数の上限と、表示する小数桁数を指定します。計算は \(n = 0, 1, 2, \ldots\) と項を加えていき、πの値がそれ以上変化しなくなった時点で自動的に打ち切ります。多くの場合、数項のうちに収束します。内部計算にはIEEE-754の倍精度浮動小数点を用いているため、表示桁数の設定にかかわらず、信頼できるのは有効数字でおよそ15〜16桁です。
公式の仕組み
ラマヌジャン第1公式は、πの逆数を組み立てる形になっています。定数の係数 \(\sqrt{8}/9801\) が無限和に掛かり、その第n項は階乗の比 \((4n)!/(4^n n!)^4\) と、一次の因子 \((1103 + 26390n)\) を \(99^{4n}\) で割ったものを組み合わせています。 $$\frac{1}{\pi} = \frac{\sqrt{8}}{9801} \sum_{n=0}^{\infty} \frac{(4n)!}{(4^n n!)^4} \cdot \frac{1103 + 26390n}{99^{4n}}$$ 和Sが求まれば、πは \(1/(\text{係数} \times S)\) として復元できます。チュドノフスキー公式も同様の仕組みですが収束はさらに速く、1項あたりおよそ14桁ずつ正しい桁が増えます。
計算例
ラマヌジャン第1公式で \(n=0\) の項だけを使ってみましょう。係数は \(\sqrt{8}/9801 = 0.000288583\ldots\)、\(n=0\) の項は \(1 \times 1103 = 1103\) です。よって $$\frac{1}{\pi} = 0.000288583 \times 1103 = 0.31831\ldots$$ となり、\(\pi = 3.14159273\)、すでに小数第6位ほどまで正しい値が得られます。さらに \(n=1\) の項を加えると \(\pi = 3.14159265358979\) となり、約16桁まで正確になります。
よくある質問
表示桁数を増やしても精度が上がらないのはなぜ？ 倍精度浮動小数点が保持できるのは有効数字でおよそ15〜16桁までです。それ以上の精度を得るには任意精度演算が必要になります。
実際に必要な項数はどのくらい？ 倍精度の限界まで求めるなら、ラマヌジャン第1公式で約2項、チュドノフスキー公式ではわずか1〜2項で十分です。
どの級数が最も速い？ 収束が最も速いのはチュドノフスキー公式で、現代のπの桁数記録更新にも実際に使われているアルゴリズムです。
