円周率を評価する

同義語:円周率の評価estimating pi

概要

円周率を評価するとは、有理数や根号で書ける数 $a<b$ について $a<\pi<b$ を証明することである。近似値 $\pi\approx\frac{22}7$ はずれの向きが分からないので評価ではなく、有限個の評価から $\pi$ が無理数だとは分からない。高校の道具では、単位円に内接・外接する正多角形の周で挟む方法($\sin x<x<\tan x$ による)、$\int_0^1\frac{x^4(1-x)^4}{1+x^2}dx=\frac{22}7-\pi$ のように「有理数 $-\pi$」を小さな正の積分として設計する方法、$\arctan$ の級数を小さな $x$ で使う Machin の公式 $\frac\pi4=4\arctan\frac15-\arctan\frac1{239}$ の 3 つがあり、この順に精度が上がる。本記事はこの 3 本の記事への案内である。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 円周率, 三角関数, 定積分

この記事について

円周率 $\pi=3.14159\ldots$ の値を知っていることと、「$3.14<\pi<3.15$ である」ことを証明できることは別である。この記事は、高校で習う道具だけで $\pi$ を有理数や根号で挟む方法を学ぶ 3 本の記事への入口である。ここでは「評価する」とは何かを例で確かめ、各記事で何が分かるかと読む順番を案内する。証明はそれぞれの記事で行う。

円周率を評価するとは

円周率を評価する

有理数や根号で書ける数 $a< b$ について、$a<\pi< b$ が成り立つことを証明することを、円周率を評価する(挟む)という。$b-a$ が小さいほど、精密な評価である。

正方形と正六角形で挟む

単位円(半径 $1$)の周の長さは $2\pi$ である。円に内接する正六角形は 1 辺 $1$ の正三角形 6 つでできているので、周の長さは $6$ である。外接する正方形の周の長さは $8$ である。「内接する多角形の周 $<$ 円周 $<$ 外接する多角形の周」から、$6<2\pi<8$、つまり $3<\pi<4$ である。外接する正方形の代わりに外接する正六角形(周 $4\sqrt3$)を使えば、$3<\pi<2\sqrt3=3.464\ldots$ になる。
この不等式がなぜ正しいかは、弦・弧・接線の長さの比較 $\sin x< x<\tan x$($0< x<\frac\pi2$)で証明できる(正多角形で円周率を挟む)。

反例:近似値は評価ではない

「$\pi\approx3.14$」や「$\pi\approx\frac{22}7$」は近似値であって、評価ではない。どちら側にずれているかが分からないからである。実際、$3.14<\pi$ だが $\frac{22}7=3.1428\ldots>\pi$ であり、ずれの向きは近似値ごとに違う。$\pi<\frac{22}7$ を示すには証明が要る(積分で円周率を評価する)。

反例:評価をいくら重ねても無理数とは分からない

$3.1415926<\pi<3.1415927$ のような評価は、$\pi$ が分母の大きな分数である可能性を否定しない。たとえば $\frac{31415926}{10000000}$ と $\frac{31415927}{10000000}$ の間には、無数の分数がある。$\pi$ が分数でない(無理数である)ことは、有限個の評価からは決して分からず、別の証明が要る(円周率の無理性)。さらに $\pi$ は、$0$ でない有理数係数の多項式の根にもならない超越数であることが知られている(証明しない。超越数)。

方法の地図

方法何を比べるかなぜ効くか誤差の目安記事
正多角形で挟む弧と弦、弧と接線$\sin x< x<\tan x$正 $n$ 角形で $\frac1{n^2}$ 程度正多角形で円周率を挟む
積分で挟む小さな正の関数の積分正の関数の積分は正設計しだいで非常に小さい積分で円周率を評価する
速く収束する級数$\arctan$ の級数の打ち切り誤差が $x^{2n+1}$ で減るMachin の公式なら 1 項で約 1.4 桁Machinの公式
無理数性の証明「整数」と「小さい正の数」$0$ と $1$ の間に整数はない(評価ではなく不可能性)円周率の無理性

図1:方法ごとに得られる、π を挟む区間の幅(対数目盛)。正六角形の 0.46 から、Machin の公式(各 8 項)の 2.5×10^(-12) まで大きく違う 図1:方法ごとに得られる、π を挟む区間の幅(対数目盛)。正六角形の 0.46 から、Machin の公式(各 8 項)の 2.5×10^(-12) まで大きく違う
図 1 の棒の長さは、各方法で得られる $a<\pi< b$ の幅 $b-a$ である。上から順に、内接・外接正六角形による $3<\pi<2\sqrt3$(ex-epi-guide-hexagon)、四分円を幅 $\frac1{10}$ の長方形で挟んだ場合(ex-epi-guide-naive)、内接・外接正九十六角形の周による場合(ex-epi-guide-polygon)、$\frac{22}7$ の積分表示による $\frac{22}7-\frac1{630}<\pi<\frac{22}7-\frac1{1260}$(ex-epi-guide-integral)である。最後の 2 本は、Machin の公式で各 3 項・各 8 項まで使った近似値から、誤差の上界を引いた数と足した数で $\pi$ を挟んだ場合で、幅は誤差の上界の 2 倍になる(ex-epi-guide-machin)。

