円周率

テンプレート:Redirect テンプレート:円周率 円周率（えんしゅうりつ、テンプレート:Lang-en-short、テンプレート:Lang-de-short、テンプレート:Lang-zh-short）とは、円の直径に対する円周の長さの比率のことをいいテンプレート:Sfn、数学定数の一つである。通常、円周率はギリシア文字であるテンプレート:Math テンプレート:Refnestで表される。円の直径から円周の長さや円の面積を求めるときに用いるテンプレート:Sfn。また、数学をはじめ、物理学、工学といった科学の様々な理論の計算式にも出現し、最も重要な数学定数とも言われる^[1]。

円周率は無理数であり、超越数でもある。

円周率の計算において功績のあったルドルフ・ファン・クーレンに因み、ルドルフ数とも呼ばれる。ルドルフは小数点以下35桁まで計算した^[2]。小数点以下35桁までの値は次の通りである。

π = 3.14159 26535 89793 23846 26433 83279 50288 \dots

基礎

表記と呼び方

円周率を表すギリシア文字テンプレート:Mvar は、ギリシア語でいずれも周辺・円周・周を意味するテンプレート:Lang テンプレート:Sfn^[3]（ペリメトロス）あるいはテンプレート:Lang^[4]（ペリペレイア）の頭文字から取られたテンプレート:Refnest。文字テンプレート:Π をウィリアム・オートレッドは1631年に著した著書において半円周の長さを表す文字として用い、アイザック・バローは論文において半径テンプレート:Mvar の円周の長さとして用いた^[5]。ウィリアム・ジョーンズ (1706) やレオンハルト・オイラーらにより（現代と同じく）円周の直径に対する比率を表す記号として用いられ、それが広まったテンプレート:Sfn^[5]。日本では「パイ」と発音する^[3]。

数テンプレート:Π を指す言葉には、日本・中国・韓国における「円周率（圓周率）」、ドイツの「Kreiszahl」（Kreis は円（周）、Zahl は数の意）の他、それを計算した人物の名前を取った「アルキメデス数」（テンプレート:Lang-en-short）、「ルドルフ数」（テンプレート:Lang-en-short、テンプレート:Lang-de-short）などがある。一般にドイツ語を除くヨーロッパの諸言語には「円周率」に対応する単語はない^[3]^[6]。

なお、「テンプレート:Π」の字体は、表示環境によってはキリル文字の [[П|テンプレート:Lang]] に近いテンプレート:Lang などと表示されることがある。また、ギリシャ文字「テンプレート:Π」は、円周率とは無関係に、素数計数関数や、基本群・ホモトピー群、ある種の写像（射影など）を表すのに用いられることもある。

定義

平面幾何学において、円周率テンプレート:Π は、円の周長の直径に対する比率として定義される。すなわち、円の周長をテンプレート:Mvar, 直径をテンプレート:Mvar としたとき、

π = \frac{C}{d}

である。全ての円は互いに相似なので、この比率は円の大きさに依らず一定である。

ところが、この定義は円の周長を用いているため、曲線の長さを最初に定義していない解析学などの分野では、テンプレート:Π が現れる際に問題となることがある。この場合、円の周長に言及せず、解析学などにおける性質の一つをテンプレート:Π の定義とすることが多いテンプレート:Sfn。この際のテンプレート:Π の定義の一般なものとして、三角関数テンプレート:Math がテンプレート:Math を取るようなテンプレート:Math2 の最小値の2倍とするもの、級数による定義、定積分による定義などがある。後述の#円周率に関する式も参照。

歴史

テンプレート:See also テンプレート:Multiple image

古代

円周の直径に対する比率は円の大きさに依らず一定であり、それは 3 より少し大きい^{[注 1]}ことは古代エジプトやバビロニア、インド、ギリシアの幾何学者たちにはすでに知られていた。また、古代インドやギリシアの数学者たちの間では半径テンプレート:Mvar の円板の面積がテンプレート:Math であることも知られていた。さらに、アルキメデスは正テンプレート:Math角形を用いて半径テンプレート:Mvar の球の体積がテンプレート:Math であることや、この球の表面積がテンプレート:Math（その球の大円による切り口の面積の4倍）であることを導き出した。

円周率を小数で最初に記述したのは小数を発明した中国である。263年に魏の劉徽が3072角形を使用し3.14159と計算し、5世紀に祖沖之が十尺もの直径の円を使用して3.1415926<テンプレート:Π<3.1415927 と求め、以後1000年間、全世界でこれ以上正確な計算はなされなかった。祖の計算が正確であったことは、1300年頃に趙友欽が16384辺の内接多角形により確かめた^[7]。

近代まで

テンプレート:出典の明記 14世紀インドの数学者・天文学者であるサンガマグラーマのマーダヴァは次のテンプレート:Π の級数表示を見いだしている（ライプニッツの公式）：

\frac{π}{4} = 1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7} + \dots = \sum_{n = 0}^{\infty} \frac{(- 1)^{n}}{2 n + 1}

これは逆正接関数テンプレート:Math のテイラー展開のテンプレート:Math2 での表式になっている。マーダヴァはまた、

π = \sqrt{12} (1 - \frac{1}{3 \cdot 3} + \frac{1}{5 \cdot 3^{2}} - \frac{1}{7 \cdot 3^{3}} + \dots)

を用いてテンプレート:Π の値を小数点以下11桁まで求めている。

17世紀、ドイツのルドルフ・ファン・コーレンが正テンプレート:Math億角形を使い、小数点以下第35位まで計算。1699年（または1706年）にエイブラハム・シャープが小数点以下第72～127位まで求めた。

18世紀フランスの数学者アブラーム・ド・モアブルは、ある定数テンプレート:Mvar を取ると、コインをテンプレート:Math回投げて表がテンプレート:Mvar回だけ出る確率は、テンプレート:Mvar が十分大きいとき

