ライプニッツの公式

ライプニッツの公式（ライプニッツのこうしき、テンプレート:Lang-en）とは円周率の値を求めるための公式の一つである。以下の級数で表される。

$1-\frac{1}{3}+\frac{1}{5}-\frac{1}{7}+\frac{1}{9}-\cdots =\frac{\pi }{4}$

これは初項が 1 で各項が奇数の逆数である交項級数が テンプレート:Math に収束することを意味する。総和の記号を用いると以下のようになる。

$\sum _{n=0}^{\mathrm{\infty }}\frac{(-1{)}^n}{2n+1}=\frac{\pi }{4}$

この公式を名付けたのはライプニッツであるが、これはすでに15世紀のインドの数学者マーダヴァがライプニッツより300年ほど前に発見していたものである。公式の発見がマーダヴァの功績であることを示すためにマーダヴァ-ライプニッツ級数と呼ばれることもある。

証明

冪級数展開を用いる証明

三角関数の一つ テンプレート:Math を テンプレート:Mvar について微分すると

$\frac{\mathrm{d}}{\mathrm{d}\theta }\tan \theta =1+{\tan }^2\theta$

となる。ここで テンプレート:Math とおくと

$\frac{\mathrm{d}x}{\mathrm{d}\theta }=1+x^2,\frac{\mathrm{d}\theta }{\mathrm{d}x}=\frac{1}{1+x^2}\cdots (1)$

が導かれる。

また以下の等比級数を考える。

$1-x^2+x^4-x^6+x^8-\cdots =\frac{1}{1+x^2}(|x|<1)\cdots (2)$

左辺は公比が テンプレート:Math であり、テンプレート:Math すなわち テンプレート:Math のとき テンプレート:Math に収束する。(1), (2)式から

$\frac{\mathrm{d}\theta }{\mathrm{d}x}=1-x^2+x^4-x^6+x^8-\cdots (|x|<1)$

が得られる。この両辺を テンプレート:Mvar について項別積分すれば

$\theta =x-\frac{x^3}{3}+\frac{x^5}{5}-\frac{x^7}{7}+\frac{x^9}{9}-\cdots (|x|<1)\cdots (3)$

となる（この時、左辺をarctan xと表すとグレゴリー級数のかたちとなる）。（テンプレート:Mathのとき テンプレート:Math であるから定数項は テンプレート:Math である。）テンプレート:Math としたので テンプレート:Math のとき テンプレート:Math である。これを利用して(3)式に テンプレート:Math と テンプレート:Math を代入すると

$\frac{\pi }{4}=1-\frac{1}{3}+\frac{1}{5}-\frac{1}{7}+\frac{1}{9}-\cdots$

という式が現れる。ただし テンプレート:Math は テンプレート:Mathの条件に反するので(3)式に テンプレート:Math を代入できるかどうかが問題になるが、この場合は代入してもよいことが分かっている（アーベルの連続性定理）。

フーリエ級数を用いた証明

方形波をフーリエ級数で表す証明法もある。方形波 テンプレート:Math を

$f(x)=\{\begin{array}{cc}-1& -\pi \le x<0\\ 1& 0\le x<\pi \end{array}$

と定義する（これは区分的に滑らかな関数で テンプレート:Math 上可積分である）と、フーリエ係数 テンプレート:Math はこの方形波が奇関数なので テンプレート:Math であり、テンプレート:Math は以下の式で表す。

$\begin{array}{cc}{b}_n& =\frac{1}{\pi }\underset{-\pi }{\overset{\pi }{\int }}f(x)\sin nxdx\\ & =\frac{1}{\pi }\{\underset{-\pi }{\overset{0}{\int }}(-1)\sin nxdx+\underset{0}{\overset{\pi }{\int }}1\cdot \sin nxdx\}\\ \end{array}$

これを計算すると以下のようになる。

$b_n=\{\begin{array}{cc}\frac{4}{n\pi }& n=2k+1(k\in \mathbb{N})\\ 0& n=2k\end{array}$

したがって方形波のフーリエ級数は

$f(x)=\frac{4}{\pi }(\sin x+\frac{1}{3}\sin 3x+\frac{1}{5}\sin 5x+\cdots )$

となり、テンプレート:Math は テンプレート:Math において連続であるから、両辺に テンプレート:Math を代入すると

$1=\frac{4}{\pi }(1-\frac{1}{3}+\frac{1}{5}-\cdots )$

であるのでライプニッツの公式が得られる。

性質

この公式は単純な形をしているが、実際の円周率の計算に用いるには収束が非常に遅いために全く適していない。10進法での正確な値 (= 3.1415926535…) を10桁分計算するだけで100億回以上の計算を要するほどである。ちなみに最初の500万項の部分和を計算すると テンプレート:Mvar の近似値として

3.1415924535897932384646433832795027841971693993873058…

が得られる。下線の引かれている桁だけ間違っているが、こういった誤差がいくらになるのか予想することは次の近似式で可能である。

$\frac{\pi }{4}-\sum _{n=0}^{\frac{N}{2}}\frac{(-1{)}^n}{2n+1}\sim \frac{1}{2}\sum _{m=0}^{\mathrm{\infty }}\frac{E_{2m}}{N^{2m+1}}$

テンプレート:Math はオイラー数、テンプレート:Mvar は 4 で割り切れる自然数である。テンプレート:Mvar に 10 の累乗数を代入すると、右辺の部分和からこの公式で求めた10進法表記での誤りが現れる桁の位置とその誤差を計算できる。

関連項目

• 円周率の歴史
• マチンの公式
• ゴットフリート・ライプニッツ