History of numerical approximations of π

From Wikipedia, the free encyclopedia

This page is about the history of numerical approximations of the mathematical constant π.

There is a summarizing table at chronology of computation of π.

See also history of π for other apects of the evolution of our knowledge about mathematical properties of π.

Contents 1 Early history 2 Middle ages 3 16th to 19th centuries 4 20th century 5 Less accurate approximations 5.1 Biblical value 5.2 The Indiana bill 6 Development of efficient formulae 6.1 Machin-like formulae 6.2 Other classical formulae 6.3 Modern algorithms 6.4 Formulae for binary digits 6.5 Miscellaneous formulæ 7 References

[edit] Early history

As early as the 19th century BC, Babylonian mathematicians were using π = 25/8, which is almost within 0.5% of the exact value.

An Egyptian scribe named Ahmes wrote the oldest known text to give an approximate value for π. The Rhind Mathematical Papyrus dates from the Egyptian Second Intermediate Period—though Ahmes stated that he copied a Middle Kingdom papyrus (i.e. from before 1650 BC) —and describes the value in such a way that the result obtained comes out to 256 divided by 81 or 3.160.

The Indian astronomer Yajnavalkya gave astronomical calculations in the Shatapatha Brahmana (c. 9th century BC) that led to a fractional approximation of π = 339/108 (which equals 3.13888…, which is correct to two decimal places when rounded).

In the third century BC, Archimedes proved the sharp inequalities 3 + 10/71 < π < 3 + 1/7, by means of regular 96-gons. (A simple modern proof that indeed 3 + 1 over 7 exceeds π may be founded on the knowledge of the differential of the arcus tangens function.)

The Chinese mathematician Liu Hui in AD 263 computed π to 3.141014, which is correct to 3 decimal places, though he suggested that 3.14 was a good enough approximation.

[edit] Middle ages

Until 1000 AD, π was known to less than 10 decimal digits only.

The Indian mathematician and astronomer Aryabhata in the 5th century gave an accurate approximation for π. He wrote "Add four to one hundred, multiply by eight and then add sixty-two thousand. The result is approximately the circumference of a circle of diameter twenty thousand. By this rule the relation of the circumference to diameter is given." In other words (4+100)×8 + 62000 is the circumference of a circle with diameter 20000. This provides a value of π = 62832/20000 = 3.1416, correct when rounded off to 4 decimal places.

The Chinese mathematician and astronomer Zu Chongzhi in the 5th century computed π between 3.1415926 and 3.1415927, which was correct to 7 decimal places. He gave two other approximations of π: 355/113 and 22/7.

In the 14th century, the Indian mathematician and astronomer Madhava of Sangamagrama gave two methods for computing the value of π. One of these methods is to obtain a rapidly converging series by transforming the original infinite series of π. By doing so, he obtained the infinite series

and used the first 21 terms to compute an approximation of π correct to 11 decimal places as 3.14159265359.

The other method he used was to add a remainder term to the original series of π. He used the remainder term

in the infinite series expansion of to improve the approximation of π to 13 decimal places of accuracy when n = 75.

The Persian Muslim mathematician and astronomer Ghyath ad-din Jamshid Kashani (1380 – 1429) correctly computed 2π to 9 sexagesimal (base 60) digits.^[1] This figure is equivalent to 16 decimal (base 10) digits as

which equates to

He achieved this level of accuracy by calculating the perimeter of a regular polygon with 3 × 2¹⁸ sides.^{[citation needed]}

[edit] 16th to 19th centuries

The German mathematician Ludolph van Ceulen (circa 1600) computed the first 32 decimal places of π. He was so proud of this accomplishment that he had them inscribed on his tombstone.

The Slovene mathematician Jurij Vega in 1789 calculated the first 140 decimal places for π of which the first 126 were correct [1] and held the world record for 52 years until 1841, when William Rutherford calculated 208 decimal places of which the first 152 were correct. Vega improved John Machin's formula from 1706 and his method is still mentioned today.

The English amateur mathematician William Shanks, a man of independent means, spent over 20 years calculating π to 707 decimal places (accomplished in 1873). His routine was as follows: he would calculate new digits all morning; and then he would spend all afternoon checking his morning's work. His work was made possible by the recent invention of the logarithm and its tables by Napier and Briggs. This was the longest expansion of π until the advent of the electronic digital computer a century later.

The Gauss-Legendre algorithm is used for calculating digits of π.

[edit] 20th century

In 1910, the Indian mathematician Srinivasa Ramanujan found several rapidly converging infinite series of π, including

which computes a further 8 decimal places of π with each term in the series. His series are now the basis for the fastest algorithms currently used to calculate π.

From the mid-20th century onwards, all calculations of π were done with the help of calculators or computers.

In 1944, D. F. Ferguson, with the aid of a mechanical desk calculator, found that William Shanks had made a mistake in the 528th decimal place, and that all succeeding digits were incorrect.

In the early years of the computer, the first expansion of π to 1,000,000 decimal places was computed by Maryland mathematician Dr. Daniel Shanks and his team at the United States Naval Research Laboratory (N.R.L.) in Washington, D.C. (Dr. Shanks's son Oliver Shanks, also a mathematician, states that there is no connection to William Shanks, and the family's roots are in Central Europe).

In 1961, Daniel Shanks and his team used two different power series for calculating the digital of π. For one it was known that any error would produce a value slightly high, and for the other, it was known that any error would produce a value slightly low. And hence, as long as the two series produced the same digits, there was a very high confidence that they were correct. The first 1,000,000 digits of π were published by the N.R.L.

