Search This Blog

Tuesday, September 9, 2014

Pi

Pi

From Wikipedia, the free encyclopedia

The number π is a mathematical constant, the ratio of a circle's circumference to its diameter, approximately equal to 3.14159. It has been represented by the Greek letter "π" since the mid-18th century though it is also sometimes spelled out as "pi" (/p/).

Being an irrational number, π cannot be expressed exactly as a common fraction, although fractions such as 22/7 and other rational numbers are commonly used to approximate π. Consequently its decimal representation never ends and never settles into a permanent repeating pattern. The digits appear to be randomly distributed although, to date, no proof of this has been discovered. Also, π is a transcendental number – a number that is not the root of any non-zero polynomial having rational coefficients. This transcendence of π implies that it is impossible to solve the ancient challenge of squaring the circle with a compass and straight-edge.

For thousands of years, mathematicians have attempted to extend their understanding of π, sometimes by computing its value to a high degree of accuracy. Before the 15th century mathematicians such as Archimedes and Liu Hui used geometrical techniques, based on polygons, to estimate the value of π. Starting around the 15th century, new algorithms based on infinite series revolutionized the computation of π. In the 20th and 21st centuries mathematicians and computer scientists discovered new approaches that, when combined with increasing computational power, extended the decimal representation of π to, as of late 2013, over 12 trillion (1013) digits.[1] Scientific applications generally require no more than 40 digits of π so the primary motivation for these computations is the human desire to break records. However, the extensive calculations involved have been used to test supercomputers and high-precision multiplication algorithms.

Because its definition relates to the circle, π is found in many formulae in trigonometry and geometry, especially those concerning circles, ellipses or spheres. It is also found in formulae used in other branches of science such as cosmology, number theory, statistics, fractals, thermodynamics, mechanics and electromagnetism. The ubiquity of π makes it one of the most widely-known mathematical constants both inside and outside the scientific community: Several books devoted to it have been published, the number is celebrated on Pi Day and record-setting calculations of the digits of π often result in news headlines. Attempts to memorize the value of π with increasing precision have led to records of over 67,000 digits.

Fundamentals

Name

The symbol used by mathematicians to represent the ratio of a circle's circumference to its diameter is the Greek letter π, sometimes spelled out as pi. In English, π is pronounced as "pie" ( /p/, /ˈpaɪ/).[2]
In mathematical use, the lower-case letter π (or π in sans-serif font) is distinguished from the capital letter Π, which denotes a product of a sequence.

The choice of the symbol π is discussed below.

Definition

A diagram of a circle, with the width labeled as diameter, and the perimeter labeled as circumference
The circumference of a circle is slightly more than three times as long as its diameter. The exact ratio is called π.

π is commonly defined as the ratio of a circle's circumference C to its diameter d:[3]
 \pi = \frac{C}{d}
The ratio C/d is constant, regardless of the circle's size. For example, if a circle has twice the diameter of another circle it will also have twice the circumference, preserving the ratio C/d. This definition of π implicitly makes use of flat (Euclidean) geometry; although the notion of a circle can be extended to any curved (non-Euclidean) geometry, these new circles will no longer satisfy the formula π = C/d.[3] There are also other definitions of π that do not immediately involve circles at all. For example, π is twice the smallest positive x for which cos(x) equals 0.[3][4]

Properties

π is an irrational number, meaning that it cannot be written as the ratio of two integers (fractions such as 22/7 are commonly used to approximate π; no common fraction (ratio of whole numbers) can be its exact value).[5] Since π is irrational, it has an infinite number of digits in its decimal representation, and it does not end with an infinitely repeating pattern of digits. There are several proofs that π is irrational; they generally require calculus and rely on the reductio ad absurdum technique. The degree to which π can be approximated by rational numbers (called the irrationality measure) is not precisely known; estimates have established that the irrationality measure is larger than the measure of e or ln(2) but smaller than the measure of Liouville numbers.[6]
A diagram of a square and circle, both with identical area; the length of the side of the square is the square root of pi
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.

π is a transcendental number, which means that it is not the solution of any non-constant polynomial with rational coefficients, such as \scriptstyle \frac{x^5}{120}\,-\,\frac{x^3}{6}\,+\,x\,=\,0.[7][8] The transcendence of π has two important consequences: First, π cannot be expressed using any combination of rational numbers and square roots or n-th roots such as \scriptstyle \sqrt[3]{31} or \scriptstyle \sqrt[2]{10}. Second, since no transcendental number can be constructed with compass and straightedge, it is not possible to "square the circle". In other words, it is impossible to construct, using compass and straightedge alone, a square whose area is equal to the area of a given circle.[9] Squaring a circle was one of the important geometry problems of the classical antiquity.[10] Amateur mathematicians in modern times have sometimes attempted to square the circle and sometimes claim success despite the fact that it is impossible.[11]

The digits of π have no apparent pattern and have passed tests for statistical randomness, including tests for normality; a number of infinite length is called normal when all possible sequences of digits (of any given length) appear equally often.[12] The hypothesis that π is normal has not been proven or disproven.[12] Since the advent of computers, a large number of digits of π have been available on which to perform statistical analysis. Yasumasa Kanada has performed detailed statistical analyses on the decimal digits of π and found them consistent with normality; for example, the frequency of the ten digits 0 to 9 were subjected to statistical significance tests, and no evidence of a pattern was found.[13] Despite the fact that π's digits pass statistical tests for randomness, π contains some sequences of digits that may appear non-random to non-mathematicians, such as the Feynman point, which is a sequence of six consecutive 9s that begins at the 762nd decimal place of the decimal representation of π.[14]

