14
大学数学基礎解説
文献あり

超幾何定理を利用した円周率の公式

666
0
$$$$

twitterで何度も擦ってる、超幾何級数から生まれた円周率の公式です。既出だったら教えてください。

目次

  1. 公式の内容
  2. 超幾何定理
  3. 証明1
  4. 証明2

公式の内容

円周率

$$ \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)} $$

これで、超幾何定理が示せました。
では、いよいよ公式の証明をしましょう。

証明1

まず、超幾何級数に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)} $$
今回は抵触しないので、範囲については割愛させていただきます。
ここまできたら、あとは簡単なお仕事です。

証明2

先ほどの式に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)!} $$
これで円周率の式が導けました!!!
ほかにも、公式2に3/2を代入すれば
$$ \pi= \frac{37}{12}+\frac{3}{4}\sum_{n=1}^{\infty}\frac{(2n-1)!!}{2^n(n+\frac{1}{2})(n+2)!} $$
という式が得られます。
(本当はこれのはずだったのに勘違いしてたなんて口が裂けても言えない)
こっちのほうが収束が速いので、もしかしたら入れる値を大きくすれば収束の速さが速くなるかもしれません。わかる方がいたら教えていただければ幸いです。

読んでいただきありがとうございました!!

参考文献

投稿日:202346
更新日:20231112

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

かくことは あれども筆は 進まぬか コマンド打つのが 割とめんどい

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中