Tuesday, January 25, 2011

The history of Pi

The name of the Greek letter π is pi.[8] The name pi is commonly used as an alternative to using the Greek letter. As a mathematical symbol, the Greek letter is not capitalized (Π) even at the beginning of a sentence, and instead the lower case (π) is used at the beginning of a sentence. When referring to this constant, the symbol π is always pronounced "pie" in English, which is the conventional English pronunciation of the Greek letter. The constant is named "π" because "π" is the first letter of the Greek word περίμετρος (perimeter), probably referring to its use in the formula perimeter/diameter which is constant for all circles, the word "perimeter" being synonymous here with "circumference."[9] William Jones was the first to use the Greek letter in this way, in 1706,[10] and it was later popularized by Leonhard Euler in 1737.[11][12] William Jones wrote:
There are various other ways of finding the Lengths or Areas of particular Curve Lines, or Planes, which may very much facilitate the Practice; as for instance, in the Circle, the Diameter is to the Circumference as 1 to  ...  3.14159, &c. = π[13]
The capital letter pi (Π) has a completely different mathematical meaning; it is used for expressing products (notice that the word "product" begins with the letter "p" just like "perimeter/diameter" does). It can also refer to the osmotic pressure of a solution.

Geometric definition

In Euclidean plane geometry, π is defined as the ratio of a circle's circumference C to its diameter d:[9]
 \pi = \frac{C}{d}.
The ratio C/d is constant, regardless of a circle's size. For example, if a circle has twice the diameter d of another circle it will also have twice the circumference C, preserving the ratio C/d.
Alternatively π can be defined as the ratio of a circle's area A to the area of a square whose side is equal to the radius r of the circle:[9][14]
 \pi = \frac{A}{r^2}.
These definitions depend on results of Euclidean geometry, such as the fact that all circles are similar, and the fact that the right-hand-sides of these two equations are equal to each other (i.e. the area of a disk is Cr/2). These two geometric definitions can be considered a problem when π occurs in areas of mathematics that otherwise do not involve geometry. For this reason, mathematicians often prefer to define π without reference to geometry, instead selecting one of its analytic properties as a definition. A common choice is to define π as twice the smallest positive x for which the trigonometric function cos(x) equals zero.[15]
Circumference = π × diameter
Area of the circle equals π times the area of the shaded square
Because π is a transcendental number, squaring the circle is not possible in a finite number of steps using the classical tools of compass and straightedge.

Irrationality and transcendence

π is an irrational number, meaning that it cannot be written as the ratio of two integers. π is also a transcendental number, meaning that there is no polynomial with rational coefficients for which π is a root.[16] An important consequence of the transcendence of π is the fact that it is not constructible. Because the coordinates of all points that can be constructed with compass and straightedge are constructible numbers, it is impossible to square the circle: that is, it is impossible to construct, using compass and straightedge alone, a square whose area is equal to the area of a given circle.[17] This is historically significant, for squaring a circle is one of the easily understood elementary geometry problems left to us from antiquity; many amateurs in modern times have attempted to solve each of these problems, and their efforts are sometimes ingenious, but in this case, doomed to failure: a fact not always understood by the amateur involved.[18][19]

Decimal representation

The decimal representation of π truncated to 50 decimal places is:[20]
π = 3.14159265358979323846264338327950288419716939937510...
Various online web sites provide π to many more digits.[21] While the decimal representation of π has been computed to more than a trillion (1012) digits,[22] elementary applications, such as estimating the circumference of a circle, will rarely require more than a dozen decimal places. For example, the decimal representation of π truncated to 11 decimal places is good enough to estimate the circumference of any circle that fits inside the Earth with an error of less than one millimetre, and the decimal representation of π truncated to 39 decimal places is sufficient to estimate the circumference of any circle that fits in the observable universe with precision comparable to the radius of a hydrogen atom.[23][24]
Because π is an irrational number, its decimal representation does not repeat, and therefore does not terminate. This sequence of non-repeating digits has fascinated mathematicians and laymen alike, and much effort over the last few centuries has been put into computing ever more of these digits and investigating π's properties.[25] Despite much analytical work, and supercomputer calculations that have determined over 1 trillion digits of the decimal representation of π, no simple base-10 pattern in the digits has ever been found.[26] Digits of the decimal representation of π are available on many web pages, and there is software for calculating the decimal representation of π to billions of digits on any personal computer.

Estimating π