Continued fractions

A photograph of the Greek letter pi, created as a large stone mosaic embedded in the ground.
The constant π is represented in this mosaic outside the mathematics building at the Technische Universität Berlin.

Like all irrational numbers, π cannot be represented as a common fraction (also known as a simple or vulgar fraction), by the very definition of "irrational". But every irrational number, including π, can be represented by an infinite series of nested fractions, called a continued fraction:

\pi=3+\textstyle \frac{1}{7+\textstyle \frac{1}{15+\textstyle \frac{1}{1+\textstyle \frac{1}{292+\textstyle \frac{1}{1+\textstyle \frac{1}{1+\textstyle \frac{1}{1+\ddots}}}}}}}
 Truncating the continued fraction at any point generates a fraction that provides an approximation for π; two such fractions (22/7 and 355/113) have been used historically to approximate the constant. Each approximation generated in this way is a best rational approximation; that is, each is closer to π than any other fraction with the same or a smaller denominator.[15] Because π is known to be transcendental, it is by definition not algebraic and so cannot be a quadratic irrational. Therefore π cannot have a periodic continued fraction. Although the simple continued fraction for π (shown above) also does not exhibit any other obvious pattern,[16] mathematicians have discovered several generalized continued fractions that do, such as:[17]
\pi=\textstyle \cfrac{4}{1+\textstyle \frac{1^2}{2+\textstyle \frac{3^2}{2+\textstyle \frac{5^2}{2+\textstyle \frac{7^2}{2+\textstyle \frac{9^2}{2+\ddots}}}}}}
=3+\textstyle \frac{1^2}{6+\textstyle \frac{3^2}{6+\textstyle \frac{5^2}{6+\textstyle \frac{7^2}{6+\textstyle \frac{9^2}{6+\ddots}}}}}
=\textstyle \cfrac{4}{1+\textstyle \frac{1^2}{3+\textstyle \frac{2^2}{5+\textstyle \frac{3^2}{7+\textstyle \frac{4^2}{9+\ddots}}}}}

Approximate value

Fractions: Approximate fractions include (in order of increasing accuracy) 22/7, 333/106, 355/113, 52163/16604, 103993/33102, and 245850922/78256779.[15] (List is selected terms from OEISA063674 and OEISA063673.)
  • Decimal: The first 50 decimal digits are 3.14159265358979323846264338327950288419716939937510...[18] OEISA000796
  • Binary: The base 2 approximation to 48 digits is 11.001001000011111101101010100010001000010110100011...
  • Hexadecimal: The base 16 approximation to 20 digits is 3.243F6A8885A308D31319...[19]
  • Sexagesimal: A base 60 approximation to four sexagesimal digits is 3;8,29,44,1

History

Antiquity

The Great Pyramid at Giza, constructed c. 2589–2566 BC, was built with a perimeter of about 1760 cubits and a height of about 280 cubits; the ratio 1760/280 ≈ 6.2857 is approximately equal to 2π ≈ 6.2832. Based on this ratio, some Egyptologists concluded that the pyramid builders had knowledge of π and deliberately designed the pyramid to incorporate the proportions of a circle.[20]
Others maintain that the suggested relationship to π is merely a coincidence, because there is no evidence that the pyramid builders had any knowledge of π, and because the dimensions of the pyramid are based on other factors.[21]

The earliest written approximations of π are found in Egypt and Babylon, both within 1 percent of the true value. In Babylon, a clay tablet dated 1900–1600 BC has a geometrical statement that, by implication, treats π as 25/8 = 3.1250.[22] In Egypt, the Rhind Papyrus, dated around 1650 BC, but copied from a document dated to 1850 BC has a formula for the area of a circle that treats π as (16/9)2 ≈ 3.1605.[22]

In India around 600 BC, the Shulba Sutras (Sanskrit texts that are rich in mathematical contents) treat π as (9785/5568)2 ≈ 3.088.[23] In 150 BC, or perhaps earlier, Indian sources treat π as \scriptstyle \sqrt{10} ≈ 3.1622.[24]

Two verses in the Hebrew Bible (written between the 8th and 3rd centuries BC) describe a ceremonial pool in the Temple of Solomon with a diameter of ten cubits and a circumference of thirty cubits; the verses imply π is about three if the pool is circular.[25][26] Rabbi Nehemiah explained the discrepancy as being due to the thickness of the vessel. His early work of geometry, Mishnat ha-Middot, was written around 150 AD and takes the value of π to be three and one seventh.[27] See Approximations of π.

Polygon approximation era

diagram of a hexagon and pentagon circumscribed outside a circle
π can be estimated by computing the perimeters of circumscribed and inscribed polygons.

The first recorded algorithm for rigorously calculating the value of π was a geometrical approach using polygons, devised around 250 BC by the Greek mathematician Archimedes.[28] This polygonal algorithm dominated for over 1,000 years, and as a result π is sometimes referred to as "Archimedes' constant".[29] Archimedes computed upper and lower bounds of π by drawing a regular hexagon inside and outside a circle, and successively doubling the number of sides until he reached a 96-sided regular polygon. By calculating the perimeters of these polygons, he proved that 223/71 < π < 22/7 (3.1408 < π < 3.1429).[30] Archimedes' upper bound of 22/7 may have led to a widespread popular belief that π is equal to 22/7.[31] Around 150 AD, Greek-Roman scientist Ptolemy, in his Almagest, gave a value for π of 3.1416, which he may have obtained from Archimedes or from Apollonius of Perga.[32] Mathematicians using polygonal algorithms reached 39 digits of π in 1630, a record only broken in 1699 when infinite series were used to reach 71 digits.[33]
A painting of a man studying
Archimedes developed the polygonal approach to approximating π.

