円周率の歴史

本記事では、数学定数のひとつである円周率の歴史について詳述する。

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

1. 凡例
-----

[学]：数学的事実に関する発見・論争等
[法]：計算法の考案・改良等
[値]：計算・値の使用
[値]：計算・値の使用
[文]：文化・社会

2. 年表
-----

2.1. 級数の発見前　— 13世紀まで —
----------------------

紀元前2000年頃[値]  1936年にスーサで発見された粘土板などから、古代バビロニアでは、正六角形の周と円周を比べ、円周率の近似値として 3}} や 3+sfrac, 3+sfrac などが使われたと考えられている。紀元前1650年頃[学][値] 既に、古代エジプトでは、円周と直径の比の値と、円の面積と半径の平方の比の値が等しいことは知られていた。神官アハメスが書き残したリンド・パピルスには、円積問題の古典的な解法の一つが記されており、円の直径からその sfrac を引いた長さを一辺とする正方形の面積と、元の円の面積が等しいとしている。これは、円周率を近似的にsfracとみなすことに相当し、それなりに精度の高い近似値であったが普及はしなかった。リンド・パピルスはアハメスによって写されたものであり、内容自体はさらに紀元前1800年頃にまで遡ると考えられている。紀元前5世紀頃[学] アナクサゴラスが、アポロンへの不敬罪で投獄されている間に、円積問題に取り組んだ。[法] ヘラクレアのアンティフォンは、円に内接する正多角形の面積を求めることにより円周率を計算する方法を編み出した。アンティフォンは、それぞれの正多角形から正方形が作図できることから、円積問題が解決できると主張した。[値] すぐに、同じヘラクレアのブリソンが、外接する正多角形の面積を求めて内側と外側の両方から円の面積を評価し近似値を得た。紀元前3世紀[法][値] アルキメデスは、円の面積が円周率と半径の平方の積に等しいことを証明した。さらに、3の平方根の最良近似分数 sfrac および sfrac  を利用して、円に外接および内接する正六角形、正十二角形、正二十四角形、正四十八角形、正九十六角形の辺の長さの上界および下界をそれぞれ計算することにより 3 + sfrac を求めた。小数だと 3.14084 < //π// < 3.14286}} である。1世紀[値] ローマ帝国の著名な建築家ウィトルウィウスは、sfrac を使った。素数の7}}よりも、2}}の3}}乗である 8}} で割ったほうが建築には便利だったためである。小数だと 3.125}} である。2世紀[値] 天文学者プトレマイオスは sfrac
を使った。小数だと約3.1417}} である。[値] 後漢の太史令だった張衡は、円に外接する正方形の周と円周を比べ、円周率を sqrt とした。約3.162}} になる。3世紀[値] 呉の王蕃は sfrac を用いた。約3.1555}} である。263年[値]  魏の劉徽は『九章算術』の注釈の中で、ブリソンと同様の方法を用い 3.14 + sfrac であることを示している。小数では 3.14102 4 < //π// < 3.14270 4}} である。さらに正3072}}角形を用いて、3.14159}} という近似値も得た。5世紀[値]  7世紀に編纂された隋書律暦志によると、天文学者の祖沖之は、当時としては非常に正確な評価 3.14159 26 < π < 3.14159 27}} を示した。ヨーロッパでこれほど正確な評価を得るには、16世紀まで待たねばならない。さらに、分数での近似値 sfracと sfracを与えている。正確な方法は伝わっていないが、九章算術の方法を踏襲したと推測すると、上記の結果を得るには少なくとも円に内接する正24576}}角形の辺の長さを計算しなければならない。隋書では現代と同じ「圓周率」という語が用いられている。祖沖之の息子のzhは、父とともに球の体積の計算方法を導き出したことで知られる。500年頃[値] インドのアリヤバータは、円に内接する正 n}} 角形と正 2//n//}} 角形の周の長さの間に成り立つ関係式を求め、正384}}角形の周の長さから sqrt  と求めた。この平方根の近似値として sfrac  を与えた。650年頃[値] インドのブラーマグプタは、正12}}角形、正24}}角形、正48}}角形、正96}}角形の周の長さから、n}} が大きくなるにつれ正 3 × 2sup 角形の周の長さは sqrt に近づくとし、これを円周率とした。1220年[値] イタリアのレオナルド・フィボナッチが円周率を sfrac と計算した。これは、約3.1418}} である。