Number system Approximation of π
Binary 11.00100100001111110110...[27]
Octal 3.11037552421026430215...
Decimal 3.14159265358979323846264338327950288...
Hexadecimal 3.243F6A8885A308D31319...[28]
Rational approximations 3, 227, 333106, 355113, 103993/33102, ...[29] (listed in order of increasing accuracy)
Continued fraction [3;7,15,1,292,1,1,1,2,1,3,1,14,2,1,1...][30] (This fraction is not periodic. Shown in linear notation)
3.1415926535 8979323846 2643383279 5028841971 6939937510 5820974944 5923078164 0628620899 8628034825 3421170679 8214808651 3282306647 0938446095 5058223172 5359408128 4811174502 8410270193 8521105559 6446229489 5493038196 4428810975 6659334461 2847564823 3786783165 2712019091 4564856692 3460348610 4543266482 1339360726 0249141273 7245870066 0631558817 4881520920 9628292540 9171536436 7892590360 0113305305 4882046652 1384146951 9415116094 3305727036 5759591953 0921861173 8193261179 3105118548 0744623799 6274956735 1885752724 8912279381 8301194912 9833673362 4406566430 8602139494 6395224737 1907021798 6094370277 0539217176 2931767523 8467481846 7669405132 0005681271 4526356082 7785771342 7577896091 7363717872 1468440901 2249534301 4654958537 1050792279 6892589235 4201995611 2129021960 8640344181 5981362977 4771309960 5187072113 4999999837 2978049951 0597317328 1609631859 5024459455 3469083026 4252230825 3344685035 2619311881 7101000313 7838752886 5875332083 8142061717 7669147303 5982534904 2875546873 1159562863 8823537875 9375195778 1857780532 1712268066 1300192787 6611195909 2164201989 3809525720 1065485863 2788659361 5338182796 8230301952 0353018529 6899577362 2599413891 2497217752 8347913151 5574857242 4541506959
An estimate of π accurate to 1120 decimal digits was obtained using a gear-driven calculator in 1948, by John Wrench and Levi Smith. This was the most accurate estimate of π before electronic computers came into use.[31]
The earliest numerical approximation of π is almost certainly the value 3.[32] In cases where little precision is required, it may be an acceptable substitute. That 3 is an underestimate follows from the fact that it is the ratio of the perimeter of an inscribed regular hexagon to the diameter of the circle.
π can be empirically estimated by drawing a large circle, then measuring its diameter and circumference and dividing the circumference by the diameter. Another geometry-based approach, attributed to Archimedes,[33] is to calculate the perimeter, Pn , of a regular polygon with n sides circumscribed around a circle with diameter d. Then compute the limit of a sequence as n increases to infinity:
\pi = \lim_{n \to \infty}\frac{P_{n}}{d}.
This limit converges because the more sides the polygon has, the closer the ratio approaches π. Archimedes determined the accuracy of this approach by comparing the perimeter of the circumscribed polygon with the perimeter of a regular polygon with the same number of sides inscribed inside the circle. Using a polygon with 96 sides, he computed the fractional range:[34]
3\tfrac{10}{71} < \pi < 3\tfrac{1}{7}.
π can also be calculated using purely mathematical methods. Due to the transcendental nature of π, there are no closed form expressions for the number in terms of algebraic numbers and functions.[16] Formulas for calculating π using elementary arithmetic typically include series or summation notation (such as "..."), which indicates that the formula is really a formula for an infinite sequence of approximations to π.[35] The more terms included in a calculation, the closer to π the result will get.
Most formulae used for calculating the value of π have desirable mathematical properties, but are difficult to understand without a background in trigonometry and calculus. However, some are quite simple, such as this form of the Gregory-Leibniz series:[36]
\pi = 4\sum^\infty_{k=0} \frac{(-1)^k}{2k+1} =
While that series is easy to write and calculate, it is not immediately obvious why it yields π. In addition, this series converges so slowly that nearly 300 terms are needed to calculate π correctly to 2 decimal places.[37] However, by computing this series in a somewhat more clever way by taking the midpoints of partial sums, it can be made to converge much faster. Let the sequence
\pi_{0,1} = \frac{4}{1},\ \pi_{0,2} 
=\frac{4}{1}-\frac{4}{3},\ \pi_{0,3} 
=\frac{4}{1}-\frac{4}{3}+\frac{4}{5},\ \pi_{0,4} 
=\frac{4}{1}-\frac{4}{3}+\frac{4}{5}-\frac{4}{7}, \cdots\!
and then define
\pi_{i,j} = 
\frac{\pi_{i-1,j}+\pi_{i-1,j+1}}{2}\text{ for all }i,j\ge 1
then computing π10,10 will take similar computation time to computing 150 terms of the original series in a brute-force manner, and \pi_{10,10}=3.141592653\ldots, correct to 9 decimal places. This computation is an example of the van Wijngaarden transformation.[38]
For many purposes, 3.14 or 22/7 is close enough, although engineers often use 3.1416 (5 significant figures) or 3.14159 (6 significant figures) for more precision. The approximations 22/7 and 355/113, with 3 and 7 significant figures respectively, are obtained from the simple continued fraction expansion of π. The approximation 355/113 (3.1415929...) is the best one that may be expressed with a three-digit or four-digit numerator and denominator; the next good approximation 52163/16604 (3.141592387...) requires much bigger numbers, due to the large number 292 in the continued fraction expansion of π.[29]