\frac{C}{\sqrt{n}} \exp {- \frac{(x - n)^{2}}{n}}

で近似できることを、テンプレート:Math2 における数値計算により見いだした。この正規分布の概念は1738年に出版されたド・モアブルの『巡り合わせの理論』に現れている。ド・モアブルの友人のジェイムズ・スターリングは後に、 $C = \frac{1}{\sqrt{2 π}}$ であることを示した。

1751年にヨハン・ハインリヒ・ランベルトは、テンプレート:Mvar がテンプレート:Math でない有理数ならば正接関数テンプレート:Math の値は無理数であることを示し、その系（の対偶）としてテンプレート:Π は無理数であることを導いた。さらに1882年にフェルディナント・フォン・リンデマンはテンプレート:Π が超越数であることを示し、円積問題（与えられた長さを半径とする円と等積の正方形を作図する問題）は解くことができないことを導いた。

1873年、ウィリアム・シャンクスは彼自身の手で小数点以下第707位までを計算した（ただしその結果は途中で生じた誤りにより小数点以下第527位までしか正しくなかった）。

和算における円周率の取り扱い

江戸時代の初期の和算家の3.16

江戸時代初期の数学書である毛利重能の『割算書』では円周率を3.16としている。その弟子の吉田光由の『塵劫記』でも3.16となっているテンプレート:Sfn。しかし、当時の先進国中国の文献にはこの3.16という数値は見られず、中国の文献の数値を引き写したとは考えにくいというテンプレート:Sfn。そのため、なぜ初期の和算家が円周率を3.16としたかの理由はよく分かっていないテンプレート:Sfn。おそらく、毛利重能とその弟子の吉田光由などの先駆者らは、円周率を実際に測定して3.14ないし3.16ほどの値を得たが、最後の桁の数字に確信が持てなかったため、「円のような美しい形を求める数値は、もっと美しい数値になっていいはずだ」と考え、「美しい理論」を求めた。その結果テンプレート:Math2 が美しい数値として採用されたと推測されているテンプレート:Sfn。その考えは日本で2番目に3.14の値を計算で求めた野沢定長の『算九回』（延宝五年:1677年）の中にも見られ、その著書の中で「忽然として円算の妙を悟った」として「円周率の値は形＝経験によって求めれば3.14であるが、理＝思弁によって求めれば3.16である」として「両方とも捨てるべきでない」としたテンプレート:Sfn。

和算家が計算した3.14

江戸初期、1600年代前半頃から、円を対象とした和算的研究である「円理」が始まる。その最初のテーマの一つが円周率を数学的に計算する努力でありテンプレート:Sfn、1663年に日本で初めて村松茂清が『算爼（さんそ）』において「円の内接多角形の周の長さを計算する方法」で3.14…という値を算出したテンプレート:Sfn。『算爼』では円に内接する正8角形から角数を順次2倍していき、内接テンプレート:Math角形の周の長さで、

3.1415 9264 8777 6988 6924 8

と小数点以下21桁まで算出している。これは実際の値と小数第7位まで一致しているテンプレート:Sfn。その後1680年代に入ると、円周率の値を3.16とする数学書はなくなり、3.14に統一されたテンプレート:Sfn。1681年頃には関孝和が内接テンプレート:Math角形（13万1072角形）の計算を工夫し、小数第16位まで現代の値と同じ数値を算出した。この計算値は関の死後1712年に刊行された『括要算法』に記されているテンプレート:Sfn。

日本の和算家に特徴的なのは、1663年に3.14が初めて導き出されても、その後1673年までの10年間に円周率の値を3.14とした算数書のいずれもが、先行者の円周率をそのまま引き継ぐことをせず、それぞれ独自の値を提出していたことであるテンプレート:Sfn。この背景には当時の遺題継承テンプレート:Refnest運動に「他人の算法をうけつぐ」と共に「自己の算法を誇る」という性格があったためだというテンプレート:Sfn。そのため古い3.16の値が疑われてから、遺題継承の際に必ずといってよいほど円周率の値が変えられているテンプレート:Sfn。しかしながら江戸時代の3大和算書『塵劫記』『改算記』『算法闕疑抄』の増補改訂版では1680年代には3.14に統一されたテンプレート:Sfn。

3.14から3.16への逆行

しかし、遺題継承運動は1641年に始まって1699年頃には終わってしまいテンプレート:Sfn、いったん3.14に統一された円周率の値は江戸時代後半になると揺らぎ始め、古い3.16に逆行するという現象が生じたテンプレート:Sfn。文政年間（1818～30年）に出版された算数書とソロバン書を悉皆調査した結果では、円周率の値を3.14とするものと、3.16とするものの2系統があることが明らかにされたテンプレート:Sfn。いくらか専門的な数学書では3.14とされているのに、大衆向けの小冊子の中では3.16の方が普通に用いられていたテンプレート:Sfn。

当時の識者である橘南谿（1754-1806年）は「いまに至り3.16あるいは3.14色々に論ずれども、なおきわめがたきところあり」と述べ、3.14はまだ確定していないとしているテンプレート:Sfn。儒学者の荻生徂徠も和算家の算出した3.14の根拠に納得しなかったテンプレート:Sfn。当時の和算家のほとんどは、円に内接する多角形の周を計算することで円周率を計算した。内接多角形の角数を増やすほど求まる円周率の桁は増えていくので、素人目にはその値が増大する一方に見える。「それがいくら増えても3.1416を超えない」ということを和算家たちはついに納得させることができなかったのであるテンプレート:Sfn。

そのような和算家以外の素人たちを納得させるには、どうしても万人に納得させる「理」に基づいて計算してみせる他はないテンプレート:Sfn。それを行うには西洋で行われたように、「円を内接多角形と外接多角形ではさんで、円周率の上限と下限を示すこと」が必要であったが、（次の鎌田による成果を例外として）和算家はついにその方法を取ることがなかったテンプレート:Sfn。