2.2. 幾何から解析へ　— 14世紀から20世紀前半 —
-----------------------------

「円に内接・外接する多角形に基づく近似」なる幾何学的な考察から「級数を利用した近似」なる解析的な考察への移行は、インドでは1400年頃から1500年代に、ヨーロッパでは1600年代に、日本では1700年代に起きた。

14世紀1400年頃インド南西部では天文学・数学が花開き、当時の世界最先端の研究が行われた。ケーララ学派と総称される学者たちは、三角関数・逆三角関数のマクローリン展開を天文計算に利用した。これらの級数はニーラカンタの時代には既に知られており、ニーラカンタの発見とされることがある。しかし、ニーラカンタの天文学書『アールヤバティーヤ・バーシャ』によると、正弦のマクローリン級数の展開式は彼より前の時代の学者の業績であるという。その学者とは、サンガマグラーマのマーダヴァである。以下の式も、マーダヴァの発見とされることが多い：この級数展開は、以下と同等である。この級数は、しばらくしてヨーロッパでもジェームズ・グレゴリーとゴットフリート・ライプニッツにより再発見され、一般的にはグレゴリー級数、もしくはライプニッツ級数などと呼ばれる。[値]  一説によると、マーダヴァは上式から直ちに得られる等式、の21項を計算し、//π// ≈ 3.14159 26535 9}} を得たという。これは小数点以下10桁目まで正しい。別の資料によると、彼の近似値は //π// ≈ sfrac で、これは //π// ≈ 3.14159 26535 92222…}} に当たる。マーダヴァが円周率10桁を得たとすると、祖沖之の7桁以来、約1000年ぶりの世界記録更新である。上記の級数を30項目まで使えば円周率の15桁が決定でき、42項目まで使えば20桁が決定できる。この他にもケーララ学派は円周率の評価に利用できるいくつもの結果を得ていて、その気になれば比較的簡単に円周率の桁数を伸ばせる立場にあった。実際、R. Gupta は、マーダヴァが約17桁まで計算したと予想している。しかし、記録は見つかっておらず、現時点では想像の域を出ない。[法] ケーララ学派による円周率の近似は級数に基づくもので、剰余項も考察している。他地域ではこの200年後にまだ正多角形の外周に基づく計算をしていることを考えると、極めて先進的だった。円周率の計算法として新しいというだけでなく、無限や極限を扱う新しい数学への大きな一歩だった。 15世紀1424年[値] 
ペルシャの天文学者・数学者ジャムシード・カーシャーニーは、当時使われていた円周率の近似値の不正確さに不満を抱き、天文計算に必要十分な精度で円周と半径の比を決定したいと考えた。1424年の『円周論』において、彼はアルキメデスの方法を拡張して正805,306,368}}  角形を用いる計算を行い、60進数による次の評価を得た。ここで、6; 16, 59, …}} は 6 + sfrac を表す。現代的な表記に直せば:彼は近似値 2//π// =}} 6; 16, 59, 28, 1, 34, 51, 46, 14, 50}} を採用し、10進表示 //π// =}} 3.14159 26535 89793 25}} も与えた。これは小数点以下16桁目まで正しく、末尾の17桁目も真の値に近い。記録に残る当時最良の円周率の近似値であり、この世界記録は1596年にルドルフ・ファン・コーレンが小数点以下20桁を示すまで172年間、破られなかった。この業績は、西洋では1920年代まで知られていなかった。1500年頃[学] ケーララの天文学者ニーラカンタが、円周率の無理性を指摘した。彼の著書『アールヤバティーヤ・バーシャ』には、こうある：「直径が何らかの長さの単位で計測されて、その単位の比として表されるなら、その同じ単位によって円周を同様に計測することはできない。よってまた同様に、円周が何らかの単位で計測可能であるのなら、直径はその同じ単位によっては計測できない。」ケーララ学派は円周率の級数表示を知っていたため、この認識は自然に生じたのだろう。[値]  ニーラカンタの『タントラ・サングラハ』には、エレガントな分数表示 //π// ≈ sfrac が含まれる。これは sfrac と同様の最良近似分数で、小数点以下9桁目まで正しい。 16世紀1503年アルキメデスの『円の計測について』と『放物線の面積について』のラテン語訳が、ベネチアで出版された。1543年ニコロ・フォンタナ・タルタリア  が、アルキメデスの一部の著作のラテン語訳をベネチアで再出版した。1544年アルキメデスの著作の原文が、初めてまとめて出版された。出版地はバーゼル
で、ラテン語訳付きだった。これによりヨーロッパでは彼の業績が広く知られるようになり、円周率の研究もこれを出発点として本格的に再開された。この時点での西洋の円周率研究は紀元前のアルキメデスの時代からあまり進歩していなかったが、これ以降は急速に発展する。1579年[値]  フランソワ・ビエタが、円に内接・外接する正393,216}}角形の周の長さから 3.14159 26535 < //π// < 3.14159 26537}} という評価をした。ビエタはさらに、無限乗積を示し の計算を試みた。1585年[値] オランダのアドリアン・アンソニスが sfrac と評価し、両端の平均に近い値として sfrac を得た。これは、約3.14159 292}} である。1593年[値]  フランドルのアドリアーン・ファン・ローメンが、『数学的観念序説：多角形法』の中で 3.14159 26535 89793 05 < //π// < 3.14159 26535 89793 15}} に当たる評価を与え、//π// ≈ 3.14159 26535 89793 1}} とした。これは小数点以下15桁目まで正しい。アル＝カーシーの世界記録16桁  にはわずかに及ばなかったが、この時点でヨーロッパ最良の近似値であり、ビエトの結果  の改良となっている。ただし、円周率の真の値は上記の区間に含まれておらず、厳密な評価ではない。計算は正 15 × 2sup角形を用いるものだった。彼は21歳年上のファン・コーレンと親交があり、円周率に興味を持ち始めたのは彼の影響らしい。1596年[値]  ルドルフ・ファン・コーレンが、『円について』で円周率の小数点以下20桁を決定した。ファン・コーレンはまず、正 5 × 2｛sup角形、正 4 × 2sup角形、正 3 × 2sup角形を用いて、円周率をそれぞれ12桁、16桁、18桁まで求めた。さらに、正 15 × 2sup  角形に基づき次の評価を与えた:上界・下界の平均を取って //π// ≈ 3.14159 26535 89793 23846}}
とすれば、結果的に全20桁が正しい。しかし、ファン・コーレンの態度は厳格で、上記の結果は19桁のみ有効であると正しく指摘した。最後に彼は の20桁を示した:この計算は、辺の数をさらに2倍にした正 15 × 2sup  角形に基づく。ファン・ローメンの15桁の計算  の改良であり、アル＝カーシーの16桁の記録  を上回る新しい世界記録の達成だった。ファン・コーレンはヒルデスハイムで生まれ、ホラントに移住した。フェンシングと数学の教師だった。高等教育は受けていなかったが、円積問題や円周率を巡る数学上の論争に巻き込まれ、1590年頃から円周率に興味を持ち始めたらしい。 17世紀1610年頃[値]  ファン・コーレンは、1610年に亡くなるまでのいずれかの時点で、正 2sup角形を使って の35桁目までを正しく評価した。この結果は、1621年、弟子のスネリウスの著書『キュクロメトリクス：円の計測について』で公表されたほか、本人の墓に刻まれた。墓石は後代に滅失したが、碑文とスケッチは残っており、2000年に復元された。かつてドイツでは、彼の名に因んで円周率をルドルフ数  と呼んだ。1621年[法][値] オランダのヴィレブロルト・スネルが、円周の長さの評価式を与える。この式と円に内接・外接する正 6}} 角形から 3.14022 < //π// < 3.14160}} と評価した。この式の証明はクリスティアーン・ホイヘンスによって与えられ、さらにホイヘンスによって改良された結果、正六角形を用いただけで 3.14159 26533 < //π// < 3.14159 26538}} と評価できるまでになった。スネリウスはファン・コーレンの弟子だった。彼の方法なら、ファン・コーレンが正 2sup 角形を使って得た35桁は、2sup 角形を考えるだけで得られるという。その気になれば、計算記録を更新できる立場だった。しかし、彼は別の分野で活躍しており、すでに35桁あった円周率の有効数字をさらに伸ばすために時間を割くことはしなかった。1630年[値] 
オーストリア出身の天文学者・数学者クリストフ・グリーンベルガーは、スネリウスの手法を用いて円周率の小数点以下39桁目までを計算し、1630年に出版された自著『三角法の基礎』の中で公表した。39桁目は 7 だが、彼はそれを 6 と 9 の間だと正しく評価した。桁数という意味では38桁目まで確定させたことになる。1655年[法] イギリスのジョン・ウォリスは無限乗積を示した。ビエタの公式のように根号が無いため計算はしやすいが、収束はとても遅い。同じくイギリスのブラウンカーが、連分数を用いた公式を示した。この公式により が無理数であることが分かる。1663年[値] 村松茂清が『算俎』を著し、円に内接する正 2sup 角形 }} の辺の長さから //π// ≒ 3.1415 92648 77769 88692 48}} とし、小数点以下7桁まで正しい値を求めた。ファン・コーレンなどの計算には遠く及ばないものの、近似値として単に 3.16}} という値を示すのみであった『塵劫記』や、中国などを通じて入ってくる算書に頼り切ってきたそれ以前の和算から一歩を踏み出し、日本で初めて数学的な方法で円周率を計算し発表した和算家が村松である。和算において、円周率をはじめとする円に関する研究は「円理」と呼ばれ、一定の発展を見せたが、例えば外接多角形との「はさみうち」によって何桁目まで正しいかを論証する、といったような基本的な数学的発展さえわずかであったのが「和算の限界」であった。1665年[学] イギリスの政治哲学者のトマス・ホッブズが円積問題の解を公表し、ウォリスとの間で論争になる。ホッブズは死ぬまで厳密解と近似解の違いを理解できずに論争を続けた。1671年[法] スコットランドのジェームス・グレゴリーにより、グレゴリー級数が発見される。これとは独立に1674年にゴットフリート・ライプニッツも同じ発見をしており、グレゴリー・ライプニッツ級数とも呼ばれる。ライプニッツは //x// =}} 1}} を代入し、マーダヴァと同じ級数を得た。1681年[値] 暦の作成にあたって円周率の近似値が必要になったため、関孝和が正 131,072}} 角形を使って小数第 16
位まで算出した。関が最終的に採用した近似値は「3.14159 26535 9微弱」というものだったが、エイトケン補外を用いた途中計算では小数第 16 位まで正確に求めている。西洋でエイトケン補外が再発見されたのは1876年、ハンス・フォン・ネーゲルスバッハによってである。1699年[値]  イギリス人のエイブラハム・シャープがグレゴリー・ライプニッツ級数に を入れ、 を小数第 72 位まで求めた。

