円周率の近似

本記事では、数学定数のひとつである円周率の近似（えんしゅうりつのきんじ）について詳述する。

円周率 π は無理数であるため、小数部分は循環せず無限に続く。さらに、円周率 π は超越数でもあるため、その連分数表示は循環しない。その近似値は何千年にも亘り世界中で計算されてきた。

 紀元前5000年ごろ～紀元前3世紀ごろまで 

 古代の円周率 
人類は少なくとも紀元前5000年ごろから、車輪のような円形の物体を使って重いものを運ぶなど、円を役立つ形として認識し、利用してきた。特に円に関して高い数学的知識を持っていたのが、現在のイラク南部に当たる地域に紀元前2500年ごろから住んでいたバビロニア人。彼らは、円の研究を進め、円周の長さが円の直径に比例すること、すなわち円周率の存在に気付いていた。

 古代エジプトでの円周率 
古代エジプト人は円周率を3.125とした。彼らの結論にたどり着くには、簡単な実験をするだけでよい。まず、ひもと木の棒2本を用意し、ひもの両端に木の棒をくくりつける。そして、コンパスのように円を描く。そのあと、そのままひもを定規の代わりにし、半径の何倍になる調べる。そうすると約6.25倍になることがわかる。あとは2で割ると直径の約3.125倍になることがわかる。

 古代バビロニアでの円周率 
古代バビロニア人は円周率を3.16049...とした。彼らの結論にたどり着く方法は、彼らが残した記録を見ることでわかる。

 紀元前3世紀～16世紀前半まで 

 アルキメデスの手法 
紀元前3世紀に入ると、円周率の真の値に限りなく近づくことができる画期的な方法を考え出した人が現れた。古代ギリシャの数学者で物理学者のアルキメデス。彼の考えだした方法は次の通り。
 円に内接する正六角形の外周(=3)<直径1の円の外周(=π)<円に外接する正六角形の外周(=3.4661...)　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　この円に内接する正多角形と円に外接する正多角形の辺の数を無限に増やせば円周率が求められる。
+アルキメデスの手法で求めた円周率の記録人物時代記録アルキメデス紀元前287年～212年小数点以下2桁祖 沖之430年～501年小数点以下7桁ルドルフ・ファン・ケーレン1540年～1610年小数点以下35桁関 孝和1642年～1708年小数点以下10桁建部 賢弘1664年～1739年小数点以下40桁
アルキメデスの方法では正96角形を使って3.14までを決定しているがこれを見るとかなり精度が悪いことがわかる。また、東京大学の「を示せ。」という問題も正8角形や正12角形を用いて示すことができる。また、周でなく面積でも求めることができるが、こちらの方が精度が悪く、面積を使ってを示すためには正12角形を用いる必要がある。を示すためには正24角形などを用いる必要がある。

 16世紀なかば～16世紀末まで 

 ヴィエトの公式 

フランソワ・ビエトによる。アルキメデスの方法と本質的に同じで三角関数の半角の公式を用いて導出することができる。

 ウォリスの公式 

正弦関数の無限乗積展開を用いて導出することができる。

 ブラウンカーの公式 

ウォリスの公式を変形することで導出できる。

 17世紀～現代 
17世紀からは無限級数(無限和)を使った円周率を求める方法が考え出された。

 1.arctanを使った公式 

 マーダヴァ・グレゴリー・ライプニッツ(Madhava-Gregory-Leibniz)級数 
ライプニッツはarctanのマクローリン展開

にを代入し、

を得た(ライプニッツの公式)。

arctanのマクローリン展開はすでにマーダヴァが、ライプニッツの公式は独立にグレゴリーが発見していたため、マーダヴァ・グレゴリー・ライプニッツ級数と呼ばれるようになった。

 マチン(Machin)の公式 
ジョン・マチンはarctanのマクローリン展開から

を得た。この公式は収束が速く、マチン自身も円周率を100桁求めている。また、ウィリアム・シャンクスは707桁目まで求めたが、その後が527桁目までしか一致しておらず、正しかったのは527桁目までということがわかった。

 クリンゲンシュティルナ(Klingenstierna)の公式 

 ガウス(Gauss)の公式 

 高野喜久雄の公式 

 オイラー(Euler)の公式 

2つ目の公式はオイラー変換と呼ばれる変換でarctanを無限乗積展開するとやが現れるため、十進数での計算がしやすい。実際、オイラー自身もこの公式を使って1時間程度で20桁求めたといわれている。

 ストーマー(Störmer)の公式 

 ヴェガ(Vega)の公式 

 クラウゼン(Clausen)の公式 

 ダース(Dahse)の公式 

 ラザフォード(Rutherford)の公式 

 2.ラマヌジャン型公式…複雑だが、収束が速い公式 

 ラマヌジャン(Ramanujan)の公式 

この級数は1項ごとに約8桁ずつ正確な桁が増える。

 チュドノフスキー(Chudnovsky)級数 

この級数は1項ごとに約14桁ずつ正確な桁が増える。

 3.反復公式 

 ガウス=ルジャンドルのアルゴリズム 
「ガウス＝ルジャンドルのアルゴリズム」も参照

はつの数値の算術幾何平均を求めるために、それぞれの数値を算術平均(相加平均)と幾何平均(相乗平均)で置き換えていくものである。

 初期値の設定 


 反復式 
 が希望する精度（桁数）になるまで以下の計算を繰り返す。小数第位まで求めるとき　 回程度の反復でよい。



  の算出 
円周率  は、 を用いて以下のように近似される。



最初の3回の反復で得られる数値（最後の桁は真値とは異なる）は以下の通りである。

 　（小数点以下桁目までが正しい）
 　（小数点以下桁目までが正しい）
 （小数点以下桁目までが正しい）

この計算過程は二次収束する。ガウス自身もこの式を用いて反復を回まで行って桁まで正しいことを確認したことが知られている。

 その他の近似方法 

を用いる近似方法もある。例えば

より

がわかる。これで四角形の周を利用したものと同じ評価が得られた。

 記録 
Google社の岩尾エマはるかが2019年3月14日に円周率を31兆4159億2653万5897桁計算したと発表した。所要時間は約120日。岩尾エマはるかはその後100兆桁計算したと発表。さらに2024年3月、StorageReviewの編集者らが105兆桁まで計算した。所要時間は約70日。いずれもChudnovskyの公式が使われている。

 関連項目 
円周率が22/7より小さいことの証明
円周率の日

Category:数学に関する記事
Category:近似アルゴリズム
Category:円周率
Category:超越数