宅間流和算の円周率

日本で唯一「円周を内接・外接多角形で挟み込んで円周率の上限と下限を示す」ことに成功したのは鎌田俊清（1678-1747年）が享保七年（1722年）に著した『宅間流テンプレート:Refnest円理』である。その値は以下の通りであるテンプレート:Sfn。

内周：3.1415 9265 3589 7932 3846 2643 3665 8

外周：3.1415 9265 3589 7932 3846 2643 4166 7

鎌田は円周率の小数点以下24桁まで正しいと確信しうる円周率の値を算出することに成功していたテンプレート:Sfn。しかし、鎌田の方法は後継者を持たず、当時の識者に知られることがなかったテンプレート:Sfn。

級数展開による算出

日本の数学史では級数による値の算出は広く一般的であった。円周率の級数による公式は多くの学者に研究されており、蜂谷定章、松永良弼、坂部広畔、川井久徳、長谷川寛らによるものがある^[8]。また、建部賢弘は円周率の二乗を求める日本初の公式を考案した^[9]。

和算の限界

日本の和算の弱点は単に理論面の弱さにとどまらず、万人が納得できる正しい円周率の教育・啓蒙への関心も失ったことであったテンプレート:Sfn。そのため和算家たちがいくら円周率は3.14…と書いたところで、『塵劫記』の古い円周率3.16の値がそのまま残存する結果となったテンプレート:Sfn。『塵劫記』の重版（1694年）などは古い円周率3.16のまま出版され続け、18世紀に大衆的な通俗算数書が大量に出版される際に、必ずというほど3.16という値を引き継ぐようになってしまったテンプレート:Sfn。

18世紀半ば以降の和算は数学的証明の概念の追求は無視され、せっかく宅間流の鎌田俊清がその独創的方法で正しい円周率を算出しても、全く継承されなかったテンプレート:Sfn。江戸時代後半の和算家は家元制度的な秘密主義と保守主義と、権威主義が在野の独創性を無視し、結果として学問の進歩を妨げることとなったテンプレート:Sfn。

コンピュータの利用

20世紀以降、計算機の発達により、計算された円周率の桁数は飛躍的に増大した。1949年に、電子計算機ENIACを使い72時間かけて、円周率は2037桁まで計算された^[10]。その後の数十年間、様々な計算機科学者や計算科学者など、あるいはコンピュータのアマチュアによって計算は進められ、1973年には100万桁を超えた。この進歩は、スーパーコンピュータの開発だけによるものではなく、効率の良いアルゴリズムが考案されたためである。そのうちの最も重要な発見の一つとして、1960年代の高速フーリエ変換がある。これにより、多倍長の演算が高速に実行できるようになった。

2022年6月9日に、Googleの技術者、岩尾エマはるかがGoogle Cloudで、チュドノフスキー級数を使い、157日23時間かけて100兆桁を計算したと発表^[11]。

性質

無理性

テンプレート:Π は無理数であるため、循環しない無限小数である。

テンプレート:For2 テンプレート:Π は無理数である。つまり、2つの整数の商で表すことはできず、小数展開は循環しない。このことは1761年にヨハン・ハインリヒ・ランベルトが証明したが、厳密性に欠けた部分があった。その部分は1806年にルジャンドルによって補われた。

したがって、円周率のコンピュータによる計算や暗唱、十進法表示での小数部分の各数字 (0, 1, …, 9) の出現頻度は、人々の興味の対象となる。

超越性

テンプレート:For2 さらに、テンプレート:Π は超越数である。つまり、有理数係数の代数方程式の解にはならない。これは1882年にフェルディナント・フォン・リンデマンによって証明された（リンデマンの定理）。これより、整数から四則演算と冪根をとる操作だけを有限回組み合わせてもけっしてテンプレート:Π の値は得られないことが分かる。

テンプレート:Π が超越数であることから直ちに、古代ギリシアの三大作図問題の内の一つである「円積問題」（与えられた長さを半径とする円と等積の正方形を目盛の無い定規とコンパスを「有限回」用いて作図すること）は不可能なことが結論される。

ランダム性

2022年10月の時点で、テンプレート:Π は小数点以下100兆桁まで計算されている^[11]。そして、分かっている限りでは 0 から 9 までの数字がランダムに現れているようには見えるが、それが乱数列といえるかどうかははっきりとは分かっていない。たとえばテンプレート:Π が正規数であるかどうかも分かっていない。正規数であればテンプレート:Π の10進表示において、各桁を順に取り出して得られる数列^[12]：

3, 1, 4, 1, 5, 9, 2, 6, 5, 3, 5, …

には、0 から 9 が均等に現れるはずだが分かっておらず、それどころか、0 から 9 がそれぞれ無数に現れるのかどうかすら分かっていない。もし仮に正規数でないとすれば、乱数列でもないということになる。

5兆桁までの数字の出現回数は以下の通りである。全てほぼ等しく（約0.0005%の違いに収まる）、最も多いのは 8 で、最も少ないのは 6 である。

0：4999億9897万6328回

1：4999億9996万6055回

2：5000億0070万5108回

3：5000億0015万1332回

4：5000億0026万8680回

5：4999億9949万4448回

6：4999億9893万6471回

7：5000億0000万4756回

8：5000億0121万8003回

9：5000億0027万8819回

連分数

テンプレート:See also 分母を整数と分数の和で表すことを続けていった表示を連分数という。「整数」を最大にしていくと、分子を全て 1 にできる：

a_{0} + \frac{1}{a_{1} + \frac{1}{a_{2} + \frac{1}{⋱ + \frac{1}{a_{n}}}}}

テンプレート:Pi は無理数であるから、円周率テンプレート:Pi の連分数展開は有限項では終わらず無限項の連分数となる：