18世紀1706年[法][値]  イギリスのジョン・マチンがマチンの公式を発見する。さらに、この関係式にグレゴリー・ライプニッツ級数を用いて小数第 100 位までの円周率を求めた。[文] ウィリアム・ジョーンズが初めて を円周率の意味で用いた。1748年にレオンハルト・オイラーも同じ記法を用いたことで円周率を と表記することが広まった。1719年[値]  フランスのトーマス・ラグニーが、シャープの方法で小数第 127 位まで計算を行う。1722年[値] 建部賢弘が『綴術算経』を著し、正 1024}} 角形を用いて小数第 42 位まで求めた。「累遍増約術」を適用し、関孝和の計算に比べて遥かに少ない計算で精度を大いに改善している。なお、ルイス・フライ・リチャードソンによる同手法の提案は1910年頃である。1761年[学] ドイツのヨハン・ハインリッヒ・ランベルトによって が有理数でないことが証明される。18世紀中頃[法] レオンハルト・オイラーによって、多くの に関する式が発見される。オイラーはを用いて、 たった1時間で円周率を小数第 20 位まで計算した。1775年[学] フランスの科学アカデミーが、ギリシアの三大作図問題と永久機関についての論文審査を拒否する決議をした。1789年[値]  スロベニアの数学者ユーリイ・ヴェガは、マチンの公式を用いて小数第 140 位まで値を求め、小数第 137 位までが正しかった。この記録はその後50年破られることがなかった。1794年[学] アドリアン＝マリ・ルジャンドルによって は有理数の平方根にならないことが証明される。 19世紀 19世紀初期 [法] カール・フリードリヒ・ガウスとアドリアン＝マリ・ルジャンドルが独立に、算術平均と幾何平均を利用した反復計算アルゴリズムを研究。1970年代に再発見され、現在ではガウス＝ルジャンドルのアルゴリズムと呼ばれる。円周率を計算するものの中では非常に収束が速い。1850年頃 - 1873年[値]  イギリスのウィリアム・ラザフォードとその弟子のウィリアム・シャンクスがマチンの公式を用いて桁数の記録を塗り替えた。1852年にラザフォードが小数第 441
位、シャンクスが小数第 530 位まで計算し、小数第 441 位までは両者の計算が一致していることでその計算の正しさを確認できた。しかし、arctan sfrac が小数第 530 位までしか正しくなく、シャンクスの計算で正しかったのは、小数第 527 位までであった。その後、シャンクスは1872年に小数第 707 位まで達したが、この誤りが最後までつきまとった。1882年[学] フェルディナント・フォン・リンデマンによって が代数的数でないことが証明される。これにより の超越性が証明され、円積問題も否定的に解決された。1896年[法] カール・ストーマーは公式を発見する。1897年[文][値] アメリカ合衆国のインディアナ州の下院で、医者のエドウィン・グッドウィンによる円積問題解決方法を盛り込んだ議案246号が満場一致で通過した。グッドウィンの方法から得られる値は π =}} 3.1604, 3.2, 3.232, 4}} であり、このうち 4 については、公式に認められた最も不正確な円周率の値としてギネスブックに記載された。この法案は各審議会を通過していき上院に承認を求める段階にまで達した。しかし世論の批判に遭い、2月12日に上院によって議論の無期限延期が決められ、法案成立目前で却下された。 20世紀1910年[法] ラマヌジャンによって、級数表示が発見される。この公式は、ジョナサン & ピーター・ボールウェイン兄弟によって1987年に厳密に証明されるが、1985年にウィリアム・ゴスパーがこの公式を用いて円周率を計算し、その正確さを示している。1945年[値]  ファーガソン  が小数第 540 位までを計算し、ウィリアム・シャンクスの誤りを指摘する。シャンクスの計算は約70年間も信用されていた。このファーガソンの計算までが手計算によるものだった。手計算の時代は誤りが起こることも多かったが、この時代の数学の成果は、現代の計算機による円周率の計算においても非常に重要な役割を果たすことになる。

