円周率が22/7より小さいことの証明

テンプレート:円周率 有名な数学的事実であるところの、円周率 テンプレート:Math が テンプレート:Math より小さいことの証明（えんしゅうりつが テンプレート:Math ぶんの テンプレート:Math よりちいさいことのしょうめい）は、古代ギリシアのアルキメデスに始まり、何通りも与えられている。本項では、そのうちの一つで、微分積分学の初等的なテクニックのみを用いる、近年に発見された証明を扱う。この証明は、その数学的な美およびディオファントス近似の理論との関係によって、現代数学においても注目されてきた。スティーヴン・ルーカスは、これを「テンプレート:Math の近似に関する最も美しい結果の一つ」と呼び、ジュリアン・ハヴィルは、円周率の連分数近似の議論を終える際に「この結果に言及せざるを得ない」と述べた上で証明を示している。

もし円周率が テンプレート:Math に近いことを知っていれば、テンプレート:Math（テンプレート:Math に近い）よりも小さいことは自明である。しかし、テンプレート:Math を示すのは、テンプレート:Math が テンプレート:Math に近いことを示すよりもずっと手間は小さい。この証明の評価方法は一般化され、円周率の値を計算する系統的な方法になっている。

背景

テンプレート:Math は、テンプレート:Math の正則な連分数展開から得られる近似値の一つであり、また、表現が簡潔であることから テンプレート:Math の近似値として広く用いられている。テンプレート:Math が テンプレート:Math よりも大きいことは、これらの値の十進法での小数展開

$\begin{array}{cc}\frac{22}{7}& =3.\overline{142857}\\ \pi & =3.14159265\cdots \end{array}$

より分かる。

この近似値は古代より知られており、アルキメデスは紀元前3世紀に、テンプレート:Math が円周を大きめに見積もった近似値であることを証明したテンプレート:Refnest。彼の証明は、円に外接する正九十六角形の周長の、円の直径に対する比が、テンプレート:Math よりも小さいことを示すものであった。これ以前の証明は知られていないが、アルキメデス以前から、テンプレート:Math は円周率の近似値として用いられていた形跡がある。

より精密な テンプレート:Math の近似値に、テンプレート:Math がある。これも テンプレート:Math の連分数展開から得られる近似値であって、円周率よりわずかに大きい。

証明

証明の概略は非常に簡潔に述べられる。

$0<\underset{0}{\overset{1}{\int }}\frac{x^4(1-x{)}^4}{1+x^2}\mathrm{d}x=\frac{22}{7}-\pi$

より、テンプレート:Math は テンプレート:Math よりも小さい。

この積分の計算は、1968年のテンプレート:仮リンクの最初の問題として出題された。 パトナムの競技会の他のほとんどの問題は、この問題よりも難しい。しかし、この競技会は、直接にはあまり知られていなくとも、実は非常によく知られたことに言及していると判明する問題を、しばしば呼び物にしている。この積分は、インド工科大学の入試問題に用いられたこともある。

積分の評価の詳細

被積分関数の分母と分子が共に非負であることから、この積分の値が正であることは直ちに分かる。よって、あとはこの積分の値を、有理関数の標準的な積分の手順に従って求めればよい。

$\begin{array}{cc}0& <\underset{0}{\overset{1}{\int }}\frac{x^4(1-x{)}^4}{1+x^2}\mathrm{d}x\\ & \phantom{<}=\underset{0}{\overset{1}{\int }}\frac{x^4-4x^5+6x^6-4x^7+x^8}{1+x^2}\mathrm{d}x\\ & \phantom{<}=\underset{0}{\overset{1}{\int }}(x^6-4x^5+5x^4-4x^2+4-\frac{4}{1+x^2})\mathrm{d}x\\ & \phantom{<}=[\frac{x^7}{7}-\frac{2x^6}{3}+x^5-\frac{4x^3}{3}+4x-4\arctan x{]}_0^1\\ & \phantom{<}=\frac{1}{7}-\frac{2}{3}+1-\frac{4}{3}+4-\pi \\ & \phantom{<}=\frac{22}{7}-\pi \end{array}$

直ちに得られる円周率の評価

Dalzell は1944年の論文で、この積分に言及し、被積分関数の分母に x = 1 を代入すると下からの評価が、x = 0 を代入すると上からの評価が得られると指摘した。すなわち、

$\frac{1}{1260}=\underset{0}{\overset{1}{\int }}\frac{x^4(1-x{)}^4}{2}dx<\underset{0}{\overset{1}{\int }}\frac{x^4(1-x{)}^4}{1+x^2}dx<\underset{0}{\overset{1}{\int }}\frac{x^4(1-x{)}^4}{1}dx=\frac{1}{630}$