π = 3 + \frac{1}{7 + \frac{1}{15 + \frac{1}{1 + \frac{1}{292 + \frac{1}{1 + \frac{1}{1 + \frac{1}{1 + ⋱}}}}}}}

上記の正則連分数展開（すべての分子が 1 である連分数）を途中で打ち切ると、テンプレート:Pi の良い有理数近似が得られる。その最初の4つはテンプレート:Math2 である。これらは古くからよく知られ使用されてきた近似値である。これらはそれぞれ分母が大きくないどの分数よりもテンプレート:Pi に近く^[13]、テンプレート:Pi の最良有理数近似である。

さらに、テンプレート:Pi は超越数である（つまり代数的数でない）ことが知られている。一般に、正則連分数の分母に現れる整数部が循環するのはテンプレート:Ill2（有理係数の二次方程式の解となるような代数的数）に限られ、テンプレート:Pi は二次無理数でないためテンプレート:Ill2として表せない。加えてテンプレート:Pi の正則連分数は規則性を示さないが^[14]、テンプレート:Pi のテンプレート:Ill2では以下の規則をもつものが知られている^[15]：

\begin{matrix} π & = \frac{4}{1 + \frac{1^{2}}{2 + \frac{3^{2}}{2 + \frac{5^{2}}{2 + \frac{7^{2}}{2 + \frac{9^{2}}{2 + ⋱}}}}}} = 3 + \frac{1^{2}}{6 + \frac{3^{2}}{6 + \frac{5^{2}}{6 + \frac{7^{2}}{6 + \frac{9^{2}}{6 + ⋱}}}}} \\ = \frac{4}{1 + \frac{1^{2}}{3 + \frac{2^{2}}{5 + \frac{3^{2}}{7 + \frac{4^{2}}{9 + ⋱}}}}} \end{matrix}

未解決問題

テンプレート:Main

次の数は超越数か？

π \pm e, π e, \frac{π}{e}, π^{π}, π^{e}, π^{\sqrt{2}}, e^{π^{2}}

ただし、 $π \pm e$ と $π e, \frac{π}{e}$ は両者少なくとも一方は超越数であることは分かっているテンプレート:要出典。

テンプレート:Π は正規数か？

円周率に関する式

テンプレート:Π についての式は非常に多い。ここではその一部を紹介する。数式によってはそれ自体がテンプレート:Π の定義になり得るし、テンプレート:Π の近似値の計算などにも使われてきた。

幾何

半径テンプレート:Mvar の円の周長： $2 π r$
半径テンプレート:Mvar の円の面積： $π r^{2}$
半径テンプレート:Mvar の球の体積： $\frac{4}{3} π r^{3}$
半径テンプレート:Mvar の球の表面積： $4 π r^{2}$
長半径テンプレート:Mvar, 短半径テンプレート:Mvar の楕円の面積： $π a b$
$18 0^{\circ} = π$ rad

解析（特殊関数と虚数を除く）

ライプニッツの公式、#2千年紀、逆正接関数も参照

\sum_{n = 0}^{\infty} \frac{(- 1)^{n}}{2 n + 1} = 1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7} + \dots = \frac{π}{4}

マーダヴァの公式、逆正接関数も参照

\sqrt{12} \sum_{n = 0}^{\infty} \frac{(- 1)^{n}}{3^{n} (2 n + 1)} = \sqrt{12} (1 - \frac{1}{3 \cdot 3} + \frac{1}{5 \cdot 3^{2}} - \frac{1}{7 \cdot 3^{3}} + \dots) = π

ウォリスの公式

\prod_{n = 1}^{\infty} \frac{(2 n)^{2}}{(2 n - 1) (2 n + 1)} = \frac{2^{2}}{1 \cdot 3} \cdot \frac{4^{2}}{3 \cdot 5} \cdot \frac{6^{2}}{5 \cdot 7} \cdot \frac{8^{2}}{7 \cdot 9} \dots = \frac{π}{2}

$\prod_{n = 1}^{\infty} \frac{n^{2} + n}{n^{2} + n + \frac{1}{4}} = \frac{8}{9} \cdot \frac{24}{25} \cdot \frac{48}{49} \cdot \frac{80}{81} \cdot \frac{120}{121} \cdot \frac{168}{169} \cdot \frac{224}{225} \cdot \frac{288}{289} \dots = \frac{π}{4}$
ビエトの公式

\prod_{n = 1}^{\infty} \cos \frac{9 0^{\circ}}{2^{n}} = \sqrt{\frac{1}{2}} \sqrt{\frac{1}{2} + \frac{1}{2} \sqrt{\frac{1}{2}}} \sqrt{\frac{1}{2} + \frac{1}{2} \sqrt{\frac{1}{2} + \frac{1}{2} \sqrt{\frac{1}{2}}}} \dots = \frac{2}{π}

$\sum_{n = 1}^{\infty} \frac{1}{n^{2} 2^{n - 1}} + (\log 2)^{2} = \frac{π^{2}}{6}$ （オイラー）

$\int_{- \infty}^{\infty} e^{- x^{2}} d x = \sqrt{π}$ （ガウス積分）
$π = 2 \int_{0}^{1} \frac{d t}{\sqrt{1 - t^{2}}}$ テンプレート:Sfn
$π = \int_{- 1}^{1} \frac{d t}{\sqrt{1 - t^{2}}}$
$π = 2 \int_{- 1}^{1} \sqrt{1 - t^{2}} d t$
$π = 4 \int_{0}^{1} \frac{d t}{1 + t^{2}}$
$π = \int_{- \infty}^{\infty} \frac{d t}{1 + t^{2}}$
逆三角関数は主値を取るものとすると

\begin{matrix} π & = 2 \arccos 0 \\ = 2 \arcsin 1 \\ = 4 \arctan 1 \end{matrix}

逆三角関数（逆正弦関数）の公式より

π = 2 \sum_{n = 0}^{\infty} \frac{(2 n - 1)!!}{(2 n + 1) (2 n)!!}