In ancient China, values for π included 3.1547 (around 1 AD), \scriptstyle \sqrt{10} (100 AD, approximately 3.1623), and 142/45 (3rd century, approximately 3.1556).[34] Around 265 AD, the Wei Kingdom mathematician Liu Hui created a polygon-based iterative algorithm and used it with a 3,072-sided polygon to obtain a value of π of 3.1416.[35][36] Liu later invented a faster method of calculating π and obtained a value of 3.14 with a 96-sided polygon, by taking advantage of the fact that the differences in area of successive polygons form a geometric series with a factor of 4.[35] The Chinese mathematician Zu Chongzhi, around 480 AD, calculated that π ≈ 355/113 (a fraction that goes by the name Milü in Chinese), using Liu Hui's algorithm applied to a 12,288-sided polygon. With a correct value for its seven first decimal digits, this value of 3.141592920... remained the most accurate approximation of π available for the next 800 years.[37]

The Indian astronomer Aryabhata used a value of 3.1416 in his Āryabhaṭīya (499 AD).[38] Fibonacci in c. 1220 computed 3.1418 using a polygonal method, independent of Archimedes.[39] Italian author Dante apparently employed the value \scriptstyle 3+\sqrt{2}/10 ≈ 3.14142.[39]

The Persian astronomer Jamshīd al-Kāshī produced 16 digits in 1424 using a polygon with 3×228 sides,[40][41] which stood as the world record for about 180 years.[42] French mathematician François Viète in 1579 achieved 9 digits with a polygon of 3×217 sides.[42] Flemish mathematician Adriaan van Roomen arrived at 15 decimal places in 1593.[42] In 1596, Dutch mathematician Ludolph van Ceulen reached 20 digits, a record he later increased to 35 digits (as a result, π was called the "Ludolphian number" in Germany until the early 20th century).[43] Dutch scientist Willebrord Snellius reached 34 digits in 1621,[44] and Austrian astronomer Christoph Grienberger arrived at 38 digits in 1630,[45] which remains the most accurate approximation manually achieved using polygonal algorithms.[44]

Infinite series

Comparison of the convergence of several historical infinite series for π. Sn is the approximation after taking n terms. Each subsequent subplot magnifies the shaded area horizontally by 10 times. (click for detail)

The calculation of π was revolutionized by the development of infinite series techniques in the 16th and 17th centuries. An infinite series is the sum of the terms of an infinite sequence.[46] Infinite series allowed mathematicians to compute π with much greater precision than Archimedes and others who used geometrical techniques.[46] Although infinite series were exploited for π most notably by European mathematicians such as James Gregory and Gottfried Wilhelm Leibniz, the approach was first discovered in India sometime between 1400 and 1500 AD.[47] The first written description of an infinite series that could be used to compute π was laid out in Sanskrit verse by Indian astronomer Nilakantha Somayaji in his Tantrasamgraha, around 1500 AD.[48] The series are presented without proof, but proofs are presented in a later Indian work, Yuktibhāṣā, from around 1530 AD. Nilakantha attributes the series to an earlier Indian mathematician, Madhava of Sangamagrama, who lived c. 1350 – c. 1425.[48] Several infinite series are described, including series for sine, tangent, and cosine, which are now referred to as the Madhava series or Gregory–Leibniz series.[48] Madhava used infinite series to estimate π to 11 digits around 1400, but that value was improved on around 1430 by the Persian mathematician Jamshīd al-Kāshī, using a polygonal algorithm.[49]
A formal portrait of a man, with long hair
Isaac Newton used infinite series to compute π to 15 digits, later writing "I am ashamed to tell you to how many figures I carried these computations".[50]

The first infinite sequence discovered in Europe was an infinite product (rather than an infinite sum, which are more typically used in π calculations) found by French mathematician François Viète in 1593:[51]
 \frac2\pi = \frac{\sqrt2}2 \cdot \frac{\sqrt{2+\sqrt2}}2 \cdot \frac{\sqrt{2+\sqrt{2+\sqrt2}}}2 \cdots
The second infinite sequence found in Europe, by John Wallis in 1655, was also an infinite product.[51] The discovery of calculus, by English scientist Isaac Newton and German mathematician Gottfried Wilhelm Leibniz in the 1660s, led to the development of many infinite series for approximating π. Newton himself used an arcsin series to compute a 15 digit approximation of π in 1665 or 1666, later writing "I am ashamed to tell you to how many figures I carried these computations, having no other business at the time."[50]

In Europe, Madhava's formula was rediscovered by Scottish mathematician James Gregory in 1671, and by Leibniz in 1674:[52][53]

\arctan z = z - \frac {z^3} {3} +\frac {z^5} {5} -\frac {z^7} {7} +\cdots
This formula, the Gregory–Leibniz series, equals \scriptstyle \pi/4 when evaluated with z = 1.[53] In 1699, English mathematician Abraham Sharp used the Gregory–Leibniz series to compute π to 71 digits, breaking the previous record of 39 digits, which was set with a polygonal algorithm.[54] The Gregory–Leibniz series is simple, but converges very slowly (that is, approaches the answer gradually), so it is not used in modern π calculations.[55]

In 1706 John Machin used the Gregory–Leibniz series to produce an algorithm that converged much faster:[56]
 \frac{\pi}{4} = 4 \, \arctan \frac{1}{5} - \arctan \frac{1}{239}