高校の計算大学での見方
微分して $\sin x< x<\tan x$ を示す弧の長さは内接する折れ線の長さの上限
多項式を $x^2+1$ で割る複素数 $i$ を代入すること
等比数列の和の公式を積分する冪級数と誤差の評価
正接の加法定理複素数の積と偏角

各記事の代表的な結果

正十二角形による評価

単位円に内接・外接する正十二角形の周の半分は $3(\sqrt6-\sqrt2)=3.1058\ldots$ と $12(2-\sqrt3)=3.2153\ldots$ であり、$3(\sqrt6-\sqrt2)<\pi<12(2-\sqrt3)$ が成り立つ。辺の数を 2 倍にしていくと、正九十六角形で $3.1410<\pi<3.1428$ となる。Archimedes はこの方法で $\frac{223}{71}<\pi<\frac{22}7$ を示した(正多角形で円周率を挟む)。

22/7 と π の差を表す積分

$$ \int_0^1\frac{x^4(1-x)^4}{1+x^2}\,dx=\frac{22}7-\pi $$
であり、被積分関数は正なので $\pi<\frac{22}7$ である。分母を $1$ と $2$ で挟むと $\frac{22}7-\frac1{630}<\pi<\frac{22}7-\frac1{1260}$、つまり $3.14126\ldots<\pi<3.14206\ldots$ となる(積分で円周率を評価する)。

反例:素朴に挟むと粗い

積分で $\pi$ を挟むとき、被積分関数をそのまま別の関数で挟むだけでは粗い。たとえば $0< x<1$ で $\frac12<\frac1{1+x^2}<1$ なので、$\frac\pi4=\int_0^1\frac{dx}{1+x^2}$ から $2<\pi<4$ しか出ない。四分円 $\int_0^1\sqrt{1-x^2}\,dx=\frac\pi4$ を幅 $\frac1{10}$ の長方形で挟んでも $2.904<\pi<3.305$ である。精度を上げるには、$\frac{22}7-\pi$ のような「小さな量」そのものを積分で表す設計が要る(積分で円周率を評価する)。

反例:Leibniz の級数は遅い

$\frac\pi4=1-\frac13+\frac15-\frac17+\cdots$(Leibniz の級数)は正しいが、1000 項まで足して 4 倍しても、$\pi$ との差は約 $0.001$ 残る。項を 1 つ増やしても誤差がほとんど減らないからである。同じ $\arctan$ の級数でも、$x=\frac15$ や $x=\frac1{239}$ のような小さな値で使うと、項を 1 つ増やすごとに誤差の上界は $\frac1{25}$ 倍以下になる(Machinの公式)。

Machin の公式で 2 項ずつ

$\frac\pi4=4\arctan\frac15-\arctan\frac1{239}$ の $\arctan$ を、それぞれ $x-\frac{x^3}3$ の 2 項で近似すると
$$ \pi\approx16\Bigl(\frac15-\frac1{375}\Bigr)-4\Bigl(\frac1{239}-\frac1{3\cdot239^3}\Bigr)=3.1405970\ldots $$
となり、実際の誤差は $0.000995\ldots$ である。項を増やしたときの誤差は、各 3 項で $2.93\times10^{-5}$ 未満、各 8 項で $1.24\times10^{-12}$ 未満であることが証明できる(Machinの公式。各 8 項での実際の誤差は $1.19\times10^{-12}$ ほど)。

読む順番

  1. 正多角形で円周率を挟む:$\sin x< x<\tan x$ を微分で証明し、内接・外接正多角形の周で $\pi$ を挟む。辺の数を 2 倍にする漸化式(調和平均と相乗平均)と、誤差が $\frac1{n^2}$ 程度でしか減らないことが分かる。
  2. 積分で円周率を評価する:正の関数の積分が正であることを使い、$\frac{22}7-\pi$ を小さな正の積分として表す。割り算の余りを $x=i$ の値で求めて、どんな被積分関数なら $\pi$ が現れるかが分かる。
  3. Machinの公式:$\arctan$ の級数の誤差を積分で評価し、$x=1$ では遅いこと、小さな $x$ を組み合わせる Machin の公式で速くなることを証明する。
    1 と 2 はどちらから読んでもよい。3 は、1 で証明する $\theta<\tan\theta$ と、2 で証明する「正の関数の積分は正」を 1 つずつ使う。1 → 2 → 3 の順に、評価の精度が上がっていく。$\pi$ が無理数であることの証明は、2 と同じ「小さな正の積分」の考え方を使う 円周率の無理性 にある。

関連項目

参考文献

[4]
Stephen K. Lucas, Integral proofs that 355/113 > π, Gazette of the Australian Mathematical Society, 2005, pp. 263–266
[7]
Johann Heinrich Lambert, Mémoire sur quelques propriétés remarquables des quantités transcendantes circulaires et logarithmiques, Histoire de l'Académie Royale des Sciences et des Belles-Lettres de Berlin, 1768, pp. 265–322

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する