The Great Pyramid at Giza, constructed c.2589–2566 BC, was built with a perimeter of 1760 cubits and a height of 280 cubits; the ratio 1760/280 ≈ 2π. The same apotropaic proportions were used earlier at the Pyramid of Meidum c.2613-2589 BC and later at the pyramid of Abysir c.2453-2422. Some Egyptologists consider this to have been the result of deliberate design proportion. Verner wrote, "We can conclude that although the ancient Egyptians could not precisely define the value of π, in practice they used it".[39] Petrie, author of Pyramids and Temples of Gizeh concluded: "but these relations of areas and of circular ratio are so systematic that we should grant that they were in the builders design".[40] Others have argued that the Ancient Egyptians had no concept of pi and would not have thought to encode it in their monuments. The creation of the pyramid may instead be based on simple ratios of the sides of right angled triangles (the seked).[41] The early history of π from textual sources roughly parallels the development of mathematics as a whole.[42]


The earliest known textually evidenced approximations date from around 1900 BC; they are 256/81 (Egypt) and 25/8 (Babylonia), both within 1% of the true value.[9] The Indian text Shatapatha Brahmana gives π as 339/108 ≈ 3.139. It has been suggested that passages in the 1 Kings 7:23 and 2 Chronicles 4:2 discussing a ceremonial pool in the temple of King Solomon with a diameter of ten cubits and a circumference of thirty cubits show that the writers considered pi to have had an approximate value of three, which various authors have tried to explain away through various suggestions such as a hexagonal pool or an outward curving rim. [43]
Archimedes used the method of exhaustion to approximate the value of π.
Archimedes (287–212 BC) was the first to estimate π rigorously. He realized that its magnitude can be bounded from below and above by inscribing circles in regular polygons and calculating the outer and inner polygons' respective perimeters:[32] By using the equivalent of 96-sided polygons, he proved that 310/71 < π < 31/7.[32] The average of these values is about 3.14185.
Ptolemy, in his Almagest, gives a value of 3.1416, which he may have obtained from Apollonius of Perga.[44]
Around AD 265, the Wei Kingdom mathematician Liu Hui provided a simple and rigorous iterative algorithm to calculate π to any degree of accuracy. He himself carried through the calculation to a 3072-gon (i.e. a 3072-sided polygon) and obtained an approximate value for π of 3.1416.[45] Later, Liu Hui invented a quick method of calculating π and obtained an approximate value of 3.14 with only a 96-gon,[45] by taking advantage of the fact that the difference in area of successive polygons forms a geometric series with a factor of 4.
Around 480, the Chinese mathematician Zu Chongzhi demonstrated that π ≈ 355/113, and showed that 3.1415926 < π < 3.1415927[45] using Liu Hui's algorithm applied to a 12288-gon. This value would remain the most accurate approximation of π available for the next 900 years.
Maimonides mentions with certainty the irrationality of π in the 12th century.[46] This was proved in 1768 by Johann Heinrich Lambert.[47] In the 20th century, proofs were found that require no prerequisite knowledge beyond integral calculus. One of those, due to Ivan Niven, is widely known.[48][49] A somewhat earlier similar proof is by Mary Cartwright.[50]

Second millennium AD

Liu Hui's π algorithm
Archimedes' π algorithm
Until the second millennium AD, estimations of π were accurate to fewer than 10 decimal digits. The next major advances in the study of π came with the development of infinite series and subsequently with the discovery of calculus, which permit the estimation of π to any desired accuracy by considering sufficiently many terms of a relevant series. Around 1400, Madhava of Sangamagrama found the first known such series:
{\pi} = 4\sum^\infty_{k=0} \frac{(-1)^k}{2k+1}
 = \frac{4}{1} - \frac{4}{3} + \frac{4}{5} - \frac{4}{7} + \frac{4}{9} -
 \frac{4}{11} + \cdots\!