Machin reached 100 digits of π with this formula.[57] Other mathematicians created variants, now known as Machin-like formulae, that were used to set several successive records for calculating digits of π.[57] Machin-like formulae remained the best-known method for calculating π well into the age of computers, and were used to set records for 250 years, culminating in a 620-digit approximation in 1946 by Daniel Ferguson – the best approximation achieved without the aid of a calculating device.[58]

A 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 German mathematician Carl Friedrich Gauss.[59] British mathematician William Shanks famously took 15 years to calculate π to 707 digits, but made a mistake in the 528th digit, rendering all subsequent digits incorrect.[59]

Rate of convergence

Some infinite series for π converge faster than others. Given the choice of two infinite series for π, mathematicians will generally use the one that converges more rapidly because faster convergence reduces the amount of computation needed to calculate π to any given accuracy.[60] A simple infinite series for π is the Gregory–Leibniz series:[61]
 \pi = \frac{4}{1} - \frac{4}{3} + \frac{4}{5} - \frac{4}{7} + \frac{4}{9} - \frac{4}{11} + \frac{4}{13} - \cdots
As individual terms of this infinite series are added to the sum, the total gradually gets closer to π, and – with a sufficient number of terms – can get as close to π as desired. It converges quite slowly, though – after 500,000 terms, it produces only five correct decimal digits of π.[62]

An infinite series for π (published by Nilakantha in the 15th century) that converges more rapidly than the Gregory–Leibniz series is:[63]
 \pi = 3 + \frac{4}{2\times3\times4} - \frac{4}{4\times5\times6} + \frac{4}{6\times7\times8} - \frac{4}{8\times9\times10} + \cdots
The following table compares the convergence rates of these two series:

Infinite series for π After 1st term After 2nd term After 3rd term After 4th term After 5th term Converges to:
\scriptstyle \pi = \frac{4}{1} - \frac{4}{3} + \frac{4}{5} - \frac{4}{7} + \frac{4}{9} - \frac{4}{11} + \frac{4}{13} \cdots. 4.0000 2.6666... 3.4666... 2.8952... 3.3396... π = 3.1415...
\scriptstyle \pi = {{3}} + \frac{{4}}{2\times3\times4} - \frac{{4}}{4\times5\times6} + \frac{{4}}{6\times7\times8} \cdots. 3.0000 3.1666... 3.1333... 3.1452... 3.1396...

After five terms, the sum of the Gregory–Leibniz series is within 0.2 of the correct value of π, whereas the sum of Nilakantha's series is within 0.002 of the correct value of π. Nilakantha's series converges faster and is more useful for computing digits of π. Series that converge even faster include Machin's series and Chudnovsky's series, the latter producing 14 correct decimal digits per term.[60]

Irrationality and transcendence

Not all mathematical advances relating to π were aimed at increasing the accuracy of approximations. When Euler solved the Basel problem in 1735, finding the exact value of the sum of the reciprocal squares, he established a connection between π and the prime numbers that later contributed to the development and study of the Riemann zeta function:[64]
 \frac{\pi^2}{6} = \frac{1}{1^2} + \frac{1}{2^2} + \frac{1}{3^2} + \frac{1}{4^2} + \cdots
Swiss scientist Johann Heinrich Lambert in 1761 proved that π is irrational, meaning it is not equal to the quotient of any two whole numbers.[5] Lambert's proof exploited a continued-fraction representation of the tangent function.[65] French mathematician Adrien-Marie Legendre proved in 1794 that π2 is also irrational. In 1882, German mathematician Ferdinand von Lindemann proved that π is transcendental, confirming a conjecture made by both Legendre and Euler.[66]

Adoption of the symbol π

Leonhard Euler popularized the use of the Greek letter π in works he published in 1736 and 1748.

The earliest known use of the Greek letter π to represent the ratio of a circle's circumference to its diameter was by mathematician William Jones in his 1706 work Synopsis Palmariorum Matheseos; or, a New Introduction to the Mathematics.[67] The Greek letter first appears there in the phrase "1/2 Periphery (π)" in the discussion of a circle with radius one. Jones may have chosen π because it was the first letter in the Greek spelling of the word periphery.[68] However, he writes that his equations for π are from the "ready pen of the truly ingenious Mr. John Machin", leading to speculation that Machin may have employed the Greek letter before Jones.[69] It had indeed been used earlier for geometric concepts.[69] William Oughtred used π and δ, the Greek letter equivalents of p and d, to express ratios of periphery and diameter in the 1647 and later editions of Clavis Mathematicae.

After Jones introduced the Greek letter in 1706, it was not adopted by other mathematicians until Euler started using it, beginning with his 1736 work Mechanica. Before then, mathematicians sometimes used letters such as c or p instead.[69] Because Euler corresponded heavily with other mathematicians in Europe, the use of the Greek letter spread rapidly.[69] In 1748, Euler used π in his widely read work Introductio in analysin infinitorum (he wrote: "for the sake of brevity we will write this number as π; thus π is equal to half the circumference of a circle of radius 1") and the practice was universally adopted thereafter in the Western world.[69]

Modern quest for more digits

Computer era and iterative algorithms

Formal photo of a balding man wearing a suit
John von Neumann was part of the team that first used a digital computer, ENIAC, to compute π.


