超幾何定理を利用した円周率の公式
twitterで何度も擦ってる、超幾何級数から生まれた円周率の公式です。既出だったら教えてください。
目次
- 公式の内容
- 超幾何定理
- 証明１
- 証明２
公式の内容
$$ \pi=3+\sum_{n=1}^{\infty}\frac{(2n-1)!!}{2^{n+1}(n+\frac{1}{2})(n+1)!} $$
この式を示すために、まず以下の式を証明します
$$ \frac{\pi}{\sin {\pi x}}=\frac{1}{\Gamma(-x)}\sum_{n=0}^{\infty}\frac{\Gamma(-x+n)}{n!(x-n)} \\ただし x>-1,x\notin \mathbb{Z} $$
まず、超幾何級数について考えます。
超幾何定理
超幾何級数は、以下の形で表される多変数関数です。
$$ F(a,b,c;z)=\sum_{n=0}^{\infty}\frac{(a)_n(b)_n}{(c)_nn!}z^n \\(a)_n=a(a+1)(a+2)...(a+n-1)=\frac{\Gamma(a+n)}{\Gamma(a)} $$
$$ F(a,b,c;1)=\frac{\Gamma(c)\Gamma(c-a-b )}{\Gamma(c-a)\Gamma(c-b)} $$
超幾何級数についての説明は このサイト に載ってますので、ここではざっくり証明します。
超幾何級数は
$$
F(a,b,c;z)=\frac{\Gamma(c)}{\Gamma(a)\Gamma(c-a)}\int_{0}^{1}t^{a-1}(1-t)^{c-a-1}(1-tz)^{-b}dt
$$
というEuler積分で表されます。ここで、z=1を代入すると
$$
F(a,b,c;1)\\=\frac{\Gamma(c)}{\Gamma(a)\Gamma(c-a)}\int_{0}^{1}t^{a-1}(1-t)^{c-a-1}(1-t)^{-b}dt
\\=\frac{\Gamma(c)}{\Gamma(a)\Gamma(c-a)}\int_{0}^{1}t^{a-1}(1-t)^{c-a-b-1}dt
$$
実はこの積分、ベータ函数といって、ガンマ関数で表せちゃうんです！
$$
B(x,y)=\int_{0}^{1}t^{x-1}(1-t)^{y-1}dt=\frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}
$$
これをさっきの式に代入すると...
$$
\frac{\Gamma(c)}{\Gamma(a)\Gamma(c-a)} \cdot B(a,c-a-b)
\\=\frac{\Gamma(c)}{\Gamma(a)\Gamma(c-a)} \cdot \frac{\Gamma(a)\Gamma(c-a-b)}{\Gamma(c-b)}
\\
\therefore F(a,b,c;1)=\frac{\Gamma(c)\Gamma(c-a-b
)}{\Gamma(c-a)\Gamma(c-b)}
$$
これで、超幾何定理が示せました。
では、いよいよ公式の証明をしましょう。
証明１
まず、超幾何級数にa=b=-x,c=1-x,z=1を代入
すると
$$
F(-x,-x,1-x;1)=\sum_{n=0}^{\infty}\frac{(-x)_n^2}{(1-x)_nn!}
\\\frac{(-x)_n^2}{(1-x)_nn!}=\frac{1}{n!}\cdot\frac{\Gamma(1-x)}{\Gamma(1-x+n)}\cdot\frac{\Gamma(-x+n)^2}{\Gamma(-x)^2}
\\ \qquad\qquad\ \ \ =\frac{-x\Gamma(-x+n)}{n!(-x+n)\Gamma(-x)}
\\F(-x,-x,1-x;1)=-\frac{x}{\Gamma(-x)}\sum_{n=0}^{\infty}\frac{\Gamma(-x+n)}{n!(-x+n)} ...①
$$
また、超幾何定理より
$$
F(-x,-x,1-x;1)=\frac{\Gamma(1-x)\Gamma(1+x
)}{\Gamma(1)\Gamma(1)}=x\Gamma(x)\Gamma(1-x
)...②
$$
ここで、相反公式を使います。
$$ \Gamma(x)\Gamma(1-x)=\frac{\pi}{\sin \pi x}\qquad(x\notin \mathbb{Z} ) $$
すると②は
$$
F(-x,-x,1-x;1)=\frac{\pi x}{\sin \pi x}...③
$$
①,③より
$$
\frac{\pi x}{\sin \pi x}=\frac{x}{\Gamma(-x)}\sum_{n=0}^{\infty}\frac{\Gamma(-x+n)}{n!(x-n)}
\\\therefore \frac{\pi}{\sin \pi x}=\frac{1}{\Gamma(-x)}\sum_{n=0}^{\infty}\frac{\Gamma(-x+n)}{n!(x-n)}
$$
今回は抵触しないので、範囲については割愛させていただきます。
ここまできたら、あとは簡単なお仕事です。
証明２
先ほどの式にx=1/2を代入すると...
$$
\pi=\frac{1}{\Gamma(-\frac{1}{2})}\sum_{n=0}^{\infty}\frac{\Gamma(-\frac{1}{2}+n)}{n!(\frac{1}{2}-n)}
$$
ガンマ関数の特殊値
$$
\Gamma\left( \frac{1}{2} \right)=\sqrt{\pi}
\\\Gamma\left( \frac{1}{2}+n \right)=\frac{(2n-1)!!}{2^n}\sqrt{\pi}
\\\Gamma\left( \frac{1}{2}-n \right)=\frac{(-2)^n}{(2n-1)!!}\sqrt{\pi}
$$
より、
$$
\pi=\frac{1}{\Gamma(-\frac{1}{2})}\sum_{n=0}^{\infty}\frac{\Gamma(-\frac{1}{2}+n)}{n!(\frac{1}{2}-n)}
\\ \ \ =\frac{1}{2\sqrt{\pi}}\sum_{n=0}^{\infty}\frac{\Gamma(\frac{1}{2}+n-1)}{n!(n-\frac{1}{2})}
\\ \ \ =\frac{1}{2\sqrt{\pi}} \left(4\sqrt{\pi}+2\sqrt{\pi}+\sum_{n=2}^{\infty}\frac{1}{n!(n-\frac{1}{2})}\cdot \frac{(2n-3)!!\sqrt{\pi}}{2^{n-1}}\right)
\\ \ \ =3+\frac{1}{2}\sum_{n=2}^{\infty}\frac{1}{n!(n-\frac{1}{2})}\cdot \frac{(2n-3)!!}{2^{n-1}}
$$
最後に、形をきれいにすると
$$
\pi=3+\sum_{n=1}^{\infty}\frac{(2n-1)!!}{2^{n+1}(n+\frac{1}{2})(n+1)!}
$$
これで円周率の式が導けました！！！
ほかにも、公式２に3/2を代入すれば
$$
\pi= \frac{37}{12}+\frac{3}{4}\sum_{n=1}^{\infty}\frac{(2n-1)!!}{2^n(n+\frac{1}{2})(n+2)!}
$$
という式が得られます。
(本当はこれのはずだったのに勘違いしてたなんて口が裂けても言えない)
こっちのほうが収束が速いので、もしかしたら入れる値を大きくすれば収束の速さが速くなるかもしれません。わかる方がいたら教えていただければ幸いです。
読んでいただきありがとうございました！！