This is now known as the Madhava–Leibniz series[51][52] or Gregory-Leibniz series since it was rediscovered by James Gregory and Gottfried Leibniz in the 17th century. Unfortunately, the rate of convergence is too slow to calculate many digits in practice; about 4,000 terms must be summed to improve upon Archimedes' estimate. However, by transforming the series into
\pi &= \sqrt{12}\sum^\infty_{k=0} \frac{(-3)^{-k}}{2k+1} 
= \sqrt{12}\sum^\infty_{k=0} \frac{(-\frac{1}{3})^k}{2k+1} \\
&= \sqrt{12}\left(1-{1\over 3\cdot3}+{1\over5\cdot 
3^2}-{1\over7\cdot 3^3}+\cdots\right),
Madhava was able to estimate π as 3.14159265359, which is correct to 11 decimal places. The record was beaten in 1424 by the Persian mathematician, Jamshīd al-Kāshī, who gave an estimate π that is correct to 16 decimal digits.
The first major European contribution since Archimedes was made by the German mathematician Ludolph van Ceulen (1540–1610), who used a geometric method to give an estimate of π that is correct to 35 decimal digits. He was so proud of the calculation, which required the greater part of his life, that he had the digits engraved into his tombstone.[53] Pi is sometimes called "Ludolph's Constant", though not as often as it is called "Archimedes' Constant."[54]
Around the same time, the methods of calculus and determination of infinite series and products for geometrical quantities began to emerge in Europe. The first such representation was the Viète's formula,
\frac2\pi = \frac{\sqrt2}2 \cdot 
\frac{\sqrt{2+\sqrt2}}2 \cdot \frac{\sqrt{2+\sqrt{2+\sqrt2}}}2 \cdot 
found by François Viète in 1593. Another famous result is Wallis' product,
\frac{\pi}{2} = \prod^\infty_{k=1} 
\frac{(2k)^2}{(2k)^2-1} = \frac{2}{1} \cdot \frac{2}{3} \cdot 
\frac{4}{3} \cdot \frac{4}{5} \cdot \frac{6}{5} \cdot \frac{6}{7} \cdot 
\frac{8}{7} \cdot \frac{8}{9} \cdots\ = \frac{4}{3} \cdot \frac{16}{15} 
\cdot \frac{36}{35} \cdot \frac{64}{63} \cdots\!
by John Wallis in 1655. Isaac Newton derived the arcsin series for π in 1665-66 and calculated 15 digits:
\pi &= 6 \arcsin \frac {1} {2} \\
& = 6 \left( \frac {1} {2^1 \cdot 1} + \left( \frac {1} {2} \right) 
\frac {1} {2^3 \cdot 3} + \left( \frac {1 \cdot 3} {2 \cdot 4} \right) 
\frac {1} {2^5 \cdot 5} + \left( \frac{1 \cdot 3 \cdot 5} {2 \cdot 4 
\cdot 6 } \right) \frac{1} {2^7 \cdot 7} + \cdots \right) \\
&= 3 \sum_{n=0}^\infty \frac {\binom {2n} n} {16^n(2n+1)} \\
& = 3 + \frac{1}{8} + \frac{9}{640} + \frac{15}{7168} + 
\frac{35}{98304} + \frac{189}{2883584} + \frac{693}{54525952} + 
\frac{429}{167772160} + \cdots 
although he later confessed: "I am ashamed to tell you to how many figures I carried these computations, having no other business at the time."[55] It converges linearly to π, with a rate of convergence μ that adds at least three decimal places for every five terms. As n approaches infinity, μ approaches 1/4 and 1/μ approaches 4:
 \mu = \frac {(2n-1)^2} {8n(2n+1)}; \frac {1} 
{\mu} = \frac {8n(2n+1)} {(2n-1)^2} .
In 1706 John Machin was the first to compute 100 decimals of π, using the arctan series in the formula
\frac{\pi}{4} = 4 \, \arctan \frac{1}{5} - 
\arctan \frac{1}{239}\!
\arctan \, x = \sum^\infty_{k=0} \frac{(-1)^k 
x^{2k+1}}{2k+1} = x - \frac{x^3}{3} + \frac{x^5}{5} - \frac{x^7}{7} + 
Formulas of this type, now known as Machin-like formulas, were used to set several successive records and remained the best known method for calculating π well into the age of computers. A remarkable record was set by the calculating prodigy Zacharias Dase, who in 1844 employed a Machin-like formula to calculate 200 decimals of π in his head at the behest of Gauss. The best value at the end of the 19th century was due to William Shanks, who took 15 years to calculate π with 707 digits, although due to a mistake only the first 527 were correct. (To avoid such errors, modern record calculations of any kind are often performed twice, with two different formulas. If the results are the same, they are likely to be correct.)
Theoretical advances in the 18th century led to insights about π's nature that could not be achieved through numerical calculation alone. Johann Heinrich Lambert proved the irrationality of π in 1761, and Adrien-Marie Legendre also proved in 1794 π2 to be irrational. When Leonhard Euler in 1735 solved the famous Basel problem, finding the exact value of the Riemann zeta function of 2,
 \zeta(2)=\sum^\infty_{k=1} \frac{1}{k^2} = 
\frac{1}{1^2} + \frac{1}{2^2} + \frac{1}{3^2} + \frac{1}{4^2} + 
which is π2/6, he established a deep connection between π and the prime numbers. Both Legendre and Euler speculated that π might be transcendental, which was finally proved in 1882 by Ferdinand von Lindemann.

Computation in the computer age

