Machin-like formula

From Wikipedia, the free encyclopedia

In mathematics, Machin-like formulas are a class of identities involving π = 3.14159... that generalize John Machin's formula from 1706:

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

which he used along with the Taylor series expansion of arctan to compute π to 100 decimal places.

Machin-like formulas have the form

\frac{\pi}{4} = \sum_{n}^N a_n \arctan\frac{1}{b_n}

with an and bn integers.

The same method is still among the most efficient known for computing a large number of digits of π with digital computers.

Contents

[edit] Derivation

To understand where this formula comes from, start with following basic ideas:

\frac{\pi}{4} = \arctan(1)
\tan(2\arctan(a)) = \frac{2 a} { 1 - a ^ 2} (tangent double angle identity)
\tan(\arctan(a)-\arctan(b)) = \frac{a-b} { 1 + a b} (tangent difference identity)
\frac{\pi}{16} = 0.196349\dots (approximately)
\arctan\left(\frac{1}{5}\right) = \arctan(0.2) = 0.197395\dots (approximately)

In other words, for small numbers, arctangent is to a good approximation just the identity function. This leads to the possibility that a number q can be found such that

\frac{\pi}{16} = \arctan(\frac{1}{5}) - \frac{1}{4} \arctan(q).

Using elementary algebra and the identities given above, we can solve for q:

q = \tan\left(4 \arctan\left(\frac{1}{5}\right) - \frac{\pi}{4}\right) = \frac{\tan\left(4 \arctan\left(\frac{1}{5}\right)\right) - 1} { 1 + \tan\left(4 \arctan\left(\frac{1}{5}\right)\right)}

Similarly, two applications of the double angle identity yields

\tan\left(4 \arctan\left(\frac{1}{5}\right)\right) = \frac{120}{119}

and so

q = \frac{\frac{120}{119} - 1}{1 +\frac{120}{119}} = \frac{1}{239}.

[edit] Two-term formulas

Machin-like formulas with two terms include Euler's

\frac{\pi}{4} = \arctan\frac{1}{2} + \arctan\frac{1}{3},

Hermann's,

\frac{\pi}{4} = 2 \arctan\frac{1}{2} - \arctan\frac{1}{7},

and Hutton's

\frac{\pi}{4} = 2 \arctan\frac{1}{3} + \arctan\frac{1}{7}.

Many such formulas can be calculated using arctan identities. For example, it holds that for any two positive integers, such that c < d:

\frac{\pi}{4} = \arctan\frac{c}{d} + \arctan\frac{d - c}{c + d}.

[edit] More terms

The current record for digits of π, 1,241,100,000,000, by Yasumasa Kanada of Tokyo University, was obtained in 2002. A 64-node Hitachi supercomputer with 1 terabyte of main memory, performing 2 trillion operations per second, was used to evaluate the following Machin-like formulas:

\frac{\pi}{4} = 12 \arctan\frac{1}{49} + 32 \arctan\frac{1}{57} - 5 \arctan\frac{1}{239} + 12 \arctan\frac{1}{110443}
Kikuo Takano (1982).
\frac{\pi}{4} = 44 \arctan\frac{1}{57} + 7 \arctan\frac{1}{239} - 12 \arctan\frac{1}{682} + 24 \arctan\frac{1}{12943}
F. C. W. Störmer (1896).

The more efficient currently known Machin-like formulas for computing:

\frac{\pi}{4} = 183 \arctan\frac{1}{239} + 32 \arctan\frac{1}{1023} - 68 \arctan\frac{1}{5832} + 12 \arctan\frac{1}{110443} -12\arctan\frac{1}{4841182} -100 \arctan\frac{1}{6826318}
黃見利(Hwang Chien-Lih) (1997).
\frac{\pi}{4} = 183 \arctan\frac{1}{239} + 32 \arctan\frac{1}{1023} - 68 \arctan\frac{1}{5832} + 12 \arctan\frac{1}{113021} -100 \arctan\frac{1}{6826318} -12 \arctan\frac{1}{33366019650} +12\arctan\frac{1}{43599522992503626068}
黃見利(Hwang Chien-Lih) (2003).

[edit] External links

In other languages