逆三角関数（逆正接関数）の公式より
- 逆正接関数のテイラー展開による： $\begin{matrix} π & = 4 \sum_{n = 0}^{\infty} \frac{(- 1)^{n}}{2 n + 1} \end{matrix}$
- オイラーによる^[16]： $\begin{matrix} π & = 2 \sum_{n = 0}^{\infty} \frac{n!}{(2 n + 1)!!} \\ = \sum_{n = 0}^{\infty} \frac{2^{n + 1} (n!)^{2}}{(2 n + 1)!} \end{matrix}$
双曲線関数（双曲線余接関数）の公式より

\frac{1}{e^{2} - 1} = \sum_{n = 1}^{\infty} \frac{1}{(n π)^{2} + 1}

ニュートンの無平方根公式^[17]

\begin{matrix} π & = 3 \sum_{n = 0}^{\infty} \frac{(2 n)!}{(2 n + 1) 1 6^{n} (n!)^{2}} \\ = 3 \sum_{n = 0}^{\infty} \frac{(n + 1) C_{n}}{(2 n + 1) 1 6^{n}} \end{matrix}

（テンプレート:Mvar はカタラン数）この式は、

π = 6 \arcsin \frac{1}{2}

のマクローリン級数となっている^[18]。

マチンの公式（マチン、1706年^[18]^[19]（1709年とも））

4 \arctan \frac{1}{5} - \arctan \frac{1}{239} = \frac{π}{4}

4 arccot 5 - arccot 239 = \frac{π}{4}

と書かれることもある。

テンプレート:Math とテンプレート:Math が二進法と相性が良く、収束も早いため、コンピュータでの円周率計算によく使われる公式の一つである。

テンプレート:Math の連分数表示

\frac{4}{π} = 1 + \frac{1}{3 + \frac{4}{5 + \frac{9}{7 + \frac{16}{9 + \frac{25}{⋱}}}}}

ガウス＝ルジャンドルのアルゴリズム

初期値の設定：

a_{0} = 1, b_{0} = \frac{1}{\sqrt{2}}, t_{0} = \frac{1}{4}, p_{0} = 1 .

反復式：テンプレート:Math2 が希望する桁数になるまで以下の計算を繰り返す。小数第テンプレート:Mvar位まで求めるときテンプレート:Math回程度の反復でよい。

\begin{matrix} a_{n + 1} & = \frac{a_{n} + b_{n}}{2}, \\ b_{n + 1} & = \sqrt{a_{n} b_{n}}, \\ t_{n + 1} & = t_{n} - p_{n} (a_{n} - a_{n + 1})^{2}, \\ p_{n + 1} & = 2 p_{n} . \end{matrix}

テンプレート:Π の算出：円周率テンプレート:Π は、テンプレート:Math2 を用いて以下のように近似される。

π \approx \frac{1}{4 t_{n}} (a_{n} + b_{n})^{2}

非常に収束が早く^{[注 2]}、金田康正が1995年に42億桁、2002年に1.24兆桁を計算した[[スーパーπ|スーパーテンプレート:Π]] に使われていた。

$n! \sim \sqrt{2 π n} {(\frac{n}{e})}^{n}$ （スターリングの近似。テンプレート:Math は $\lim_{n \to \infty} \frac{f (n)}{g (n)} = 1$ を表す）
$\frac{1}{π} = \frac{2 \sqrt{2}}{9 9^{2}} \sum_{n = 0}^{\infty} \frac{(26390 n + 1103) \cdot (4 n)!}{(4^{n} 9 9^{n} \cdot n!)^{4}}$ （ラマヌジャン）
$\frac{4}{π} = \sum_{n = 0}^{\infty} \frac{(- 1)^{n} (21460 n + 1123) \cdot (4 n)!}{88 2^{2 n + 1} (4^{n} n!)^{4}}$ （ラマヌジャン）
$\sum_{n = 1}^{\infty} \frac{n}{e^{2 π n} - 1} = \frac{1}{24} - \frac{1}{8 π}$ （ラマヌジャン）
$\frac{1}{π} = \frac{12}{C_{0} \sqrt{C_{0}}} \sum_{n = 0}^{\infty} \frac{(C_{2} n + C_{1}) \cdot (6 n)!}{(- C_{0})^{3 n} \cdot (3 n)! \cdot (n!)^{3}}$ （テンプレート:仮リンク）^[20]

（各定数と、その素因数分解：

テンプレート:Math

テンプレート:Math）

$\frac{1}{π} = \frac{12}{C_{2} \sqrt{C_{2}}} \sum_{n = 0}^{\infty} \frac{(- 1)^{n} (6 n)! (C_{0} + C_{1} n)}{(3 n)! (n!)^{3} C_{2}^{n}}$ (テンプレート:仮リンク, テンプレート:仮リンク)^[20]

（各定数の値：

テンプレート:Math

テンプレート:Math）

David Bailey, Peter Borwein, およびサイモン・プラウフによるもの（ベイリー＝ボールウェイン＝プラウフの公式・俗称 "BBP"）、Adamchik と Wagon によるもの、Fabrice Bellard によるもの^[21]^[18]等については、あまりに高度なため割愛する。

複素解析

$e^{i π} + 1 = 0$ （オイラーの等式）
$\sum_{k = 0}^{n - 1} e^{i \cdot \frac{2 π k}{n}} = 0$ （テンプレート:Mvar はテンプレート:Math 以上の整数）

後者はオイラーの等式の一般化であり、[[1の冪根|テンプレート:Math のテンプレート:Mvar乗根]]の総和はテンプレート:Math になることを示している。テンプレート:Math2 とするとオイラーの等式になる。