であり、これより円周率 テンプレート:Π の評価として

$\frac{22}{7}-\frac{1}{630}<\pi <\frac{22}{7}-\frac{1}{1260}$

を得る。十進小数で表現すると 3.1412… < テンプレート:Π < 3.1420… となる。これらの評価値の誤差は、0.015% 未満である。

より良い評価

テンプレート:Π のより良い上限としてよく知られた テンプレート:Sfrac は、

$0<\underset{0}{\overset{1}{\int }}\frac{x^8(1-x{)}^8(25+816x^2)}{3164(1+x^2)}dx=\frac{355}{113}-\pi$

として与えられる。この値は十進小数で表すと

$\frac{355}{113}=3.14159292\cdots$

であり、小数点以下6桁まで テンプレート:Π と一致する。被積分関数の分母に x = 1 を代入すると、積分値の下限

$\underset{0}{\overset{1}{\int }}\frac{x^8(1-x{)}^8(25+816x^2)}{6328}dx=\frac{911}{5261111856}=0.000000173\cdots$

を得る。x = 0 を代入すると、積分値の上限として、上記の数の2倍を得るので、円周率の評価として

$\frac{355}{113}-\frac{911}{2630555928}<\pi <\frac{355}{113}-\frac{911}{5261111856}$

を得る。この評価は、十進小数として表すと 3.141 592 57… < テンプレート:Π < 3.141 592 74… であり、太文字で表した小数点以下6桁まで円周率と一致する。

一般化

テンプレート:Harvtxt や テンプレート:Harvtxt で示唆されているように、ここまでの議論は一般化される（ただし、どちらの文献も詳しい計算は含んでいない）。任意の正整数 n に対して

$\frac{1}{2^{2n-1}}\underset{0}{\overset{1}{\int }}{x}^{4n}(1-x{)}^{4n}dx<\frac{1}{2^{2n-2}}\underset{0}{\overset{1}{\int }}\frac{x^{4n}(1-x{)}^{4n}}{1+x^2}dx<\frac{1}{2^{2n-2}}\underset{0}{\overset{1}{\int }}{x}^{4n}(1-x{)}^{4n}dx,$

が成り立ち、第二の積分は

$\begin{array}{cc}& \frac{1}{2^{2n-2}}\underset{0}{\overset{1}{\int }}\frac{x^{4n}(1-x{)}^{4n}}{1+x^2}dx\\ & =\sum _{j=0}^{2n-1}\frac{(-1{)}^j}{2^{2n-j-2}(8n-j-1)(\genfrac{}{}{0ex}{}{8n-j-2}{4n+j})}+(-1{)}^n(\pi -4\sum _{j=0}^{3n-1}\frac{(-1{)}^j}{2j+1})\end{array}$

となって テンプレート:Π を含む式を得る。この式の最後の和は、ライプニッツの公式に現れる級数を有限で切ったものになっている。戻って第一の積分は

$\frac{1}{2^{2n-1}}\underset{0}{\overset{1}{\int }}{x}^{4n}(1-x{)}^{4n}dx=\frac{1}{2^{2n-1}(8n+1)(\genfrac{}{}{0ex}{}{8n}{4n})}$

であり、第三の積分はこの2倍である。この値は、n が大きくなるにつれて急激に 0 に近付くので、ここで議論している方法が、テンプレート:Π の近似値を求めるのに適していることが分かる。

これらの積分の計算

任意の整数 k ≥ 0 と l ≥ 2 に対して、

$\begin{array}{cc}{x}^k(1-x{)}^l& =(1-2x+x^2)x^k(1-x{)}^{l-2}\\ & =(1+x^2)x^k(1-x{)}^{l-2}-2x^{k+1}(1-x{)}^{l-2}\end{array}$

であるから、この変形を 2n 回繰り返すことにより

$x^{4n}(1-x{)}^{4n}=(1+x^2)\sum _{j=0}^{2n-1}(-2{)}^j{x}^{4n+j}(1-x{)}^{4n-2(j+1)}+(-2{)}^{2n}{x}^{6n}.$

を得る。また、

$\begin{array}{cc}{x}^{6n}-(-1{)}^{3n}& =\sum _{j=1}^{3n}(-1{)}^{3n-j}{x}^{2j}-\sum _{j=0}^{3n-1}(-1{)}^{3n-j}{x}^{2j}\\ & =\sum _{j=0}^{3n-1}((-1{)}^{3n-(j+1)}{x}^{2(j+1)}-(-1{)}^{3n-j}{x}^{2j})\\ & =-(1+x^2)\sum _{j=0}^{3n-1}(-1{)}^{3n-j}{x}^{2j},\\ \end{array}$