2.3. 計算機による計算の時代　— 20世紀後半以後 —
-----------------------------

1947–1948年[値]  ファーガソンは、卓上計算機を使用して808桁まで求めた。この計算は、レビ・スミスとジョン・レンチによっても検算され、シャンクスの計算が間違いであることが繰り返し確認された。1949年[値]  ライトウィーズナーが ENIAC を用いてマチンの公式により 2037桁を 70時間かけて計算した。1954年[値]  S・C・ニコルソンとJ・ジーネルが、IBM NORC を用いて3089桁を13分で計算した。1958年[値]フランソワ・ジェニューイが、IBM 704 を用いて 1万桁まで計算した。1961年[値]ジョン・レンチとダニエル・シャンクスが IBM 7090 を用いて 10万桁まで計算した。計算にはストーマーの公式を使用した。検算にはガウスの公式を使用した。1966年[値]パリの原子力エネルギー委員会にある IBM 7030 を用いて25万桁まで計算した。1967年[値]パリの原子力エネルギー委員会にある CDC 6600 を用いて50万桁まで計算した。1973年[値]ジャン・ギューとマルティーヌ・ブイエが CDC 7600 を用いて 100万1250桁まで計算した。[法] ユージン・サラミンとリチャード・ブレントが独立に、算術幾何平均を用いたアルゴリズムを発見する。現在ではガウス＝ルジャンドルのアルゴリズムと呼ばれる。1982年[値]田村良明が MELCOM COSMO 900 Ⅱ を用いてサラミンとブレントのアルゴリズムにより209万7144桁まで計算。[値]田村良明と金田康正が HITAC M-280H を用いて419万4288桁、ついで838万8576桁まで計算。1983年[値]金田康正と吉野さやかが HITAC M-280H を用いて1677万7206桁まで計算。[値]後保範と金田康正が HITAC S-810/20 を用いて1001万3395桁まで計算。アルゴリズムはガウスの公式による。[値]若松登志樹がシャープのパソコン MZ-80B を用いてガウスの公式により7万1508桁まで計算。