The development of computers in the mid-20th century again revolutionized the hunt for digits of π. American mathematicians John Wrench and Levi Smith reached 1,120 digits in 1949 using a desk calculator.[70] Using an inverse tangent (arctan) infinite series, a team led by George Reitwiesner and John von Neumann that same year achieved 2,037 digits with a calculation that took 70 hours of computer time on the ENIAC computer.[71] The record, always relying on an arctan series, was broken repeatedly (7,480 digits in 1957; 10,000 digits in 1958; 100,000 digits in 1961) until 1 million digits were reached in 1973.[72]

Two additional developments around 1980 once again accelerated the ability to compute π. First, the discovery of new iterative algorithms for computing π, which were much faster than the infinite series; and second, the invention of fast multiplication algorithms that could multiply large numbers very rapidly.[73] Such algorithms are particularly important in modern π computations, because most of the computer's time is devoted to multiplication.[74] They include the Karatsuba algorithm, Toom–Cook multiplication, and Fourier transform-based methods.[75]

The iterative algorithms were independently published in 1975–1976 by American physicist Eugene Salamin and Australian scientist Richard Brent.[76] These avoid reliance on infinite series. An iterative algorithm repeats a specific calculation, each iteration using the outputs from prior steps as its inputs, and produces a result in each step that converges to the desired value. The approach was actually invented over 160 years earlier by Carl Friedrich Gauss, in what is now termed the arithmetic–geometric mean method (AGM method) or Gauss–Legendre algorithm.[76] As modified by Salamin and Brent, it is also referred to as the Brent–Salamin algorithm.

The iterative algorithms were widely used after 1980 because they are faster than infinite series algorithms: whereas infinite series typically increase the number of correct digits additively in successive terms, iterative algorithms generally multiply the number of correct digits at each step. For example, the Brent-Salamin algorithm doubles the number of digits in each iteration. In 1984, the Canadian brothers John and Peter Borwein produced an iterative algorithm that quadruples the number of digits in each step; and in 1987, one that increases the number of digits five times in each step.[77] Iterative methods were used by Japanese mathematician Yasumasa Kanada to set several records for computing π between 1995 and 2002.[78] This rapid convergence comes at a price: the iterative algorithms require significantly more memory than infinite series.[78]

Motivations for computing π

As mathematicians discovered new algorithms, and computers became available, the number of known decimal digits of π increased dramatically. Note that the vertical scale is logarithmic.

For most numerical calculations involving π, a handful of digits provide sufficient precision. According to Jörg Arndt and Christoph Haenel, thirty-nine digits are sufficient to perform most cosmological calculations, because that is the accuracy necessary to calculate the volume of the known universe with a precision of one atom.[79] Despite this, people have worked strenuously to compute π to thousands and millions of digits.[80] This effort may be partly ascribed to the human compulsion to break records, and such achievements with π often make headlines around the world.[81][82] They also have practical benefits, such as testing supercomputers, testing numerical analysis algorithms (including high-precision multiplication algorithms); and within pure mathematics itself, providing data for evaluating the randomness of the digits of π.[83]

Rapidly convergent series

Photo portrait of a man
Srinivasa Ramanujan, working in isolation in India, produced many innovative series for computing π.

Modern π calculators do not use iterative algorithms exclusively. New infinite series were discovered in the 1980s and 1990s that are as fast as iterative algorithms, yet are simpler and less memory intensive.[78] The fast iterative algorithms were anticipated in 1914, when the Indian mathematician Srinivasa Ramanujan published dozens of innovative new formulae for π, remarkable for their elegance, mathematical depth, and rapid convergence.[84] One of his formulae, based on modular equations, is
\frac{1}{\pi} = \frac{2 \sqrt 2}{9801} \sum_{k=0}^\infty \frac{(4k)!(1103+26390k)}{k!^4(396^{4k})}.
This series converges much more rapidly than most arctan series, including Machin's formula.[85] Bill Gosper was the first to use it for advances in the calculation of π, setting a record of 17 million digits in 1985.[86] Ramanujan's formulae anticipated the modern algorithms developed by the Borwein brothers and the Chudnovsky brothers.[87] The Chudnovsky formula developed in 1987 is
\frac{1}{\pi} = \frac{12}{640320^{3/2}} \sum_{k=0}^\infty \frac{(6k)! (13591409 + 545140134k)}{(3k)!(k!)^3 (-640320)^{3k}}.
It produces about 14 digits of π per term,[88] and has been used for several record-setting π calculations, including the first to surpass 1 billion (109) digits in 1989 by the Chudnovsky brothers, 2.7 trillion (2.7×1012) digits by Fabrice Bellard in 2009, and 10 trillion (1013) digits in 2011 by Alexander Yee and Shigeru Kondo.[89][1] For similar formulas, see also the Ramanujan–Sato series.