Computer animation used to demonstrate the geometric concept of π.
Practically, a physicist needs only 39 digits of π to make a circle the size of the observable universe accurate to the size of a hydrogen atom.[23]
The advent of digital computers in the 20th century led to an increased rate of new π calculation records. John von Neumann et al. used ENIAC to compute 2037 digits of π in 1949, a calculation that took 70 hours.[56][57] Additional thousands of decimal places were obtained in the following decades, with the million-digit milestone passed in 1973. Progress was not only due to faster hardware, but also new algorithms. One of the most significant developments was the discovery of the fast Fourier transform (FFT) in the 1960s, which allows computers to perform arithmetic on extremely large numbers quickly.
In the beginning of the 20th century, the Indian mathematician Srinivasa Ramanujan found many new formulas for π, some remarkable for their elegance, mathematical depth and rapid convergence.[58] One of his formulas is the series,
\frac{1}{\pi} = \frac{2 \sqrt 2}{9801} 
\sum_{k=0}^\infty \frac{(4k)!(1103+26390k)}{(k!)^4 396^{4k}}\!
where k! is the factorial of k.
A collection of some others are in the table below:[59][60]
\pi=\frac{1}{Z}\! Z=\sum_{n=0}^{\infty } \frac{((2n)!)^3(42n+5)}
\pi=\frac{4}{Z}\! Z=\sum_{n=0}^{\infty } 
\frac{(-1)^n(4n)!(21460n+1123)} {(n!)^4{441}^{2n+1}{2}^{10n+1}}
\pi=\frac{4}{Z}\! Z=\sum_{n=0}^{\infty } \frac{(6n+1)\left ( 
\frac{1}{2} \right )^3_n} {{4^n}(n!)^3}\!
\pi=\frac{32}{Z}\! Z=\sum_{n=0}^{\infty } \left 
(\frac{\sqrt{5}-1}{2}  \right )^{8n} \frac{(42n\sqrt{5} +30n + 
5\sqrt{5}-1) \left ( \frac{1}{2} \right )^3_n} {{64^n}(n!)^3}\!
\pi=\frac{27}{4Z}\! Z=\sum_{n=0}^{\infty } \left (\frac{2}{27}  
\right )^n \frac{(15n+2)\left ( \frac{1}{2} \right )_n \left ( 
\frac{1}{3} \right )_n \left ( \frac{2}{3} \right )_n} {(n!)^3}\!
\pi=\frac{15\sqrt{3}}{2Z}\! Z=\sum_{n=0}^{\infty } \left ( \frac{4}{125} 
\right )^n \frac{(33n+4)\left ( \frac{1}{2} \right )_n \left ( 
\frac{1}{3} \right )_n \left ( \frac{2}{3} \right )_n} {(n!)^3}\!
\pi=\frac{85\sqrt{85}}{18\sqrt{3}Z}\! Z=\sum_{n=0}^{\infty } \left ( \frac{4}{85} 
\right )^n \frac{(133n+8)\left ( \frac{1}{2} \right )_n \left ( 
\frac{1}{6} \right )_n \left ( \frac{5}{6} \right )_n} {(n!)^3}\!
\pi=\frac{5\sqrt{5}}{2\sqrt{3}Z} \! Z=\sum_{n=0}^{\infty } \left ( \frac{4}{125} 
\right )^n \frac{(11n+1)\left ( \frac{1}{2} \right )_n \left ( 
\frac{1}{6} \right )_n \left ( \frac{5}{6} \right )_n} {(n!)^3}\!
\pi=\frac{2\sqrt{3}}{Z} \! Z=\sum_{n=0}^{\infty } \frac{(8n+1)\left ( 
\frac{1}{2} \right )_n \left ( \frac{1}{4} \right )_n \left ( 
\frac{3}{4} \right )_n} {(n!)^3{9}^{n}}\!
\pi=\frac{\sqrt{3}}{9Z} \! Z=\sum_{n=0}^{\infty } \frac{(40n+3)\left ( 
\frac{1}{2} \right )_n \left ( \frac{1}{4} \right )_n \left ( 
\frac{3}{4} \right )_n} {(n!)^3{49}^{2n+1}}\!
\pi=\frac{2\sqrt{11}}{11Z} \! Z=\sum_{n=0}^{\infty } \frac{(280n+19)\left ( 
\frac{1}{2} \right )_n \left ( \frac{1}{4} \right )_n \left ( 
\frac{3}{4} \right )_n} {(n!)^3{99}^{2n+1}}\!
\pi=\frac{\sqrt{2}}{4Z} \! Z=\sum_{n=0}^{\infty } \frac{(10n+1) \left ( 
\frac{1}{2} \right )_n \left ( \frac{1}{4} \right )_n \left ( 
\frac{3}{4} \right )_n} {(n!)^3{9}^{2n+1}}\!
\pi=\frac{4\sqrt{5}}{5Z}  \! Z=\sum_{n=0}^{\infty } \frac{(644n+41) \left (
 \frac{1}{2} \right )_n \left ( \frac{1}{4} \right )_n \left ( 
\frac{3}{4} \right )_n} {(n!)^35^n{72}^{2n+1}}\!
\pi=\frac{4\sqrt{3}}{3Z} \! Z=\sum_{n=0}^{\infty } \frac{(-1)^n(28n+3) 
\left ( \frac{1}{2} \right )_n \left ( \frac{1}{4} \right )_n \left ( 
\frac{3}{4} \right )_n} { (n!)^3{3^n}{4}^{n+1}}\!
 \pi=\frac{4}{Z}\! Z=\sum_{n=0}^{\infty } \frac{(-1)^n(20n+3) 
\left ( \frac{1}{2} \right )_n \left ( \frac{1}{4} \right )_n \left ( 
\frac{3}{4} \right )_n} { (n!)^3{2}^{2n+1}}\!
\pi=\frac{72}{Z} \! Z=\sum_{n=0}^{\infty } 
\pi=\frac{3528}{Z} \! Z=\sum_{n=0}^{\infty } 
(x)_n \!
is the pochhammer symbol for the falling factorial.
The related one found by the Chudnovsky brothers in 1987 is
\frac{426880 \sqrt{10005}}{\pi} = 
\sum_{k=0}^\infty \frac{(6k)! (13591409 + 545140134k)}{(3k)!(k!)^3 
which delivers 14 digits per term.[58] The Chudnovskys used this formula to set several π computing records in the end of the 1980s, including the first calculation of over one billion (1,011,196,691) decimals in 1989. It remains the formula of choice for π calculating software that runs on personal computers, as opposed to the supercomputers used to set modern records. On August 5, 2010, PhysOrg.com reported that Japanese and American computer experts Shigeru Kondo and Alexander Yee said they've calculated the value of π to 5 trillion decimal places on a personal computer, double the previous record.[61]
Whereas series typically increase the accuracy with a fixed amount for each added term, there exist iterative algorithms that multiply the number of correct digits at each step, with the downside that each step generally requires an expensive calculation. A breakthrough was made in 1975, when Richard Brent and Eugene Salamin independently discovered the Brent–Salamin algorithm, which uses only arithmetic to double the number of correct digits at each step.[62] The algorithm consists of setting
a_0 = 1 \quad \quad \quad b_0 = \frac{1}{\sqrt
 2} \quad \quad \quad t_0 = \frac{1}{4} \quad \quad \quad p_0 = 1\!
and iterating
a_{n+1} = \frac{a_n+b_n}{2} \quad \quad \quad 
b_{n+1} = \sqrt{a_n b_n}\!
t_{n+1} = t_n - p_n (a_n-a_{n+1})^2 \quad 
\quad \quad p_{n+1} = 2 p_n\!
until an and bn are close enough. Then the estimate for π is given by
\pi \approx \frac{(a_n + b_n)^2}{4 t_n}.\!
Using this scheme, 25 iterations suffice to reach 45 million correct decimals. A similar algorithm that quadruples the accuracy in each step has been found by Jonathan and Peter Borwein.[63] The methods have been used by Yasumasa Kanada and team to set most of the π calculation records since 1980, up to a calculation of 206,158,430,000 decimals of π in 1999. As of January 2010, the record was almost 2.7 trillion digits.[64] This beats the previous record of 2,576,980,370,000 decimals, set by Daisuke Takahashi on the T2K-Tsukuba System, a supercomputer at the University of Tsukuba northeast of Tokyo.[65]
An important recent development was the Bailey–Borwein–Plouffe formula (BBP formula), discovered by Simon Plouffe and named after the authors of the paper in which the formula was first published, David H. Bailey, Peter Borwein, and Simon Plouffe.[66] The formula,
\pi = \sum_{k=0}^\infty \frac{1}{16^k} \left( 
\frac{4}{8k + 1} - \frac{2}{8k + 4} - \frac{1}{8k + 5} - \frac{1}{8k + 
is remarkable because it allows extracting any individual hexadecimal or binary digit of π without calculating all the preceding ones.[66] Between 1998 and 2000, the distributed computing project PiHex used a modification of the BBP formula due to Fabrice Bellard to compute the quadrillionth (1,000,000,000,000,000:th) bit of π, which turned out to be 0.[67]
If a formula of the form
\pi = \sum_{k=0}^\infty \frac{1}{b^{ck}} 
were found where b and c are positive integers and p and q are polynomials with fixed degree and integer coefficients (as in the BPP formula above), this would be one the most efficient ways of computing any digit of π at any position in base bc without computing all the preceding digits in that base, in a time just depending on the size of the integer k and on the fixed degree of the polynomials. Plouffe also describes such formulas as the interesting ones for computing numbers of class SC*, in a logarithmically polynomial space and almost linear time, depending only on the size (order of magnitude) of the integer k, and requiring modest computing resources. The previous formula (found by Plouffe for π with b = 2 and c = 4, but also found for log(9/10) and for a few other irrational constants), implies that π is a SC* number.
In 2006, Simon Plouffe, using the integer relation algorithm PSLQ, found a series of formulas.[68] Let q = eπ (Gelfond's constant), then
\frac{\pi}{24} = \sum_{n=1}^\infty \frac{1}{n}
 \left(\frac{3}{q^n-1} - \frac{4}{q^{2n}-1} + \frac{1}{q^{4n}-1}\right)
\frac{\pi^3}{180} = \sum_{n=1}^\infty 
\frac{1}{n^3} \left(\frac{4}{q^n-1} - \frac{5}{q^{2n}-1} + 
and others of form,
\pi^k = \sum_{n=1}^\infty \frac{1}{n^k} 
\left(\frac{a}{q^n-1} + \frac{b}{q^{2n}-1} + \frac{c}{q^{4n}-1}\right)
where k is an odd number, and abc are rational numbers.
In the previous formula, if k is of the form 4m + 3, then the formula has the particularly simple form,
p\pi^k = \sum_{n=1}^\infty \frac{1}{n^k} 
\left(\frac{2^{k-1}}{q^n-1} - \frac{2^{k-1}+1}{q^{2n}-1} + 
for some rational number p where the denominator is a highly factorable number. General expressions for these kinds of sums are given in [69]

Pi and continued fraction

The sequence of partial denominators of the simple continued fraction of π does not show any obvious pattern:[30]
\pi=3+\textstyle \cfrac{1}{7+\textstyle 
\cfrac{1}{15+\textstyle \cfrac{1}{1+\textstyle \cfrac{1}{292+\textstyle 
\cfrac{1}{1+\textstyle \cfrac{1}{1+\textstyle \cfrac{1}{1+\ddots}}}}}}}
However, there are generalized continued fractions for π with a perfectly regular structure, such as:[70]
\pi=\textstyle \cfrac{4}{1+\textstyle 
\cfrac{1^2}{2+\textstyle \cfrac{3^2}{2+\textstyle 
\cfrac{5^2}{2+\textstyle \cfrac{7^2}{2+\textstyle 
=3+\textstyle \cfrac{1^2}{6+\textstyle \cfrac{3^2}{6+\textstyle 
\cfrac{5^2}{6+\textstyle \cfrac{7^2}{6+\textstyle 
=\textstyle \cfrac{4}{1+\textstyle \cfrac{1^2}{3+\textstyle 
\cfrac{2^2}{5+\textstyle \cfrac{3^2}{7+\textstyle 

Memorizing digits

Recent decades have seen a surge in the record for number of digits memorized.
Well before computers were used in calculating π, memorizing a record number of digits had become an obsession for some people. In 2006, Akira Haraguchi, a retired Japanese engineer, claimed to have recited 100,000 decimal places.[71] This, however, has yet to be verified by Guinness World Records. The Guinness-recognized record for remembered digits of π is 67,890 digits, held by Lu Chao, a 24-year-old graduate student from China.[72] It took him 24 hours and 4 minutes to recite to the 67,890th decimal place of π without an error.[73]
There are many ways to memorize π, including the use of "piems", which are poems that represent π in a way such that the length of each word (in letters) represents a digit. Here is an example of a piem, originally devised by Sir James Jeans: How I want a drink, alcoholic of course, after the heavy lectures involving quantum mechanics.[74][75] Notice how the first word has 3 letters, the second word has 1, the third has 4, the fourth has 1, the fifth has 5, and so on. The Cadaeic Cadenza contains the first 3835 digits of π in this manner.[76] Piems are related to the entire field of humorous yet serious study that involves the use of mnemonic techniques to remember the digits of π, known as piphilology. In other languages there are similar methods of memorization. However, this method proves inefficient for large memorizations of π. Other methods include remembering patterns in the numbers and the method of loci.[77][78]

Open questions

One open question about π is whether it is a normal number—whether any digit block occurs in the expansion of π just as often as one would statistically expect if the digits had been produced completely "randomly", and that this is true in every integer base, not just base 10.[79] Current knowledge on this point is very weak; e.g., it is not even known which of the digits 0,...,9 occur infinitely often in the decimal expansion of π,[80] although it is clear that at least two such digits must occur infinitely often, since otherwise π would be rational, which it is not.
Bailey and Crandall showed in 2000 that the existence of the above mentioned Bailey-Borwein-Plouffe formula and similar formulas imply that the normality in base 2 of π and various other constants can be reduced to a plausible conjecture of chaos theory.[81]
It is also unknown whether π and e are algebraically independent, although Yuri Nesterenko proved the algebraic independence of {π, eπ, Γ(1/4)} in 1996.[82]

Use in mathematics and science

π is ubiquitous in mathematics, science, and engineering.[83] It appears even in places that lack an obvious connection to the circles of Euclidean geometry.[84]

Geometry and trigonometry

For any circle with radius r and diameter d = 2r, the circumference is πd and the area is πr2. Further, π appears in formulas for areas and volumes of many other geometrical shapes based on circles, such as ellipses, spheres, cones, and tori.[85] Accordingly, π appears in definite integrals that describe circumference, area or volume of shapes generated by circles. In the basic case, half the area of the unit disk is given by the integral:[86] \int_{-1}^1 \sqrt{1-x^2}\,dx = \frac{\pi}{2} and
\int_{-1}^1\frac{1}{\sqrt{1-x^2}}\,dx = \pi
gives half the circumference of the unit circle.[85] More complicated shapes can be integrated as solids of revolution.[87]
From the unit-circle definition of the trigonometric functions also follows that the sine and cosine have period 2π. That is, for all x and integers n, sin(x) = sin(x + 2πn) and cos(x) = cos(x + 2πn). Because sin(0) = 0, sin(2πn) = 0 for all integers n. Also, the angle measure of 180° is equal to π radians. In other words, 1° = (π/180) radians.
In modern mathematics, π is often defined using trigonometric functions, for example as the smallest positive x for which sin x = 0, to avoid unnecessary dependence on the subtleties of Euclidean geometry and integration. Equivalently, π can be defined using the inverse trigonometric functions, for example as π = 2 arccos(0) or π = 4 arctan(1). Expanding inverse trigonometric functions as power series is the easiest way to derive infinite series for π.

Complex numbers and calculus

Euler's formula depicted on the complex plane. Increasing the angle φ to π radians (180°) yields Euler's identity.
A complex number z can be expressed in polar coordinates as follows:
z = r\cdot(\cos\varphi + i\sin\varphi)
The frequent appearance of π in complex analysis can be related to the behavior of the exponential function of a complex variable, described by Euler's formula
e^{i\varphi} = \cos \varphi + i\sin \varphi 
where i is the imaginary unit satisfying i2 = −1 and e ≈ 2.71828 is Euler's number. This formula implies that imaginary powers of e describe rotations on the unit circle in the complex plane; these rotations have a period of 360° = 2π. In particular, the 180° rotation φ = π results in the remarkable Euler's identity
e^{i \pi} = -1.\!
e^{i \pi} + 1 = 0.\!
Euler's identity is famous for linking several basic mathematical concepts in one concise, elegant expression.
There are n different n-th roots of unity
e^{2 \pi i k/n} \qquad (k = 0, 1, 2, \dots, n -
The Gaussian integral
A consequence is that the gamma function of a half-integer is a rational multiple of √π.


Although not a physical constant, π appears routinely in equations describing fundamental principles of the Universe, due in no small part to its relationship to the nature of the circle and, correspondingly, spherical coordinate systems. Using units such as Planck units can sometimes eliminate π from formulae.
 \Delta x\, \Delta p \ge \frac{h}{4\pi} = 
\frac \hbar 2
 R_{ik} - {g_{ik} R \over 2} + \Lambda g_{ik} =
 {8 \pi G \over c^4} T_{ik}
\Lambda = 8\pi G\ \rho_\text{vac}
 F = \frac{\left|q_1q_2\right|}{4 \pi 
\varepsilon_0 r^2}
 \mu_0 = 4 \pi \cdot 
\left(\frac{2\pi}{P}\right)^2 \, a^3 =\omega^2
  a^3=G (M+m)

Probability and statistics

In probability and statistics, there are many distributions whose formulas contain π, including:
f(x) = {1 \over \sigma\sqrt{2\pi} 
}\,e^{-(x-\mu )^2/(2\sigma^2)}
f(x) = \frac{1}{\pi (1 + x^2)}.
Note that since \int_{-\infty}^{\infty} 
f(x)\,dx = 1 for any probability density function f(x), the above formulas can be used to produce other integral formulas for π.[95]
Buffon's needle problem is sometimes quoted as an empirical approximation of π in "popular mathematics" works. Consider dropping a needle of length L repeatedly on a surface containing parallel lines drawn S units apart (with S > L). If the needle is dropped n times and x of those times it comes to rest crossing a line (x > 0), then one may approximate π using the Monte Carlo method:[96][97][98][99]
\pi \approx \frac{2nL}{xS}.
Though this result is mathematically impeccable, it cannot be used to determine more than very few digits of π by experiment. Reliably getting just three digits (including the initial "3") right requires millions of throws,[96] and the number of throws grows exponentially with the number of digits desired. Furthermore, any error in the measurement of the lengths L and S will transfer directly to an error in the approximated π. For example, a difference of a single atom in the length of a 10-centimeter needle would show up around the 9th digit of the result. In practice, uncertainties in determining whether the needle actually crosses a line when it appears to exactly touch it will limit the attainable accuracy to much less than 9 digits.[citation needed]

Geomorphology and chaos theory

Under ideal conditions (uniform gentle slope on an homogeneously erodible substrate), the ratio between the actual length of a river and its straight-line from source to mouth length tends to approach π.[100] Albert Einstein was the first to suggest that rivers have a tendency towards an ever more loopy path because the slightest curve will lead to faster currents on the outer side, which in turn will result in more erosion and a sharper bend. The sharper the bend, the faster the currents on the outer edge, the more the erosion, the more the river will twist and so on. However, increasing loopiness will result in rivers doubling back on themselves and effectively short-circuiting, creating an ox-bow lake. The balance between these two opposing factors leads to an average ratio of π between the actual length and the direct distance between source and mouth.[101]

In popular culture

A "Pi pie" to celebrate Pi Day
Nobel prize winning poet Wisława Szymborska wrote a poem about π, and here is an excerpt:[102]
The caravan of digits that is pi
does not stop at the edge of the page,
but runs off the table and into the air,
over the wall, a leaf, a bird's nest, the clouds, straight into the sky,
through all the bloatedness and bottomlessness.
Oh how short, all but mouse-like is the comet's tail!
Many schools around the world observe Pi Day (March 14, from 3.14).[103] At least one cheer at the Massachusetts Institute of Technology includes "3.14159!"[104]
On November 7, 2005, alternative musician Kate Bush released the album Aerial. The album contains the song "π" whose lyrics consist principally of Bush singing the digits of π to music, beginning with "3.14".[105]
In Carl Sagan's novel Contact, π played a key role in the story. The novel suggested that there was a message buried deep within the digits of π placed there by the creator of the universe.[106] This part of the story was omitted from the film adaptation of the novel.
Darren Aronofsky's film Pi deals with a number theorist.
The Wheel of Dublin (Ferris wheel) has been nicknamed "the pi in the sky".[107]
In the fictional movie, Night at the Museum: Battle of the Smithsonian, pi is the answer to the combination that will allow the Tablet of Akh-man-Ra to open the gates to the Underworld.