マチンの公式の概念図。逆正接関数 arctan x は偏角として考えることができるため、マチンの公式は上図のように解釈することができる。
マチンの公式 （英 : 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}}}
複素平面上における、逆正接関数による偏角の表現。
a b 実数 とし、a > 0 i 虚数単位 とする。
複素数 z = a + bi
arg
z
=
arctan
b
a
{\displaystyle \arg z=\arctan {\frac {b}{a}}}
である。
複素数の偏角の範囲は arctan の主値と同じ範囲に取るものとする。
n 整数 とする。ド・モアブルの定理 によると z n 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 まで計算すると次表のようになる。桁数の欄は実際の円周率の値と一致している小数点以下の桁数である。参考までに末尾に π の値も載せた。
未知の円周率を計算するときには、誤差を評価し、有効な桁数を調べなければならないが、ここでは既に知られている円周率の値と比べて、一致することを確認するだけにとどめる。
m d (m )桁数
項数
0
3.141 62 … 3
4
1
3.14159 265 26 … 8
8
2
3.14159 26535 89 83 … 12
12
3
3.14159 26535 89793 23 63 … 17
16
4
3.14159 26535 89793 23846 2 75 … 21
20
5
3.14159 26535 89793 23846 26433 77 … 25
24
6
3.14159 26535 89793 23846 26433 83279 81 … 30
28
7
3.14159 26535 89793 23846 26433 83279 5028 66 … 34
32
8
3.14159 26535 89793 23846 26433 83279 50288 419 81 … 38
36
9
3.14159 26535 89793 23846 26433 83279 50288 41971 693 41 … 43
40
10
3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37 84 … 47
44
π 3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37510 58209 7494 …
（参考）
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
π
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 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没; 初等幾何学のシムソン線 の発見者）による公式[ 1]
π
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による。