In 2006, Canadian mathematician Simon Plouffe used the PSLQ integer relation algorithm[90] to generate several new formulas for π, conforming to the following template:
\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 \mathit{q} is eπ (Gelfond's constant),  \mathit{k} is an odd number, and \mathit{a, b, c} are certain rational numbers that Plouffe computed.[91]

Spigot algorithms

Two algorithms were discovered in 1995 that opened up new avenues of research into π. They are called spigot algorithms because, like water dripping from a spigot, they produce single digits of π that are not reused after they are calculated.[92][93] This is in contrast to infinite series or iterative algorithms, which retain and use all intermediate digits until the final result is produced.[92]

American mathematicians Stan Wagon and Stanley Rabinowitz produced a simple spigot algorithm in 1995.[93][94][95] Its speed is comparable to arctan algorithms, but not as fast as iterative algorithms.[94]

Another spigot algorithm, the BBP digit extraction algorithm, was discovered in 1995 by Simon Plouffe:[96][97]
 \pi = \sum_{i=0}^\infty \frac{1}{16^i} \left( \frac{4}{8i + 1} - \frac{2}{8i + 4} - \frac{1}{8i + 5} - \frac{1}{8i + 6}\right)
This formula, unlike others before it, can produce any individual hexadecimal digit of π without calculating all the preceding digits.[96] Individual binary digits may be extracted from individual hexadecimal digits, and octal digits can be extracted from one or two hexadecimal digits. Variations of the algorithm have been discovered, but no digit extraction algorithm has yet been found that rapidly produces decimal digits.[98] An important application of digit extraction algorithms is to validate new claims of record π computations: After a new record is claimed, the decimal result is converted to hexadecimal, and then a digit extraction algorithm is used to calculate several random hexadecimal digits near the end; if they match, this provides a measure of confidence that the entire computation is correct.[1]

Between 1998 and 2000, the distributed computing project PiHex used Bellard's formula (a modification of the BBP algorithm) to compute the quadrillionth (1015th) bit of π, which turned out to be 0.[99] In September 2010, a Yahoo! employee used the company's Hadoop application on one thousand computers over a 23-day period to compute 256 bits of π at the two-quadrillionth (2×1015th) bit, which also happens to be zero.[100]

Use

Because π is closely related to the circle, it is found in many formulae from the fields of geometry and trigonometry, particularly those concerning circles, spheres, or ellipses. Formulae from other branches of science also include π in some of their important formulae, including sciences such as statistics, fractals, thermodynamics, mechanics, cosmology, number theory, and electromagnetism.

Geometry and trigonometry

A diagram of a circle with a square coving the circle's upper right quadrant.
The area of the circle equals π times the shaded area.

π appears in formulae for areas and volumes of geometrical shapes based on circles, such as ellipses, spheres, cones, and tori. Below are some of the more common formulae that involve π.[101]
  • The circumference of a circle with radius r is  2 \pi r.
  • The area of a circle with radius r is  \pi r^2.
  • The volume of a sphere with radius r is  \tfrac43\pi r^3.
  • The surface area of a sphere with radius r is  4 \pi r^2.
The formulae above are special cases of the surface area S_n(r) and volume V_n(r) of an n-dimensional sphere.

S_n(r) = \frac{n\pi^{n/2}}{\Gamma(\frac{n}{2}+1)}r^{n-1}
V_n(r) = \frac{\pi^{n/2}}{\Gamma(\frac{n}{2}+1)}r^n

π appears in definite integrals that describe circumference, area, or volume of shapes generated by circles. For example, an integral that specifies half the area of a circle of radius one is given by:[102]
\int_{-1}^1 \sqrt{1-x^2}\,dx = \frac{\pi}{2}.
In that integral the function \scriptstyle \sqrt{1-x^2} represents the top half of a circle (the square root is a consequence of the Pythagorean theorem), and the integral \scriptstyle \int_{-1}^1 computes the area between that half of a circle and the x axis.
Diagram showing graphs of functions
Sine and cosine functions repeat with period 2π.

The trigonometric functions rely on angles, and mathematicians generally use radians as units of measurement. π plays an important role in angles measured in radians, which are defined so that a complete circle spans an angle of 2π radians.[103] The angle measure of 180° is equal to π radians, and 1° = π/180 radians.[103]

Common trigonometric functions have periods that are multiples of π; for example, sine and cosine have period 2π,[104] so for any angle θ and any integer k, \scriptstyle \sin\theta = \sin\left(\theta + 2\pi k \right) and \scriptstyle \cos\theta = \cos\left(\theta + 2\pi k \right).[104]

Monte Carlo methods

Needles of length ℓ scattered on stripes with width t
Buffon's needle. Needles a and b are dropped randomly.
Thousands of dots randomly covering a square and a circle inscribed in the square.
Random dots are placed on the quadrant of a square with a circle inscribed in it.
Monte Carlo methods, based on random trials, can be used to approximate π.

Monte Carlo methods, which evaluate the results of multiple random trials, can be used to create approximations of π.[105] Buffon's needle is one such technique: If a needle of length is dropped n times on a surface on which parallel lines are drawn t units apart, and if x of those times it comes to rest crossing a line (x > 0), then one may approximate π based on the counts:[106]
\pi \approx \frac{2n\ell}{xt}
Another Monte Carlo method for computing π is to draw a circle inscribed in a square, and randomly place dots in the square. The ratio of dots inside the circle to the total number of dots will approximately equal \scriptstyle \pi/4.[107]

Monte Carlo methods for approximating π are very slow compared to other methods, and are never used to approximate π when speed or accuracy are desired.[108]

Complex numbers and analysis

A diagram of a unit circle centered at the origin in the complex plane, including a ray from the center of the circle to its edge, with the triangle legs labeled with sine and cosine functions.
The association between imaginary powers of the number e and points on the unit circle centered at the origin in the complex plane given by Euler's formula.

Any complex number, say z, can be expressed using a pair of real numbers. In the polar coordinate system, one number (radius or r) is used to represent z's distance from the origin of the complex plane and the other (angle or φ) to represent a counter-clockwise rotation from the positive real line as follows:[109]
z = r\cdot(\cos\varphi + i\sin\varphi),
where i is the imaginary unit satisfying i2 = −1. 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:[110]
e^{i\varphi} = \cos \varphi + i\sin \varphi,
where the constant e is the base of the natural logarithm. This formula establishes a correspondence between imaginary powers of e and points on the unit circle centered at the origin of the complex plane. Setting φ = π in Euler's formula results in Euler's identity, celebrated by mathematicians because it contains the five most important mathematical constants:[110][111]
e^{i \pi} + 1 = 0.
There are n different complex numbers z satisfying  z^n = 1, and these are called the "n-th roots of unity".[112] They are given by this formula:
e^{2 \pi i k/n} \qquad (k = 0, 1, 2, \dots, n - 1).
Cauchy's integral formula governs complex analytic functions and establishes an important relationship between integration and differentiation, including the fact that the values of a complex function within a closed boundary are entirely determined by the values on the boundary:[113][114]
f (z_{0}) = \frac{1}{ 2\pi i } \oint_\gamma { f(z) \over z-z_0 }\,dz
An complex black shape on a blue background.
π can be computed from the Mandelbrot set, by counting the number of iterations required before point (−0.75, ε) diverges.

An occurrence of π in the Mandelbrot set fractal was discovered by American David Boll in 1991.[115] He examined the behavior of the Mandelbrot set near the "neck" at (−0.75, 0). If points with coordinates (−0.75, ε) are considered, as ε tends to zero, the number of iterations until divergence for the point multiplied by ε converges to π. The point (0.25, ε) at the cusp of the large "valley" on the right side of the Mandelbrot set behaves similarly: the number of iterations until divergence multiplied by the square root of ε tends to π.[115][116]

The gamma function extends the concept of factorial – which is normally defined only for whole numbers – to all real numbers. When the gamma function is evaluated at half-integers, the result contains π; for example \scriptstyle \Gamma(1/2) = \sqrt{\pi} and \scriptstyle\Gamma(5/2) = \frac {3 \sqrt{\pi}} {4} .[117] The gamma function can be used to create a simple approximation to \scriptstyle n! for large \scriptstyle n: \scriptstyle n! \sim \sqrt{2 \pi n} \left(\frac{n}{e}\right)^n which is known as Stirling's approximation.[118]

Number theory and Riemann zeta function

The Riemann zeta function ζ(s) is used in many areas of mathematics. When evaluated at \scriptstyle s = 2 it can be written as
 \zeta(2) = \frac{1}{1^2} + \frac{1}{2^2} + \frac{1}{3^2} + \cdots
Finding a simple solution for this infinite series was a famous problem in mathematics called the Basel problem. Leonhard Euler solved it in 1735 when he showed it was equal to \scriptstyle \pi^2/6.[64] Euler's result leads to the number theory result that the probability of two random numbers being relatively prime (that is, having no shared factors) is equal to \scriptstyle 6/ \pi^2.[119][120] This probability is based on the observation that the probability that any number is divisible by a prime \scriptstyle p is \scriptstyle 1/p (for example, every 7th integer is divisible by 7.) Hence the probability that two numbers are both divisible by this prime is \scriptstyle 1/p^2, and the probability that at least one of them is not is \scriptstyle 1-1/p^2. For distinct primes, these divisibility events are mutually independent; so the probability that two numbers are relatively prime is given by a product over all primes:[121]
\prod_p^{\infty} \left(1-\frac{1}{p^2}\right) = \left( \prod_p^{\infty} \frac{1}{1-p^{-2}} \right)^{-1} = \frac{1}{1 + \frac{1}{2^2} + \frac{1}{3^2} + \cdots } = \frac{1}{\zeta(2)} = \frac{6}{\pi^2} \approx 61\%
This probability can be used in conjunction with a random number generator to approximate π using a Monte Carlo approach.[122]

Probability and statistics

A graph of the Gaussian function
ƒ(x) = ex2. The colored region between the function and the x-axis has area  \scriptstyle\sqrt{\pi} .

The fields of probability and statistics frequently use the normal distribution as a simple model for complex phenomena; for example, scientists generally assume that the observational error in most experiments follows a normal distribution.[123] π is found in the Gaussian function (which is the probability density function of the normal distribution) with mean μ and standard deviation σ:[124]
f(x) = {1 \over \sigma\sqrt{2\pi} }\,e^{-(x-\mu )^2/(2\sigma^2)}
The area under the graph of the normal distribution curve is given by the Gaussian integral:[124]
\int_{-\infty}^\infty e^{-x^2} \, dx=\sqrt{\pi},
while the related integral for the Cauchy distribution is
\int_{-\infty }^{\infty } \frac{1}{x^2+1} \, dx = \pi.

Outside mathematics

Describing physical phenomena

Although not a physical constant, π appears routinely in equations describing fundamental principles of the universe, often because of π's relationship to the circle and to spherical coordinate systems. A simple formula from the field of classical mechanics gives the approximate period T of a simple pendulum of length L, swinging with a small amplitude (g is the earth's gravitational acceleration):[125]
T \approx 2\pi \sqrt\frac{L}{g}.
One of the key formulae of quantum mechanics is Heisenberg's uncertainty principle, which shows that the uncertainty in the measurement of a particle's position (Δx) and momentump) cannot both be arbitrarily small at the same time (where h is Planck's constant):[126]
 \Delta x\, \Delta p \ge \frac{h}{4\pi}.
In the domain of cosmology, π appears in Einstein's field equation, a fundamental formula which forms the basis of the general theory of relativity and describes the fundamental interaction of gravitation as a result of spacetime being curved by matter and energy:[127]
 R_{ik} - {g_{ik} R \over 2} + \Lambda g_{ik} = {8 \pi G \over c^4} T_{ik},
where R_{ik}\, is the Ricci curvature tensor, R\, is the scalar curvature, g_{ik}\, is the metric tensor, \Lambda\, is the cosmological constant, G\, is Newton's gravitational constant, c\, is the speed of light in vacuum, and T_{ik}\, is the stress–energy tensor.

Coulomb's law, from the discipline of electromagnetism, describes the electric field between two electric charges (q1 and q2) separated by distance r (with ε0 representing the vacuum permittivity of free space):[128]
 F = \frac{\left|q_1q_2\right|}{4 \pi \varepsilon_0 r^2}.
The fact that π is approximately equal to 3 plays a role in the relatively long lifetime of orthopositronium. The inverse lifetime to lowest order in the fine structure constant \alpha is[129]
\frac{1}{\tau} = 2\frac{\pi^2 - 9}{9\pi}m\alpha^{6},
where m is the mass of the electron.

π is present in some structural engineering formulae, such as the buckling formula derived by Euler, which gives the maximum axial load F that a long, slender column of length L, modulus of elasticity E, and area moment of inertia I can carry without buckling:[130]
F =\frac{\pi^2EI}{L^2}.
The field of fluid dynamics contains π in Stokes' law, which approximates the frictional force F exerted on small, spherical objects of radius R, moving with velocity v in a fluid with dynamic viscosity η:[131]
F =6 \, \pi \, \eta \, R \, v .
The Fourier transform, defined below, is a mathematical operation that expresses time as a function of frequency, known as its frequency spectrum. It has many applications in physics and engineering, particularly in signal processing.[132]
\hat{f}(\xi) = \int_{-\infty}^{\infty} f(x)\ e^{- 2\pi i x \xi}\,dx
Under ideal conditions (uniform gentle slope on an homogeneously erodible substrate), the sinuosity of a meandering river approaches π. The sinuosity is the ratio between the actual length and the straight-line distance from source to mouth. Faster currents along the outside edges of a river's bends cause more erosion than along the inside edges, thus pushing the bends even farther out, and increasing the overall loopiness of the river. However, that loopiness eventually causes the river to double back on itself in places and "short-circuit", creating an ox-bow lake in the process. 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.[133][134]

Memorizing digits

Many persons have memorized large numbers of digits of π, a practice called piphilology.[135] One common technique is to memorize a story or poem in which the word lengths represent the digits of π: The first word has three letters, the second word has one, the third has four, the fourth has one, the fifth has five, and so on. An early example of a memorization aid, originally devised by English scientist James Jeans, is "How I want a drink, alcoholic of course, after the heavy lectures involving quantum mechanics."[135] When a poem is used, it is sometimes referred to as a piem. Poems for memorizing π have been composed in several languages in addition to English.[135]
The record for memorizing digits of π, certified by Guinness World Records, is 67,890 digits, recited in China by Lu Chao in 24 hours and 4 minutes on 20 November 2005.[136][137] In 2006, Akira Haraguchi, a retired Japanese engineer, claimed to have recited 100,000 decimal places, but the claim was not verified by Guinness World Records.[138] Record-setting π memorizers typically do not rely on poems, but instead use methods such as remembering number patterns and the method of loci.[139]

A few authors have used the digits of π to establish a new form of constrained writing, where the word lengths are required to represent the digits of π. The Cadaeic Cadenza contains the first 3835 digits of π in this manner,[140] and the full-length book Not a Wake contains 10,000 words, each representing one digit of π.[141]

In popular culture

Pi Pie at Delft University
A pi pie. The circular shape of pie makes it a frequent subject of pi puns.

Perhaps because of the simplicity of its definition and its ubiquitous presence in formulae, π has been represented in popular culture more than other mathematical constructs.[142]

In the Palais de la Découverte (a science museum in Paris) there is a circular room known as the pi room. On its wall are inscribed 707 digits of π. The digits are large wooden characters attached to the dome-like ceiling. The digits were based on an 1853 calculation by English mathematician William Shanks, which included an error beginning at the 528th digit. The error was detected in 1946 and corrected in 1949.[143]

In Carl Sagan's novel Contact it is suggested that the creator of the universe buried a message deep within the digits of π.[144] The digits of π have also been incorporated into the lyrics of the song "Pi" from the album Aerial by Kate Bush,[145] and a song by Hard 'n Phirm.[146]

Many schools in the United States observe Pi Day on 14 March (March is the third month, hence the date is 3/14).[147] π and its digital representation are often used by self-described "math geeks" for inside jokes among mathematically and technologically minded groups. Several college cheers at the Massachusetts Institute of Technology include "3.14159".[148]

During the 2011 auction for Nortel's portfolio of valuable technology patents, Google made a series of unusually specific bids based on mathematical and scientific constants, including π.[149]

In 1958 Albert Eagle proposed replacing π by τ = π/2 to simplify formulas.[150] However, no other authors are known to use tau in this way. Some people use a different value for tau, τ = 6.283185... = 2π,[151] arguing that τ, as the ratio of a circle's circumference to its radius rather than its diameter, is more natural than π and simplifies many formulas.[152][153] Celebrations of this number, because it approximately equals 6.28, by making 28 June "Tau Day" and eating "twice the pie",[154] have been reported in the media. However this use of τ has not made its way into mainstream mathematics.[155]

In 1897, an amateur mathematician attempted to persuade the Indiana legislature to pass the Indiana Pi Bill, which described a method to square the circle and contained text that implied various incorrect values for π, including 3.2. The bill is notorious as an attempt to establish scientific truth by legislative fiat. The bill was passed by the Indiana House of Representatives, but rejected by the Senate.[156]

Smart city

From Wikipedia, the free encyclopedia https://en.wikipedia.org/wiki/Smart_city Possible scenario of smar...