いぬおさんのおもしろ数学実験室

おいしい紅茶でも飲みながら数学、物理、工学、プログラミング、そして読書を楽しみましょう

ライプニッツの公式で円周率を計算できるか?

 次のような級数を考えましょう。

f:id:Inuosann:20200306213425p:plain

この級数は収束することを示してみます。証明の流れが分かればよしとします。一般的な証明はこれから書くことが理解できれば難しくありません。

第n部分和を次のようにおきます。

f:id:Inuosann:20200306213531p:plain

偶数項目までの和を考えます。例えば、

f:id:Inuosann:20200306213717p:plain

ですし、一方では

f:id:Inuosann:20200306213931p:plain

ですから、

f:id:Inuosann:20200306214209p:plain

f:id:Inuosann:20200306214954p:plain

f:id:Inuosann:20200306223922p:plain

実は、不思議なことにこの級数の値の4倍がπに等しいのです。ライプニッツの公式と呼ぶらしいです(ウィキペディア)。微積分の知識を用いると明らかになるのですが、ここではこの事実を認めておきましょう。さて、

f:id:Inuosann:20200306215133p:plain

f:id:Inuosann:20200306215204p:plain

 等が成り立ちますから、

f:id:Inuosann:20200306215634p:plain

をn項で打ち切ってπを求めたときの誤差の絶対値は最大でも

f:id:Inuosann:20200306224109p:plain

以下です。ライプニッツの公式は不思議だしきれいだし素晴らしいのですが、実際にπを求めるとなるとあまりうまくありません。例えば誤差を0.001以下に抑えようとすると(小数第3位まで求めようとすると)

f:id:Inuosann:20200306224148p:plain

が成立しなければなりません。これを満たす最小のnはおよそ2000。つまり2000項まで計算しないと誤差が0.001以下にならないということです。ライプニッツの公式はπを求めるのには使いにくいのです。たくさんの桁数を計算したければ、もっと収束の速い公式(マチンの公式など)があり、そちらを使います。

 

  前、かなり丁寧に書きました。こちらもご参考に。誤差の評価はしていませんが。

www.omoshiro-suugaku.com