In 1989, the Chudnovsky brothers correctly computed π to over a billion decimal places on the supercomputer IBM 3090 using the following variation of Ramanujan's infinite series of π:

In 1999, Yasumasa Kanada and his team at the University of Tokyo correctly computed π to over 200 billion decimal places on the supercomputer HITACHI SR8000/MPP (128 nodes) using another variation of Ramanujan's infinite series of π. In October 2005 they claimed to have calculated it to 1.24 trillion places. [2]

[edit] Less accurate approximations

Some approximations which have been given for π are notable in that they were less precise than previously known values.

[edit] Biblical value

It is sometimes claimed that the Bible states that π is 3, based on a passage in 1 Kings 7:23 giving measurements for the round basin located in front of the Temple in Jerusalem as having a 10 cubit diameter and a 30 cubit circumference. Rabbi Nehemiah explained this in his Mishnat ha-Middot (the earliest known Hebrew text on geometry, ca. AD 150) by saying that the diameter was measured from the outside of the brim while the circumference was measured along the inner rim. The stated dimensions would be exact if measured this way on a brim about four inches wide.

This is disputed, however, and other explanations have been offered, including that the measurements are given in round numbers (as the Hebrews tended to round off measurements to whole numbers), or that the basin may not have been exactly circular, or that the brim was wider than the bowl itself. Many reconstructions of the basin, in fact, show a wider brim extending outward from the bowl itself by several inches, which might account for the discrepancy.

[edit] The Indiana bill

The "Indiana Pi Bill" of 1897, which never passed out of committee, has been claimed to imply a number of different values for π, although the closest it comes to explicitly asserting one is the wording "the ratio of the diameter and circumference is as five-fourths to four", which would make π = 3.2.

[edit] Development of efficient formulae

[edit] Machin-like formulae

For fast calculations, one may use formulæ such as Machin's:

together with the Taylor series expansion of the function arctan(x). This formula is most easily verified using polar coordinates of complex numbers, starting with

Formulæ of this kind are known as Machin-like formulae.

[edit] Other classical formulae

Other formulæ that have been used to compute estimates of π include:

Madhava.

Euler.

Newton.

Ramanujan.

This converges extraordinarily rapidly. Ramanujan's work is the basis for the fastest algorithms used, as of the turn of the millennium, to calculate π:

David Chudnovsky and Gregory Chudnovsky.

Many other expressions for π were developed and published by the incredibly-intuitive Indian mathematician Srinivasa Ramanujan. He worked with mathematician Godfrey Harold Hardy in England for a number of years.

[edit] Modern algorithms

Extremely long decimal expansions of π are typically computed with iterative formulae like the Gauss-Legendre algorithm and Borwein's algorithm. The Salamin-Brent algorithm which was invented in 1976 is an example of the former.

Borwein's algorithm, found in 1985 by Jonathan and Peter Borwein, converges extremely fast: For and

where f(y) = (1 − y⁴)^{1 / 4}, the sequence 1 / a_k converges quartically to π, giving about 100 digits in three steps and over a trillion digits after 20 steps.

The first one million digits of π and 1/π are available from Project Gutenberg (see external links below). The current record (December 2002) by Yasumasa Kanada of Tokyo University stands at 1,241,100,000,000 digits, which were computed in September 2002 on a 64-node Hitachi supercomputer with 1 terabyte of main memory, which carries out 2 trillion operations per second, nearly twice as many as the computer used for the previous record (206 billion digits). The following Machin-like formulæ were used for this:

K. Takano (1982).

F. C. W. Störmer (1896).

These approximations have so many digits that they are no longer of any practical use, except for testing new supercomputers. (Normality of π will always depend on the infinite string of digits on the end, not on any finite computation.)

[edit] Formulae for binary digits

Main article: Bailey-Borwein-Plouffe formula

In 1997, David H. Bailey, Peter Borwein and Simon Plouffe published a paper (Bailey, 1997) on a new formula for π as an infinite series:

This formula permits one to easily compute the k^th binary or hexadecimal digit of π, without having to compute the preceding k − 1 digits. Bailey's website contains the derivation as well as implementations in various programming languages. The PiHex project computed 64-bits around the quadrillionth bit of π (which turns out to be 0).

[edit] Miscellaneous formulæ

Historically, for a long time the base 60 was used for calculations. In this base, π can be approximated to eight (decimal!) significant figures as

(The next sexagesimal digits would be '0 47', which motivates breaking the expression at this point.)

In addition, the following expressions can be used to estimate π:

• accurate to 9 digits:

• accurate to 9 places:

• This is from Ramanujan, who claimed the goddess Namagiri appeared to him in a dream and told him the true value of π.^{[citation needed]}

• accurate to 4 digits:

• accurate to 3 digits:

• Karl Popper conjectured that Plato knew this expression, that he believed it to be exactly π, and that this is responsible for some of Plato's confidence in the omnicompetence of mathematical geometry — and Plato's repeated discussion of special right triangles that are either isosceles or halves of equilateral triangles.

• The continued fraction representation of π can be used to generate successive best rational approximations. These approximations are the best possible rational approximations of π relative to the size of their denominators. Here is a list of the first of these, punctuated at significant steps:

  

[edit] References

• Bailey, David H., Borwein, Peter B., and Plouffe, Simon (April 1997). "On the Rapid Computation of Various Polylogarithmic Constants". Mathematics of Computation 66 (218): 903–913.
• Joseph, George G. (2000). The Crest of the Peacock: Non-European Roots of Mathematics.