円周率の近似

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

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

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

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

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

\pi=3.125

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

\pi=3.16049\cdots

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

### アルキメデスの手法
紀元前3世紀に入ると、円周率の真の値に限りなく近づくことができる画期的な方法を考え出した人が現れた。古代ギリシャの数学者で物理学者のアルキメデス。彼の考えだした方法は次の通り。
 円に内接する正六角形の外周(=3)\pi>3.05を示せ。」という問題も正8角形や正12角形を用いて示すことができる。また、周でなく面積でも求めることができるが、こちらの方が精度が悪く、面積を使って\pi>3を示すためには正12角形を用いる必要がある。\pi>3.05を示すためには正24角形などを用いる必要がある。

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

### ヴィエトの公式
\frac{2}{\pi} = \prod_{n=1}^{\infty} \cos \frac{90^\circ}{2^n} = \sqrt{\frac{1}{2}}\sqrt{\frac{1}{2}+\sqrt{\frac{1}{2}}}\sqrt{\frac{1}{2}+\sqrt{\frac{1}{2}+\sqrt{\frac{1}{2}}}}\cdots

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

### ウォリスの公式
\frac{\pi}{2}=\frac{2\times2\times4\times4\times6\times6\cdots}{1\times3\times3\times5\times5\times7\cdots}

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

### ブラウンカーの公式
### <math>\pi={\cfrac  {4}{1+{\cfrac  {1^{2}}{3+{\cfrac  {2^2 }{5 +{\cfrac  {\cdots}{\cdots+{\cfrac  {n^{2} }{(2m+1)+\cdots}} }}}}}}}}</math>
ウォリスの公式を変形することで導出できる。

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

### 1.arctanを使った公式

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

\arctan{x}=\sum^\infty_{k=0}\frac{x^{2k+1}} {2k+1}=x-\frac{x^3} {3}+\frac{x^5} {5}-\cdots

にx=1を代入し、

\arctan{1}=\frac{\pi}{4}=\frac{1}{1}-\frac{1}{3}+\frac{1}{5}-\frac{1}{7}+\frac{1}{9}-\frac{1}{11}+\frac{1}{13}\cdots

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

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

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

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

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

#### クリンゲンシュティルナ(Klingenstierna)の公式
\frac{\pi}{4}=8\arctan{\frac{1}{10}}-\arctan{\frac{1}{239}}-4\arctan{\frac{1}{515}}

#### ガウス(Gauss)の公式
\frac{\pi}{4}=12\arctan{\frac{1}{18}}+8\arctan{\frac{1}{57}}-\arctan{\frac{1}{239}}

#### 高野喜久雄の公式
\frac{\pi}{4}=12\arctan{\frac{1}{49}}+32\arctan{\frac{1}{57}}-5\arctan{\frac{1}{239}}+12\arctan{\frac{1}{110443}}

#### オイラー(Euler)の公式
\frac{\pi}{4}=\arctan{\frac{1}{2}}+\arctan{\frac{1}{3}}

\frac{\pi}{4}=5\arctan{\frac{1}{7}}+2\arctan{\frac{3}{79}}

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

#### ストーマー(Störmer)の公式
\frac{\pi}{4}=6\arctan{\frac{1}{8}}+2\arctan{\frac{1}{57}}+\arctan{\frac{1}{239}}

\frac{\pi}{4}=44\arctan{\frac{1}{57}}+7\arctan{\frac{1}{239}}-12\arctan{\frac{1}{682}}+24\arctan{\frac{1}{12943}}

#### ヴェガ(Vega)の公式
\frac{\pi}{4}=2\arctan{\frac{1}{2}}-\arctan{\frac{1}{7}}

\frac{\pi}{4}=4\arctan{\frac{1}{5}}-2\arctan{\frac{1}{408}}+\arctan{\frac{1}{1393}}

#### クラウゼン(Clausen)の公式
\frac{\pi}{4}=2\arctan{\frac{1}{3}}+\arctan{\frac{1}{7}}

#### ダース(Dahse)の公式
\frac{\pi}{4}=\arctan{\frac{1}{2}}+\arctan{\frac{1}{5}}+\arctan{\frac{1}{8}}

#### ラザフォード(Rutherford)の公式
\frac{\pi}{4}=4\arctan{\frac{1}{5}}-\arctan{\frac{1}{70}}+\arctan{\frac{1}{99}}

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

#### ラマヌジャン(Ramanujan)の公式
\frac{1}{\pi}=\frac{2\sqrt{2}}{99^2}\sum_{n=0}^\infty \frac{(4n)!(1103+26390n)}{(4^n99^nn!)^4}

\frac{4}{\pi}=\sum_{n=0}^\infty \frac{{(-1)}^n(21460n+1123)(4n!)}

より

2\sqrt{2}

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

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

## 出典



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