マチンの公式のソースを表示

[[file:Machin Formula.svg|thumb|right|マチンの公式の概念図。逆正接関数 {{math|arctan ''x''}} は偏角として考えることができるため、マチンの公式は上図のように解釈することができる。]]
{{出典の明記|date=2015年12月}}
'''マチンの公式'''（{{Lang-en-short|Machin's formula}}）とは、[[1706年]]にイギリスの[[天文学者]][[ジョン・マチン]]によって発見された[[逆三角関数|逆正接関数]] {{math|arctan ''x''}} を用いた[[円周率]]を計算するための公式、すなわち
:<math> \dfrac\pi4 = 4\arctan\dfrac15 -\arctan\dfrac{1}{239} </math>
なる公式である。

== 概要 ==
[[ジェームス・グレゴリー#グレゴリー級数|グレゴリー級数]]すなわち逆正接関数 {{math|arctan ''x''}} の[[テイラー展開|マクローリン展開]]：
:<math>\arctan x=\sum_{n=0}^{\infin} \frac{(-1)^n}{2n+1} \, x^{2n+1} = x - \frac13\,x^3 + \frac15\,x^5 -\frac17\,x^7 + \frac19\,x^9 - \cdots \qquad{\rm where}\quad \vert x \vert \leq 1</math> 
に {{math|''x'' {{=}} 1}} を代入して得られる級数：
:<math>\frac{\pi}{4} =\sum_{n=0}^{\infin} \frac{(-1)^n}{2n+1} =1-\frac{1}{3} +\frac{1}{5} -\frac{1}{7} +\frac{1}{9} -\cdots</math> 
は[[ライプニッツの公式]]と呼ばれ、見た目は綺麗な公式であるものの、収束が非常に遅いことで知られる。しかしながら、{{math|''x''}} を十分小さく取れば見た目の綺麗さは多少損なわれるが、それなりに速く収束する級数を得ることができる。実際、[[エイブラハム・シャープ]]は {{math|''x'' {{=}} 1/{{sqrt|3}}}} を用い、[[円周率]]を小数点以下71桁まで計算した。

[[ジョン・マチン]]は、さらに収束性をよくするために[[三角関数|逆正接関数]] {{math|arctan ''x''}} の関係式を考え、これを[[ジェームス・グレゴリー#グレゴリー級数|グレゴリー級数]]と結び付けることにより、非常に収束速度が速い級数を得た。さらに、この公式を発見したマチン自身も円周率を100桁まで求めることに成功した。マチンの公式や、似たような {{math|arctan ''x''}} を用いた公式は、[[1970年代]]に[[算術幾何平均]]などが用いられるようになるまでは円周率の計算に用いられ計算競争に貢献した。その後しばらくは新しい[[アルゴリズム]]による円周率の計算が続いたが、[[2002年]]に[[金田康正]]によって[[高野喜久雄]]の公式が用いられ円周率を1兆2411億桁まで計算するという記録に結び付いた。

== 公式 ==
等式：
:<math> \dfrac\pi4 = 4\arctan\frac15 - \arctan\frac{1}{239}</math>
を'''マチンの公式'''という。

この項目では {{math|arctan ''x''}} は[[主値]]
:<math>-\frac{\pi}{2} <\arctan x<\frac{\pi}{2}</math>
を取るものとする。

同じことであるが、[[逆三角関数|逆余接関数]] {{math|arccot ''x''}} を用いて、
:<math> \dfrac\pi4 = 4\arccot 5 - \arccot 239 </math>
と書かれることもある。

== 主な証明 ==
=== 三角関数の公式による証明 ===
マチンの公式は三角関数の公式を用いて証明できる。

二倍角公式を2回用いて、
:<math>\tan \left( 2\arctan \frac{1}{5} \right) = \frac{2\tan (\arctan \frac{1}{5})}{1-\tan^2 (\arctan \frac{1}{5})} = \frac{2  \cdot \frac{1}{5}}{1 - (\frac{1}{5})^2} =\frac{5}{12}</math>
:<math>\tan \left( 4\arctan \frac{1}{5} \right) = \tan \left( 2 \cdot 2\arctan \frac{1}{5} \right) = \frac{2\tan (2\arctan \frac{1}{5})}{1-\tan^2 (2\arctan \frac{1}{5})} = \frac{2 \cdot \frac{5}{12}}{1- (\frac{5}{12})^2} = \frac{120}{119}</math>
加法定理により、
:<math>\tan \left( 4\arctan \frac{1}{5}-\frac{\pi}{4} \right) = \frac{\tan (4\arctan \frac{1}{5}) - \tan \frac{\pi}{4}}{1 + \tan (4\arctan \frac{1}{5}) \tan \frac{\pi}{4}} = \frac{\frac{120}{119}-1}{1 + 1 \cdot \frac{120}{119}} = \frac{1}{239}</math>
逆関数をとって、
:<math>\arctan \frac{1}{239} = 4\arctan \frac{1}{5} - \frac{\pi}{4}</math>
:ただし
<math>- \frac{\pi}{2}<4\arctan \frac{1}{5} - \frac{\pi}{4}<\frac{\pi}{2}</math>
が前提ではある。

したがって、
:<math>4\arctan \frac{1}{5} -\arctan \frac{1}{239} = \frac{\pi}{4}</math>

=== 複素数を用いた証明 ===
[[file:Arctan x on Complex Plain.svg|thumb|right|複素平面上における、逆正接関数による偏角の表現。]]
{{math|''a'' }}、{{math|''b''}} を[[実数]]とし、{{math|''a''  > 0 }}、{{math|''i''}} を[[虚数単位]]とする。

[[複素数]] {{math|''z'' {{=}} ''a'' + ''bi''}} の偏角は、
:<math>\arg z = \arctan\frac{b}{a}</math>
である。
:複素数の偏角の範囲は {{math|arctan}} の主値と同じ範囲に取るものとする。
{{math|''n''}} を[[整数]]とする。[[ド・モアブルの定理]]によると{{math| ''z''{{sup|''n''}} {{=}} (''a'' + ''bi''){{sup|''n''}}}} の偏角は、
:<math>\arg z^n = n\arctan \frac{b}{a}</math>
である。この式を利用すると、マチンの公式の左辺は
:<math>(5+i)^4 \, (239+i)^{-1} =\frac{(5+i)^4}{239+i} =2+2i</math>
の左辺の式の偏角に等しいと分かる。この式の右辺の偏角は {{math|{{pi}} / 4}} であるためマチンの公式が示される。

== マチンの公式による計算 ==
マチンの公式を
:<math>\pi =16\arctan \frac{1}{5} -4\arctan \frac{1}{239}</math>
の形にし、arctan ''x'' をグレゴリー級数に直して、それぞれ最初の方の項だけを計算して、部分和
:<math>d(m)=16\sum_{n=0}^{3m+2} \frac{(-1)^n}{2n+1} \left( \frac{1}{5} \right)^{2n+1} -4\sum_{n=0}^m \frac{(-1)^n}{2n+1} \left( \frac{1}{239} \right)^{2n+1}</math>
を取る。
:和を取る項数がそれぞれ 3''m'' + 3 項と ''m'' + 1 項であり異なっている。これは、{{sfrac|1|5}} と {{sfrac|1|239}} の値が大きく異なるので、計算する項の値の大きさを近付けるために項数を補正しているのである。''m'' が 1 増えるたびに、計算すべき項数は 4 増える。

''m'' = 1 から ''m'' = 10 まで計算すると次表のようになる。桁数の欄は実際の円周率の値と一致している小数点以下の桁数である。参考までに末尾に {{π}} の値も載せた。
:未知の円周率を計算するときには、誤差を評価し、有効な桁数を調べなければならないが、ここでは既に知られている円周率の値と比べて、一致することを確認するだけにとどめる。
{|class="wikitable"
!''m''!!''d''(''m'')!!style="white-space:nowrap"|桁数!!style="white-space:nowrap"|項数
|-
|align="right"|0
|'''3.141'''62 …
|align="right"|3
|align="right"|4
|-
|align="right"|1
|'''3.14159 265'''26 …
|align="right"|8
|align="right"|8
|-
|align="right"|2
|'''3.14159 26535 89'''83 …
|align="right"|12
|align="right"|12
|-
|align="right"|3
|'''3.14159 26535 89793 23'''63 …
|align="right"|17
|align="right"|16
|-
|align="right"|4
|'''3.14159 26535 89793 23846 2'''75 …
|align="right"|21
|align="right"|20
|-
|align="right"|5
|'''3.14159 26535 89793 23846 26433'''77 …
|align="right"|25
|align="right"|24
|-
|align="right"|6
|'''3.14159 26535 89793 23846 26433 83279'''81 …
|align="right"|30
|align="right"|28
|-
|align="right"|7
|'''3.14159 26535 89793 23846 26433 83279 5028'''66 …
|align="right"|34
|align="right"|32
|-
|align="right"|8
|'''3.14159 26535 89793 23846 26433 83279 50288 419'''81 …
|align="right"|38
|align="right"|36
|-
|align="right"|9
|'''3.14159 26535 89793 23846 26433 83279 50288 41971 693'''41 …
|align="right"|43
|align="right"|40
|-
|align="right"|10
|'''3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37'''84 …
|align="right"|47
|align="right"|44
|-
|align="right"|'''{{π}}'''||3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37510 58209 7494 …
|colspan="2" align="center"|（参考）
|}
''d''(''m'') は 4(''m'' + 1) 個の項の[[加法|足し算]]または[[減法|引き算]]によって計算されるので、''m'' = 10 のときは 44 項の和や差を計算していることになる。ここで普通のグレゴリー級数を用いた場合の値を見てみると

''x'' = 1 のときのグレゴリー級数
:<math>\pi =4\arctan 1=4\sum_{n=0}^{\infin} \frac{(-1)^n}{2n+1}</math>
は、非常に収束が遅く、''n'' = 50 までで打ち切って計算してみると
:<math>4\sum_{n=0}^{50} \frac{(-1)^n}{2n+1} =3.1611\cdots</math>
となり小数点以下 1 桁までしか円周率と一致していない。

シャープの用いた ''x'' = {{sfrac|1|{{sqrt|3}}}} の場合のグレゴリー級数
:<math>\pi =6\arctan \frac{1}{\sqrt{3}} =6\sum_{n=0}^{\infin} \frac{(-1)^n}{2n+1} \left( \frac{1}{\sqrt{3}} \right)^{2n+1}</math>
で同じように  ''n'' = 50 までで打ち切って計算すると
:<math>6\sum_{n=0}^{50} \frac{(-1)^n}{2n+1} \left( \frac{1}{\sqrt{3}} \right)^{2n+1} =3.14159 26535 89793 23846 26433 95 \cdots</math>
となり、円周率の実際の値とは小数点以下 25 桁まで一致している。上の表で見るとマチンの公式では ''d''(5) で 25 桁まで一致しており、そのときの計算に用いた項の数は 4 × (5 + 1) = 24 項であるので、シャープによる計算のほぼ半分の項数によって、小数点以下 25 桁までの円周率が得られている。

== マチンの公式の類似 ==
マチンの公式に類似した式は比較的探しやすいため、非常に多くの形の式が見つかっている。この節では、その中のほんの一部を紹介する。複素数を用いたマチンの公式の証明と同様の計算を用いるなどして、[[計算機]]を用いて公式を探索していくことも可能である。

=== 2項よりなる公式 ===
[[レオンハルト・オイラー|オイラー]]による公式（[[1748年]]）
:<math>\frac{\pi}{4} =\arctan \frac{1}{2} +\arctan \frac{1}{3}</math>
:<math>\arctan \frac{1}{p} =\arctan \frac{1}{p+q} +\arctan \frac{q}{p^2 +pq+1}</math>
:ただし ''p'', ''q'' は正の実数
:<math>\frac{\pi}{4} =5\arctan \frac{1}{7} +2\arctan \frac{3}{79}</math>
[[ヤコブ・ハーマン]] (''[[:en:Jacob Hermann (mathematician)|Jacob Hermann]]'',1678 - 1733) による式
:<math>\frac{\pi}{4} =2\arctan \frac{1}{2} -\arctan \frac{1}{7}</math>
ハットン（''Charles Hutton'',1737 - 1823）による式（1776年）
:<math>\frac{\pi}{4} =3\arctan \frac{1}{4} +\arctan \frac{5}{99}</math>
:<math>\frac{\pi}{4} =2\arctan \frac{1}{3} +\arctan \frac{1}{7}</math>
:下の式は、1779年にオイラーも独立に再発見している。

=== 3項以上よりなる公式 ===
[[カール・フリードリヒ・ガウス|ガウス]]による公式（1863年）
:<math>\frac{\pi}{4} =12\arctan \frac{1}{18} +8\arctan \frac{1}{57} -5\arctan \frac{1}{239}</math>
ストーマー ([[:en:Carl Størmer|''Fredrik Carl Mulertz Stormer'']], 1874-1957) による公式 （[[1896年]]）
:<math>\frac{\pi}{4} =44\arctan \frac{1}{57} +7\arctan \frac{1}{239} -12\arctan \frac{1}{682} +24\arctan \frac{1}{12943}</math>
:<math>\frac{\pi}{4} =6\arctan \frac{1}{8} +2\arctan \frac{1}{57} +\arctan \frac{1}{239}</math>
[[高野喜久雄]]による公式（[[1982年]]）
:<math>\frac{\pi}{4} =12\arctan \frac{1}{49} +32\arctan \frac{1}{57} -5\arctan \frac{1}{239} +12\arctan \frac{1}{110443}</math>
シムソン（Robert Simson;1687生, 1768没; 初等幾何学の[[シムソン線]]の発見者）による公式<ref>マチン・シムソンの論文 I.Tweddle(1991), John Machin and Robert Simson on Inverse-tangent series for π, Arch. Hist. Exact Sci. の'''42''',p.1~14による。</ref>
:<math>\frac{\pi}{4} =8\arctan \frac{1}{10} -4\arctan \frac{1}{515} -\arctan \frac{1}{239}</math>

== 注釈 ==
<references/>

== 関連項目 ==
* [[円周率の歴史]]

== 外部リンク ==
* {{MathWorld|title=Machin-Like Formulas|urlname=Machin-LikeFormulas}}
* [https://drive.google.com/file/d/1DwdVhIUhUlIIZFgGwx03CGVrbmmqPxvJ/view?usp=drive_link 逆正接関数による円周率計算]

{{デフォルトソート:まちんのこうしき}}
[[Category:アルゴリズム]]
[[Category:人名を冠した数式]]
[[Category:円周率]]
[[Category:解析学]]
[[Category:数学に関する記事]]
[[Category:数学のエポニム]]