$π = 4 + \sum_{n = 1}^{\infty} \frac{1}{n!} (\frac{1}{n + λ} - \frac{4}{2 n + 1}) {(\frac{(2 n + 1)^{2}}{4 (n + λ)} - n)}_{n - 1}$ （テンプレート:仮リンク、テンプレート:仮リンク）^[22]^[23]

テンプレート:Mvar は特定の条件を満たす適当な複素数であり、テンプレート:Mvar の値を変えることで収束スピードを調整できる。

超弦理論の研究に必要な計算を簡略化する最適化問題を解く過程で偶然に発見されたものだが、少ない計算量で真値に収束する数学的に優れた式でもある。

この式の背景となったモデルにはベータ関数とファインマン・ダイアグラムが用いられており、量子の相互作用について記述することが可能である。

特殊関数

$ζ (2) = \frac{1}{1^{2}} + \frac{1}{2^{2}} + \frac{1}{3^{2}} + \frac{1}{4^{2}} + \dots = \frac{π^{2}}{6}$ （1735年：オイラー、バーゼル問題、ゼータ関数）
$ζ (4) = \frac{1}{1^{4}} + \frac{1}{2^{4}} + \frac{1}{3^{4}} + \frac{1}{4^{4}} + \dots = \frac{π^{4}}{90}$
$ζ (2 n) = \frac{(- 1)^{n - 1} 2^{2 n - 1} B_{2 n}}{(2 n)!} π^{2 n} (n \in ℕ)$ （テンプレート:Mvar はベルヌーイ数）
$Γ (\frac{1}{2}) = (- \frac{1}{2})! = \sqrt{π}$ （ガンマ関数）

数論

整数全体から無作為に2つ取り出すとき、その2つが互いに素である確率はテンプレート:Math である（互いに素_(整数論)#互いに素である確率を参照）。

力学系・エルゴード理論

ロジスティック写像テンプレート:Math2 により帰納的に定まる数列テンプレート:Math を考える。初期値テンプレート:Math をテンプレート:Math 以上テンプレート:Math 以下に取るとき、そのほとんど全てで、次が成り立つ。

$\lim_{n \to \infty} \frac{1}{n} \sum_{i = 1}^{n} \sqrt{x_{i}} = \frac{2}{π}$

統計

$f (x) = \frac{1}{σ \sqrt{2 π}} \exp [- \frac{(x - μ)^{2}}{2 σ^{2}}]$ （正規分布の確率密度関数）
幅テンプレート:Math の無数の平行線の上から長さテンプレート:Math の針を落とすとき、その針が直線と共有点を持つ確率はテンプレート:Math である（ビュフォンの針）。

その他

河川の長さの水源・河口間の直線距離に対する比率は、平均すると円周率に近い^[24]。

暗唱

語呂合わせ

テンプレート:Main テンプレート:Π の桁を記憶術に頼らずに暗記する方法が各種存在している。

日本語では、語呂合わせにより、長い桁を暗記するのも比較的簡単である。有名なものとして、以下がある。テンプレート:Quote

（	産	医	師	異	国	ニ	向	コー、	産	後	厄	無	ク	産	婦	御	社	ニ	虫	サン	ザン	闇	ニ	鳴	ク	）
	3.	1	4	1	59	2	6	5	3	5	89	7	9	3	2	3	846	2	64	3	3	83	2	7	9	（小数以下30桁）

産	医	師	異	国	に	向	かう。	産	後	厄	な	く、	産	に	産	婆、	四	郎	二	郎	死	産、	産	婆	産	に	泣	く。	困る	に	母	よ	行	く	な。	一	郎	苦	産で	苦が続き、	美	奈	子	一人	を	小	屋	に	置	く。
3.	1	4	1	59	2	6	5	3	5	89	7	9	3	2	3	8	4	6	2	6	4	3	3	8	3	2	7	9	50	2	88	4	1	9	7	1	6	9	3	99	3	7	5	1	0	5	8	2	0	9	（小数以下55桁）^[25]

全く傾向が異なるものとして、

身	一つ、	宵、	獄	に	向	こう。	惨たるかな	医	薬	な	く
3.	1	41	59	2	6	5	3	5	89	7	9		（小数以下14桁）^[26]

身	ひとつ	よ	人の、	いづこに	婿見、	いつ、	厄なく	見つ、	文や	読むらん
3.	1	4	1	592	653	5	8979	3	238	46	（小数以下20桁）^[27]

英語圏では語呂合わせがうまくいかないため、単語の文字数で覚える方法がある。

Yes,	I	have	a	number.
3.	1	4	1	6	（小数以下4桁までで四捨五入）

Can	I	find	a	trick	recalling	Pi	easily?
3.	1	4	1	5	9	2	6	（7桁、また「テンプレート:Π を簡単に思い出せるトリックってある？」という文章自体がその質問の答えにもなっている）

How	I	like	a	drink,	alcoholic	of	course,	after	the	heavy	lectures	involving	quantum	mechanics!	^[28]
3.	1	4	1	5	9	2	6	5	3	5	8	9	7	9	（小数以下14桁）

3桁目の like を want としたものもある（出典は不明）。

And	if	the	lectures	were	boring	or	tiring,	then	any	odd	thinking	was	on	quartic	equations	again.
3	2	3	8	4	6	2	6	4	3	3	8	3	2	7	9	5	（上に続けて、31桁）S. ボトムリー

これらのような覚え方は多くあり、日本語では上記のものの改編で90桁までのものや、歌に合わせたもの、数値を文字に置き換えて1,000桁近く覚える方法などがある。

暗唱記録

2004年 9月25日、原口證が8時間45分かけて円周率5万4000桁の暗唱に成功し、従来の世界記録を更新した。しかしながら、実際はより多くの桁を覚えていたため、2005年 7月1日 - 7月2日に再挑戦し、8万3431桁までの暗唱に成功した。2006年 10月3日午前9時 - 10月4日午前1時30分（16時間30分）の挑戦で円周率10万桁の暗唱に成功したテンプレート:要出典。

