マチンの公式

マチンの公式（英: Machin's formula）とは、1706年にイギリスの天文学者ジョン・マチンによって発見された逆正接関数 arctan x を用いた円周率を計算するための公式、すなわち π 4 = 4 arctan ⁡ 1 5 − arctan ⁡ 1 239 {\displaystyle {\dfrac {\pi }{4}}=4\arctan {\dfrac {1}{5}}-\arctan {\dfrac {1}{239}}} なる公式である。 グレゴリー級数すなわち逆正接関数 arctan x のマクローリン展開： arctan ⁡ x = ∑ n = 0 ∞ ( − 1 ) n 2 n + 1 x 2 n + 1 = x − 1 3 x 3 + 1 5 x 5 − 1 7 x 7 + 1 9 x 9 − ⋯ w h e r e | x | ≤ 1 {\displaystyle \arctan x=\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{2n+1}}\,x^{2n+1}=x-{\frac {1}{3}}\,x^{3}+{\frac {1}{5}}\,x^{5}-{\frac {1}{7}}\,x^{7}+{\frac {1}{9}}\,x^{9}-\cdots \qquad {\rm {where}}\quad \vert x\vert \leq 1} に x = 1 を代入して得られる級数： π 4 = ∑ n = 0 ∞ ( − 1 ) n 2 n + 1 = 1 − 1 3 + 1 5 − 1 7 + 1 9 − ⋯ {\displaystyle {\frac {\pi }{4}}=\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{2n+1}}=1-{\frac {1}{3}}+{\frac {1}{5}}-{\frac {1}{7}}+{\frac {1}{9}}-\cdots } はライプニッツの公式と呼ばれ、見た目は綺麗な公式であるものの、収束が非常に遅いことで知られる。しかしながら、x を十分小さく取れば見た目の綺麗さは多少損なわれるが、それなりに速く収束する級数を得ることができる。実際、エイブラハム・シャープは x = 1/√3 を用い、円周率を小数点以下71桁まで計算した。 ジョン・マチンは、さらに収束性をよくするために逆正接関数 arctan x の関係式を考え、これをグレゴリー級数と結び付けることにより、非常に収束速度が速い級数を得た。さらに、この公式を発見したマチン自身も円周率を100桁まで求めることに成功した。マチンの公式や、似たような arctan x を用いた公式は、1970年代に算術幾何平均などが用いられるようになるまでは円周率の計算に用いられ計算競争に貢献した。その後しばらくは新しいアルゴリズムによる円周率の計算が続いたが、2002年に金田康正によって高野喜久雄の公式が用いられ円周率を1兆2411億桁まで計算するという記録に結び付いた。 等式： π 4 = 4 arctan ⁡ 1 5 − arctan ⁡ 1 239 {\displaystyle {\dfrac {\pi }{4}}=4\arctan {\frac {1}{5}}-\arctan {\frac {1}{239}}} をマチンの公式という。 この項目では arctan x は主値 − π 2 < arctan ⁡ x < π 2 {\displaystyle -{\frac {\pi }{2}}<\arctan x<{\frac {\pi }{2}}} を取るものとする。 同じことであるが、逆余接関数 arccot x を用いて、 π 4 = 4 arccot ⁡ 5 − arccot ⁡ 239 {\displaystyle {\dfrac {\pi }{4}}=4\operatorname {arccot} 5-\operatorname {arccot} 239} と書かれることもある。 マチンの公式は三角関数の公式を用いて証明できる。 二倍角公式を2回用いて、 tan ⁡ ( 2 arctan ⁡ 1 5 ) = 2 tan ⁡ ( arctan ⁡ 1 5 ) 1 − tan 2 ⁡ ( arctan ⁡ 1 5 ) = 2 ⋅ 1 5 1 − ( 1 5 ) 2 = 5 12 {\displaystyle \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}}} tan ⁡ ( 4 arctan ⁡ 1 5 ) = tan ⁡ ( 2 ⋅ 2 arctan ⁡ 1 5 ) = 2 tan ⁡ ( 2 arctan ⁡ 1 5 ) 1 − tan 2 ⁡ ( 2 arctan ⁡ 1 5 ) = 2 ⋅ 5 12 1 − ( 5 12 ) 2 = 120 119 {\displaystyle \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}}} 加法定理により、 tan ⁡ ( 4 arctan ⁡ 1 5 − π 4 ) = tan ⁡ ( 4 arctan ⁡ 1 5 ) − tan ⁡ π 4 1 + tan ⁡ ( 4 arctan ⁡ 1 5 ) tan ⁡ π 4 = 120 119 − 1 1 + 1 ⋅ 120 119 = 1 239 {\displaystyle \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}}} 逆関数をとって、 arctan ⁡ 1 239 = 4 arctan ⁡ 1 5 − π 4 {\displaystyle \arctan {\frac {1}{239}}=4\arctan {\frac {1}{5}}-{\frac {\pi }{4}}} したがって、 4 arctan ⁡ 1 5 − arctan ⁡ 1 239 = π 4 {\displaystyle 4\arctan {\frac {1}{5}}-\arctan {\frac {1}{239}}={\frac {\pi }{4}}} a 、b を実数とし、a > 0 、i を虚数単位とする。 複素数 z = a + bi の偏角は、 arg ⁡ z = arctan ⁡ b a {\displaystyle \arg z=\arctan {\frac {b}{a}}} である。 複素数の偏角の範囲は arctan の主値と同じ範囲に取るものとする。 n を整数とする。ド・モアブルの定理によると zn = (a + bi)n の偏角は、 arg ⁡ z n = n arctan ⁡ b a {\displaystyle \arg z^{n}=n\arctan {\frac {b}{a}}} である。この式を利用すると、マチンの公式の左辺は ( 5 + i ) 4 ( 239 + i ) − 1 = ( 5 + i ) 4 239 + i = 2 + 2 i {\displaystyle (5+i)^{4}\,(239+i)^{-1}={\frac {(5+i)^{4}}{239+i}}=2+2i} の左辺の式の偏角に等しいと分かる。この式の右辺の偏角は π / 4 であるためマチンの公式が示される。 マチンの公式を π = 16 arctan ⁡ 1 5 − 4 arctan ⁡ 1 239 {\displaystyle \pi =16\arctan {\frac {1}{5}}-4\arctan {\frac {1}{239}}} の形にし、arctan x をグレゴリー級数に直して、それぞれ最初の方の項だけを計算して差 d ( m ) = 16 ∑ n = 0 3 m + 2 ( − 1 ) n 2 n + 1 ( 1 5 ) 2 n + 1 − 4 ∑ n = 0 m ( − 1 ) n 2 n + 1 ( 1 239 ) 2 n + 1 {\displaystyle 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}} を取る。 和を取る項数がそれぞれ 3m + 3 項と m + 1 項であり異なっている。これは、1/5 と 1/239 の値が大きく異なるので、計算する項の値の大きさを近付けるために項数を補正しているのである。m が 1 増えるたびに、計算すべき項数は 4 増える。 m = 1 から m = 10 まで計算すると次表のようになる。桁数の欄は実際の円周率の値と一致している小数点以下の桁数である。参考までに末尾に π の値も載せた。 未知の円周率を計算するときには、誤差を評価し、有効な桁数を調べなければならないが、ここでは既に知られている円周率の値と比べて、一致することを確認するだけにとどめる。 d(m) は 4(m + 1) 個の項の足し算または引き算によって計算されるので、m = 10 のときは 44 項の和や差を計算していることになる。ここで普通のグレゴリー級数を用いた場合の値を見てみると x = 1 のときのグレゴリー級数 π = 4 arctan ⁡ 1 = 4 ∑ n = 0 ∞ ( − 1 ) n 2 n + 1 {\displaystyle \pi =4\arctan 1=4\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{2n+1}}} は、非常に収束が遅く、n = 50 までで打ち切って計算してみると 4 ∑ n = 0 50 ( − 1 ) n 2 n + 1 = 3.1611 ⋯ {\displaystyle 4\sum _{n=0}^{50}{\frac {(-1)^{n}}{2n+1}}=3.1611\cdots } となり小数点以下 1 桁までしか円周率と一致していない。 シャープの用いた x = 1/√3 の場合のグレゴリー級数 π = 6 arctan ⁡ 1 3 = 6 ∑ n = 0 ∞ ( − 1 ) n 2 n + 1 ( 1 3 ) 2 n + 1 {\displaystyle \pi =6\arctan {\frac {1}{\sqrt {3}}}=6\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{2n+1}}\left({\frac {1}{\sqrt {3}}}\right)^{2n+1}} で同じように n = 50 までで打ち切って計算すると 6 ∑ n = 0 50 ( − 1 ) n 2 n + 1 ( 1 3 ) 2 n + 1 = 3.141592653589793238462643395 ⋯ {\displaystyle 6\sum _{n=0}^{50}{\frac {(-1)^{n}}{2n+1}}\left({\frac {1}{\sqrt {3}}}\right)^{2n+1}=3.141592653589793238462643395\cdots } となり、円周率の実際の値とは小数点以下 25 桁まで一致している。上の表で見るとマチンの公式では d(5) で 25 桁まで一致しており、そのときの計算に用いた項の数は 4 × (5 + 1) = 24 項であるので、シャープによる計算のほぼ半分の項数によって、小数点以下 25 桁までの円周率が得られている。 マチンの公式に類似した式は比較的探しやすいため、非常に多くの形の式が見つかっている。この節では、その中のほんの一部を紹介する。複素数を用いたマチンの公式の証明と同様の計算を用いるなどして、計算機を用いて公式を探索していくことも可能である。 オイラーによる公式（1748年） π 4 = arctan ⁡ 1 2 + arctan ⁡ 1 3 {\displaystyle {\frac {\pi }{4}}=\arctan {\frac {1}{2}}+\arctan {\frac {1}{3}}} arctan ⁡ 1 p = arctan ⁡ 1 p + q + arctan ⁡ q p 2 + p q + 1 {\displaystyle \arctan {\frac {1}{p}}=\arctan {\frac {1}{p+q}}+\arctan {\frac {q}{p^{2}+pq+1}}} ただし p, q は正の実数 π 4 = 5 arctan ⁡ 1 7 + 2 arctan ⁡ 3 79 {\displaystyle {\frac {\pi }{4}}=5\arctan {\frac {1}{7}}+2\arctan {\frac {3}{79}}} ヤコブ・ハーマン (Jacob Hermann,1678 - 1733) による式 π 4 = 2 arctan ⁡ 1 2 − arctan ⁡ 1 7 {\displaystyle {\frac {\pi }{4}}=2\arctan {\frac {1}{2}}-\arctan {\frac {1}{7}}} ハットン（Charles Hutton,1737 - 1823）による式（1776年） π 4 = 3 arctan ⁡ 1 4 + arctan ⁡ 5 99 {\displaystyle {\frac {\pi }{4}}=3\arctan {\frac {1}{4}}+\arctan {\frac {5}{99}}} π 4 = 2 arctan ⁡ 1 3 + arctan ⁡ 1 7 {\displaystyle {\frac {\pi }{4}}=2\arctan {\frac {1}{3}}+\arctan {\frac {1}{7}}} 下の式は、1779年にオイラーも独立に再発見している。 ガウスによる公式（1863年） π 4 = 12 arctan ⁡ 1 18 + 8 arctan ⁡ 1 57 − 5 arctan ⁡ 1 239 {\displaystyle {\frac {\pi }{4}}=12\arctan {\frac {1}{18}}+8\arctan {\frac {1}{57}}-5\arctan {\frac {1}{239}}} ストーマー (Fredrik Carl Mulertz Stormer, 1874-1957) による公式 （1896年） π 4 = 44 arctan ⁡ 1 57 + 7 arctan ⁡ 1 239 − 12 arctan ⁡ 1 682 + 24 arctan ⁡ 1 12943 {\displaystyle {\frac {\pi }{4}}=44\arctan {\frac {1}{57}}+7\arctan {\frac {1}{239}}-12\arctan {\frac {1}{682}}+24\arctan {\frac {1}{12943}}} π 4 = 6 arctan ⁡ 1 8 + 2 arctan ⁡ 1 57 + arctan ⁡ 1 239 {\displaystyle {\frac {\pi }{4}}=6\arctan {\frac {1}{8}}+2\arctan {\frac {1}{57}}+\arctan {\frac {1}{239}}} 高野喜久雄による公式（1982年） π 4 = 12 arctan ⁡ 1 49 + 32 arctan ⁡ 1 57 − 5 arctan ⁡ 1 239 + 12 arctan ⁡ 1 110443 {\displaystyle {\frac {\pi }{4}}=12\arctan {\frac {1}{49}}+32\arctan {\frac {1}{57}}-5\arctan {\frac {1}{239}}+12\arctan {\frac {1}{110443}}} シムソン（Robert Simson;1687生, 1768没; 初等幾何学のシムソン線の発見者）による公式 π 4 = 8 arctan ⁡ 1 10 − 4 arctan ⁡ 1 515 − arctan ⁡ 1 239 {\displaystyle {\frac {\pi }{4}}=8\arctan {\frac {1}{10}}-4\arctan {\frac {1}{515}}-\arctan {\frac {1}{239}}} ^ マチン・シムソンの論文 I.Tweddle(1991), John Machin and Robert Simson on Inverse-tangent series for π, Arch. Hist. Exact Sci. の42,p.1~14による。 円周率の歴史 Weisstein, Eric W. "Machin-Like Formulas". mathworld.wolfram.com (英語).