である。ここに、第一の等式が成り立つのは j = 1, …, 3n − 1 の項が打ち消しあって消えるためであり、第二の等式では、第一の和において j = jテンプレート:Sub + 1 とおいて jテンプレート:Sub を改めて j としている。

これらの結果を用いて

$\begin{array}{cc}\frac{x^{4n}(1-x{)}^{4n}}{2^{2n-2}(1+x^2)}& =\sum _{j=0}^{2n-1}\frac{(-1{)}^j}{2^{2n-j-2}}{x}^{4n+j}(1-x{)}^{4n-2j-2}\\ & -4\sum _{j=0}^{3n-1}(-1{)}^{3n-j}{x}^{2j}+(-1{)}^{3n}\frac{4}{1+x^2}(*)\end{array}$

を得る。

次に、任意の整数 k, l ≥ 0 に対して、部分積分を l 回適用することにより、

$\begin{array}{cc}\underset{0}{\overset{1}{\int }}{x}^k(1-x{)}^{l}dx& =\frac{l}{k+1}\underset{0}{\overset{1}{\int }}{x}^{k+1}(1-x{)}^{l-1}dx\\ & =\cdots \\ & =\frac{l}{k+1}\frac{l-1}{k+2}\cdots \frac{1}{k+l}\underset{0}{\overset{1}{\int }}{x}^{k+l}dx\\ & =\frac{1}{(k+l+1)(\genfrac{}{}{0ex}{}{k+l}{k})}(**)\end{array}$

を得る。特に、k = l = 4n とおけば

$\underset{0}{\overset{1}{\int }}{x}^{4n}(1-x{)}^{4n}dx=\frac{1}{(8n+1)(\genfrac{}{}{0ex}{}{8n}{4n})}$

を得る。

(*) を 0 から 1 まで積分し、(**) と arctan 1 = テンプレート:Sfrac を用いると、テンプレート:Π を含んだ目的の積分値を得る。

例

n = 1 の場合はすでに見た。n = 2 の場合は、

$\frac{1}{4}\underset{0}{\overset{1}{\int }}\frac{x^8(1-x{)}^8}{1+x^2}dx=\pi -\frac{47171}{15015}$

および

$\frac{1}{8}\underset{0}{\overset{1}{\int }}{x}^8(1-x{)}^{8}dx=\frac{1}{1750320}$

より

$\frac{47171}{15015}+\frac{1}{1750320}<\pi <\frac{47171}{15015}+\frac{2}{1750320}$

すなわち 3.141 592 31… < テンプレート:Π < 3.141 592 88… となり、テンプレート:Π の近似値を小数点以下第6位まで得る。

同様に n = 3 の場合は、

$\frac{1}{16}\underset{0}{\overset{1}{\int }}\frac{x^{12}(1-x{)}^{12}}{1+x^2}dx=\frac{431302721}{137287920}-\pi$

および

$\frac{1}{32}\underset{0}{\overset{1}{\int }}{x}^{12}(1-x{)}^{12}dx=\frac{1}{2163324800}$

より

$\frac{431302721}{137287920}-\frac{2}{2163324800}<\pi <\frac{431302721}{137287920}-\frac{1}{2163324800}$

すなわち 3.141 592 653 40… < テンプレート:Π < 3.141 592 653 86… となり、テンプレート:Π の近似値を小数点以下第9位まで得る。

次に n = 4 の場合は、

$\frac{1}{64}\underset{0}{\overset{1}{\int }}\frac{x^{16}(1-x{)}^{16}}{1+x^2}dx=\pi -\frac{741269838109}{235953517800}$

および

$\frac{1}{128}\underset{0}{\overset{1}{\int }}{x}^{16}(1-x{)}^{16}dx=\frac{1}{2538963567360}$

より

$\frac{741269838109}{235953517800}+\frac{1}{2538963567360}<\pi <\frac{741269838109}{235953517800}+\frac{2}{2538963567360}$

すなわち 3.141 592 653 589 55… < テンプレート:Π < 3.141 592 653 589 95… となり、テンプレート:Π の近似値を小数点以下第12位まで得る。

関連文献

• テンプレート:Cite book
  • テンプレート:Cite book
• Lucas, Stephen (2007), "Integral approximations to テンプレート:Π with nonnegative integrands"
• テンプレート:Citation

関連項目

• 円周率の無理性の証明
• 円周率の歴史
• 世界の名著
• リンデマンの定理 -- 円周率が超越数であることを含意する定理

外部リンク

• テンプレート:MathWorld
• Jonathan Borwein, The Life of Pi