2022年2月現在で『ギネス世界記録』に認定されている円周率暗唱の世界記録は、2015年 3月21日にRajveer Meenaが10時間近くかけて暗唱した7万桁である^[29]。

文化的影響

テンプレート:節スタブ

円という日常でもよく知られた図形についての単純な定義でありながら、小数部分が循環せずに無限に続くという不可思議さから、数学における概念の中で最もよく知られたものの一つである。

3月14日は円周率の日および数学の日^{[注 3]}である。小数点以下が「永遠に続く」という意味にあやかり、3月14日に結婚するカップルもいる^[30]。また、テンプレート:Pi (テンプレート:En) とパイ (テンプレート:En) は同音異義語であること^[31]、パイが円形であることから、アメリカ合衆国など複数の国で「パイの日」として祝われ^[32]、パイ焼きやパイ食のほか、数学に関係した活動が行われる^[33]。
7月22日は円周率近似値の日とされている（[[22/7|テンプレート:Math]] は円周率の近似値）。
1999年の学習指導要領の改訂により「小学校の算数で円周率は3で計算することになる」との噂が世間に広まった^[34]が、実際には必要に応じて3で計算することも可能にするための措置であった^[35]。
2012年 8月14日、米国勢調査局が、米国の人口が円周率と同じ並びの3億1415万9265人に達したと発表した。アメリカには円周率の曲を作る人もいる^[36]。
組版処理ソフトウェア [[TeX|テンプレート:TeX]] のバージョン番号は、テンプレート:Math2 というように、更新の度に円周率に近づいていくように一桁ずつ増やされる。

実務上の近似値

円弧の長さの計算など、実務上の数値計算では、その用途に応じて必要な桁数の円周率が計算に用いられる。例として、

テンプレート:要出典範囲
公認陸上競技場の曲走路の計算では、3.1416を用いている^[37]。
NASAのJPLの惑星間航行システムにおける最高精度の計算では、小数点以下15桁までのテンプレート:Valを用いているテンプレート:Sfn。
観測可能な宇宙が球体だとして、その円周の誤差が水素原子程度になるためには、小数点以下40桁程度を使えば足りる^[38]。

値

小数点以下1000桁までの値テンプレート:Π = 3.
1415926535 8979323846 2643383279 5028841971 6939937510 5820974944 5923078164 0628620899 8628034825 3421170679 8214808651 3282306647 0938446095 5058223172 5359408128 4811174502 8410270193 8521105559 6446229489 5493038196 4428810975 6659334461 2847564823 3786783165 2712019091 4564856692 3460348610 4543266482 1339360726 0249141273 7245870066 0631558817 4881520920 9628292540 9171536436 7892590360 0113305305 4882046652 1384146951 9415116094 3305727036 5759591953 0921861173 8193261179 3105118548 0744623799 6274956735 1885752724 8912279381 8301194912 9833673362 4406566430 8602139494 6395224737 1907021798 6094370277 0539217176 2931767523 8467481846 7669405132 0005681271 4526356082 7785771342 7577896091 7363717872 1468440901 2249534301 4654958537 1050792279 6892589235 4201995611 2129021960 8640344181 5981362977 4771309960 5187072113 4999999837 2978049951 0597317328 1609631859 5024459455 3469083026 4252230825 3344685035 2619311881 7101000313 7838752886 5875332083 8142061717 7669147303 5982534904 2875546873 1159562863 8823537875 9375195778 1857780532 1712268066 1300192787 6611195909 2164201989 … （テンプレート:OEIS）^[39]

十進記数法以外の表記法による表現

二進記数法（テンプレート:仮リンク2）による最初の48桁（48ビットとも呼ばれる）はテンプレート:Gaps（テンプレート:OEIS）
三進記数法（基数3）による最初の38桁はテンプレート:Gaps（テンプレート:OEIS）
十六進記数法（基数16）による最初の20桁はテンプレート:Gaps テンプレート:Sfn（テンプレート:OEIS）
六十進記数法（基数60）による最初の5桁は 3;8,29,44,0,47^[40]（テンプレート:OEIS）

注釈

テンプレート:Notelist2

出典

テンプレート:Reflist

参考文献

テンプレート:参照方法

外部リンク

テンプレート:Commonscat テンプレート:Wiktionary

算数用語集「円周率」 - 新興出版社啓林館
円周率.jp
江戸の数学「コラム円周率」 - 国立国会図書館
円周率ナビ
円周率計算｜ラマヌジャン型公式
円周率100万桁表
寒川光：「完全楕円積分とガウス・ルジャンドル法によるπの計算」、平成29年（西暦2017年）4月14日
テンプレート:YouTube time
- ※ 円周率30兆桁をGoogleCloudを使って計算。
Nayandeep Deka Baruah, Bruce C. Berndt and Heng Huat Chan: "Ramanujan's Series for 1/π: A Survey", The American Mathematical Monthly, v116(7) (Aug-Sep.,2009),pp.567-587.

テンプレート:Normdaten

