239：Machinの公式と高野喜久雄の公式

$239$ は円周率に関するMachinの公式や高野喜久雄の公式に現れる素数です。

この記事では逆正接関数 $\arctan x$ の主値は $\displaystyle -\frac{\pi}{2} < \arctan x < \frac{\pi}{2}$ で考えることにします。

逆正接関数のTaylor展開

$\displaystyle \arctan x = \sum_{n=0}^{\infty}\frac{(-1)^n}{2n+1}x^{2n+1}, \ \ \left|x\right| < 1$

において $x=1$ を代入することによって有名なLeibnizの公式

$\displaystyle \frac{\pi}{4} = 1-\frac{1}{3}+\frac{1}{5}-\frac{1}{7}+\cdots$

が得られます。収束半径上の代入で条件収束ですが、この公式が正しいことは実際に確認することが出来ます*1。また、Leibnizは1673年にこの公式を得ていますが、実際には同時代のGregory、更に遡って14世紀にはMadhavaが発見していたようです*2。

この公式は円周率を奇数を用いた単純な級数で表す美しいものですが*3、円周率の計算に関して言えば非常に収束が遅くて使いものになりません。

一方、1706年にMachinが発見したMachinの公式

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

は非常に収束が速く、円周率の計算に実際に用いられていました。

$\displaystyle \pi = \lim_{n \to \infty} \left( \sum_{k=0}^{3n+2}\frac{(-1)^k}{2k+1}\cdot \frac{16}{5^{2k+1}} - \sum_{k=0}^n\frac{(-1)^k}{2n+1}\cdot \frac{4}{239^{2k+1}} \right)$

と書き直すと、 $n=0$ で円周率の小数点以下第3位、 $n=1$ で第7位、 $n=10$ で第43位まで一致します*4。

少し脱線しますが、高木貞治先生のエッセイ『昔と今（円周率をめぐって）』の文章がとても印象に残っているので引用します。

$\pi$ を $707$ 桁計算したのはシャンクス(Shanks)という男で、その結果は1873-75年の英国王立学士院記事（英名略）に載っているそうであるが、予はまだそれを見たことがない。それを見るひまがない。ひまがあっても興味がない。 $707$ 桁の数字そのものに、なんら科学的の価値を認め得ないのである（ここに実際は第528位に誤りがあることの注釈）。
世には計算狂－といっては、ひどすぎるならば、計算家、それは少しもったいない気がするから、まあ中を取って、計算屋とでも言うべきものがある。多くは一種変態性の持ち主である。中でも低級なのは、（西洋の話だけれども）寄席などへ出て、演技をする。５桁６桁といった大きな数を２人の客から随時に提出させて、俗歌など歌いながら、暗算で掛算をして見せるといった連中である。こういうのとは少し違うが、円周率の計算家というものが近世には多数出ている。

これは『数学の自由性』（ちくま学芸文庫）で読むことができます。

話を戻すと、Machinの公式に似た公式は大量に知られています。

例えば、Eulerは

$\displaystyle \frac{\pi}{4} = 5\arctan \frac{1}{7} - 2\arctan \frac{3}{79}$

を発見しており、Gaussは

$\displaystyle \frac{\pi}{4} = 12\arctan \frac{1}{18} +8\arctan \frac{1}{57} -5\arctan \frac{1}{239}$

を発見しています。

Störmerが1896年に発見した

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

および、高野喜久雄が1982年に発見した

$\displaystyle \frac{\pi}{4} = 12\arctan \frac{1}{49} +32\arctan \frac{1}{57} -5\arctan \frac{1}{239}+12\arctan \frac{1}{110443}$

は金田康正による円周率の記録： 1兆2411億7730万桁（2002/12/6）に用いられました。

ちなみに高野喜久雄（1927/11/20-2006/5/1）は詩人です。

Machinの公式の証明

$\displaystyle [ x ] := \arctan x$ とするとき、次が成り立ちます：

命題 $xy \neq 1$ なる $x, y \in \mathbb{R}$ に対して

$\displaystyle [ x] + [ y] = \left[ \frac{x+y}{1-xy} \right]$

が成り立つ。

これは正接の加法定理の単なる書き換えです。

$-[x] =[-x]$

および

$\displaystyle 2[ x] = \left[ \frac{2x}{1-x^2} \right], \ \ (x \neq \pm 1)$

が成り立つので、

$\displaystyle 4\left[ \frac{1}{5} \right] = 2\left[ \frac{5}{12} \right] = \left[ \frac{120}{119} \right],$

$\displaystyle 4\left[ \frac{1}{5} \right] - \left[ \frac{1}{239} \right] = \left[ \frac{120}{119} \right] + \left[ -\frac{1}{239} \right] = [ 1]$

と計算されます。よって、Marchinの公式が証明されました。

高野喜久雄の公式を証明したければ、 $[ \ ]$ に対する計算規則を用いて

$\displaystyle 12\left[ \frac{1}{39} \right]+32\left[ \frac{1}{57} \right]-5\left[ \frac{1}{239} \right]+12\left[ \frac{1}{110443} \right] = [ 1]$

を確認すればよいわけです。つまり、Machin型の公式の証明はいつだって機械的に行えることがわかります*5。

一方、新しい公式を発見するのは簡単ではありません。

*1: $x$ についての恒等式

$\displaystyle \frac{1}{1+x^2}=\sum_{k=0}^n(-1)^kx^{2k}+\frac{(-1)^{n+1}x^{2n+2}}{1+x^2}$

の両辺を区間 $[0, 1]$ で積分する。左辺は

$\displaystyle \int_0^1\frac{dx}{1+x^2} = \left[ \arctan x\right]_0^1 = \frac{\pi}{4}$

であり、右辺は

$\displaystyle \sum_{k=0}^n\frac{(-1)^k}{2k+1}+(-1)^{n+1}\int_0^1\frac{x^{2n+2}}{1+x^2}dx$

となる。

$\displaystyle 0 < \int_0^1\frac{x^{2n+2}}{1+x^2}dx < \int_0^1x^{2n+2}dx = \frac{1}{2n+3} \xrightarrow{n \to \infty} 0$

によってLeibnizの公式が得られる。

*2:Stiglerの法則。 integers.hatenablog.com

*3:もちろん、数学的にも重要です。

*4:この和の取り方の若干の工夫はWikipediaの記事を参照しました。別に工夫しなくてもよいですが、数値を参照するためにWikipediaに合わせました。

*5:タンジェントの加法定理とその拡張 | 高校数学の美しい物語で紹介されているように

$\displaystyle \sum_{k=1}^n[ x_k ] = \sum_{0\leq i \leq n, i:\text{奇数}}e_i \left/ \sum_{0 \leq i \leq n, i:\text{偶数}}e_i \right.$

が成り立ちます。ただし、 $e_i$ は $x_k$ 達の $i$ 次の基本対称式（ $e_0=1$ ）で、右辺の分母が $0$ にならない場合を考えます。

INTEGERS

数、特に整数に関する記事。

239：Machinの公式と高野喜久雄の公式

Machinの公式の証明