↑ テンプレート:Cite web
↑ テンプレート:Cite book
↑ ^3.0 ^3.1 ^3.2 テンプレート:Cite Kotobank
↑ テンプレート:Cite book
↑ ^5.0 ^5.1 引用エラー: 無効な <ref> タグです。「πの文字使用について」という名前の注釈に対するテキストが指定されていません
↑ テンプレート:Cite book
↑ テンプレート:Cite book
↑ テンプレート:Cite journal
↑ テンプレート:Cite web
↑ "An {ENIAC} Determination of pi and e to more than 2000 Decimal Places", Mathematical Tables and Other Aids to Computation, 4 (29), pp.11-15.（1950年1月）
↑ ^11.0 ^11.1 テンプレート:Cite web
↑ テンプレート:OEIS
↑ テンプレート:Harvnb
↑ テンプレート:Cite OEIS
↑ テンプレート:Cite journal
↑ テンプレート:Cite journal
↑ テンプレート:Cite web
↑ ^18.0 ^18.1 ^18.2 テンプレート:Cite web
↑ テンプレート:Cite web
↑ ^20.0 ^20.1 テンプレート:Cite web
↑ テンプレート:PDFlink
↑ テンプレート:Cite web
↑ テンプレート:Cite web
↑ テンプレート:Cite book
↑ テンプレート:Cite book 小泉が見聞した一番長いものとしている。
↑ テンプレート:Cite book 小泉は「どれも陰惨な文章なのは妙だが、･･･」と書いている。
↑ 難かしい公式も樂に覺えられる算術うた繪本（わかもと物識繪本第2輯） 3ページ写真円周率、1937年4月
↑ IUPAC 物理化学で用いられる量・単位・記号第3版 p.137、5 基礎物理定数、よく使われる数学定数の値、ISBN 978-4-06-154359-1、講談社サイエンティフィク、2009年4月20日第1刷
↑ テンプレート:Cite web
↑ 安田美沙子3・14結婚は『円周率＝永遠』の意味だったテンプレート:Webarchive スポニチアネックス 2014年3月16日（日）12時17分配信
↑ テンプレート:Cite web
↑ テンプレート:Cite web
↑ テンプレート:Cite news
↑ テンプレート:Cite news
↑ テンプレート:Cite news
↑ 米国の人口が円周率と「同じ」に 3億1415万9265人 CNN 2012.08.15 Wed posted at 12:42 JST
↑ 陸上競技場公認に関する細則、競技場に関する規程、細則第3条（距離計測）第1項（5）「曲走路の計算法は、前号の方法によって算出した実長の平均（実測半径という）に300㎜を加えて（計算半径という）円周率（3.1416）を掛けて計算する。」p.403, ルール・ハンドブック、陸上競技ルールブック2022、JAAF 日本陸上競技連盟公式サイト
↑ テンプレート:Cite web
↑ Decimal expansion of Pi (or digits of Pi). Table of n, a(n) for n = 1..20000
↑ テンプレート:Citation. Ptolemy used a three-sexagesimal-digit approximation, and Jamshīd al-Kāshī expanded this to nine digits; see テンプレート:Citation

引用エラー: 「注」という名前のグループの <ref> タグがありますが、対応する <references group="注"/> タグが見つかりません

[1] テンプレート:Cite web

[2] テンプレート:Cite book

[perimetros-3] 3.0 ^3.1 ^3.2 テンプレート:Cite Kotobank

[4] テンプレート:Cite book

[πの文字使用について-5] 5.0 ^5.1 引用エラー: 無効な <ref> タグです。「πの文字使用について」という名前の注釈に対するテキストが指定されていません

[6] テンプレート:Cite book

[8] テンプレート:Cite book

[9] テンプレート:Cite journal

[10] テンプレート:Cite web

[11] "An {ENIAC} Determination of pi and e to more than 2000 Decimal Places", Mathematical Tables and Other Aids to Computation, 4 (29), pp.11-15.（1950年1月）

[:1-12] 11.0 ^11.1 テンプレート:Cite web

[13] テンプレート:OEIS

[14] テンプレート:Harvnb

[15] テンプレート:Cite OEIS

[16] テンプレート:Cite journal

[17] テンプレート:Cite journal

[18] テンプレート:Cite web

[piqt-19] 18.0 ^18.1 ^18.2 テンプレート:Cite web

[20] テンプレート:Cite web

[:0-22] 20.0 ^20.1 テンプレート:Cite web

[23] テンプレート:PDFlink

[24] テンプレート:Cite web

[25] テンプレート:Cite web

[26] テンプレート:Cite book

[27] テンプレート:Cite book 小泉が見聞した一番長いものとしている。

[28] テンプレート:Cite book 小泉は「どれも陰惨な文章なのは妙だが、･･･」と書いている。

[29] 難かしい公式も樂に覺えられる算術うた繪本（わかもと物識繪本第2輯） 3ページ写真円周率、1937年4月

[30] IUPAC 物理化学で用いられる量・単位・記号第3版 p.137、5 基礎物理定数、よく使われる数学定数の値、ISBN 978-4-06-154359-1、講談社サイエンティフィク、2009年4月20日第1刷

[31] テンプレート:Cite web

[33] 安田美沙子3・14結婚は『円周率＝永遠』の意味だったテンプレート:Webarchive スポニチアネックス 2014年3月16日（日）12時17分配信

[34] テンプレート:Cite web

[35] テンプレート:Cite web

[36] テンプレート:Cite news

[37] テンプレート:Cite news

[38] テンプレート:Cite news

[39] 米国の人口が円周率と「同じ」に 3億1415万9265人 CNN 2012.08.15 Wed posted at 12:42 JST

[40] 陸上競技場公認に関する細則、競技場に関する規程、細則第3条（距離計測）第1項（5）「曲走路の計算法は、前号の方法によって算出した実長の平均（実測半径という）に300㎜を加えて（計算半径という）円周率（3.1416）を掛けて計算する。」p.403, ルール・ハンドブック、陸上競技ルールブック2022、JAAF 日本陸上競技連盟公式サイト

[41] テンプレート:Cite web

[42] Decimal expansion of Pi (or digits of Pi). Table of n, a(n) for n = 1..20000

[43] テンプレート:Citation. Ptolemy used a three-sexagesimal-digit approximation, and Jamshīd al-Kāshī expanded this to nine digits; see テンプレート:Citation

[1]

[2]

[3]

[4]

[5]

[6]

[注 1]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16]

[17]

[18]

[19]

[注 2]

[20]

[21]

[22]

[23]

[24]

[25]

[26]

[27]

[28]

[29]

[注 3]

[30]

[31]

[32]

[33]

[34]

[35]

[36]

[37]

[38]

[39]

[40]

円周率

基礎