Viète's formula

From HandWiki
Short description: Infinite product converging to 2/π
Viète's formula, as printed in Viète's Variorum de rebus mathematicis responsorum, liber VIII (1593)

In mathematics, Viète's formula is the following infinite product of nested radicals representing twice the reciprocal of the mathematical constant π: [math]\displaystyle{ \frac2\pi = \frac{\sqrt 2}2 \cdot \frac{\sqrt{2+\sqrt 2}}2 \cdot \frac{\sqrt{2+\sqrt{2+\sqrt 2}}}2 \cdots }[/math] It can also be represented as: [math]\displaystyle{ \frac2\pi = \prod_{n=1}^{\infty} \cos \frac{\pi}{2^{n+1}} }[/math]

The formula is named after François Viète, who published it in 1593.[1] As the first formula of European mathematics to represent an infinite process,[2] it can be given a rigorous meaning as a limit expression,[3] and marks the beginning of mathematical analysis. It has linear convergence, and can be used for calculations of π,[4] but other methods before and since have led to greater accuracy. It has also been used in calculations of the behavior of systems of springs and masses,[5] and as a motivating example for the concept of statistical independence.

The formula can be derived as a telescoping product of either the areas or perimeters of nested polygons converging to a circle. Alternatively, repeated use of the half-angle formula from trigonometry leads to a generalized formula, discovered by Leonhard Euler, that has Viète's formula as a special case. Many similar formulas involving nested roots or infinite products are now known.

Significance

François Viète (1540–1603) was a French lawyer, privy councillor to two French kings, and amateur mathematician. He published this formula in 1593 in his work Variorum de rebus mathematicis responsorum, liber VIII. At this time, methods for approximating π to (in principle) arbitrary accuracy had long been known. Viète's own method can be interpreted as a variation of an idea of Archimedes of approximating the circumference of a circle by the perimeter of a many-sided polygon,[1] used by Archimedes to find the approximation[6]

[math]\displaystyle{ \frac{223}{71} \lt \pi \lt \frac{22}{7}. }[/math]

By publishing his method as a mathematical formula, Viète formulated the first instance of an infinite product known in mathematics,[7][8] and the first example of an explicit formula for the exact value of π.[9][10] As the first representation in European mathematics of a number as the result of an infinite process rather than of a finite calculation,[11] Eli Maor highlights Viète's formula as marking the beginning of mathematical analysis[2] and Jonathan Borwein calls its appearance "the dawn of modern mathematics".[12]

Using his formula, Viète calculated π to an accuracy of nine decimal digits.[4] However, this was not the most accurate approximation to π known at the time, as the Persian mathematician Jamshīd al-Kāshī had calculated π to an accuracy of nine sexagesimal digits and 16 decimal digits in 1424.[12] Not long after Viète published his formula, Ludolph van Ceulen used a method closely related to Viète's to calculate 35 digits of π, which were published only after van Ceulen's death in 1610.[12]

Beyond its mathematical and historical significance, Viète's formula can be used to explain the different speeds of waves of different frequencies in an infinite chain of springs and masses, and the appearance of π in the limiting behavior of these speeds.[5] Additionally, a derivation of this formula as a product of integrals involving the Rademacher system, equal to the integral of products of the same functions, provides a motivating example for the concept of statistical independence.[13]

Interpretation and convergence

Viète's formula may be rewritten and understood as a limit expression[3]

[math]\displaystyle{ \lim_{n \rightarrow \infty} \prod_{i=1}^n \frac{a_i}{2} = \frac2\pi }[/math]

where

[math]\displaystyle{ \begin{align} a_1 &= \sqrt{2} \\ a_n &= \sqrt{2+a_{n-1}}. \end{align} }[/math]

For each choice of [math]\displaystyle{ n }[/math], the expression in the limit is a finite product, and as [math]\displaystyle{ n }[/math] gets arbitrarily large these finite products have values that approach the value of Viète's formula arbitrarily closely. Viète did his work long before the concepts of limits and rigorous proofs of convergence were developed in mathematics; the first proof that this limit exists was not given until the work of Ferdinand Rudio in 1891.[1][14]

Comparison of the convergence of Viète's formula (×) and several historical infinite series for π. Sn is the approximation after taking n terms. Each subsequent subplot magnifies the shaded area horizontally by 10 times.

The rate of convergence of a limit governs the number of terms of the expression needed to achieve a given number of digits of accuracy. In Viète's formula, the numbers of terms and digits are proportional to each other: the product of the first n terms in the limit gives an expression for π that is accurate to approximately 0.6n digits.[4][15] This convergence rate compares very favorably with the Wallis product, a later infinite product formula for π. Although Viète himself used his formula to calculate π only with nine-digit accuracy, an accelerated version of his formula has been used to calculate π to hundreds of thousands of digits.[4]

Related formulas

Viète's formula may be obtained as a special case of a formula for the sinc function that has often been attributed to Leonhard Euler[16], more than a century later:[1]

[math]\displaystyle{ \frac{\sin x}{x} = \cos\frac{x}{2} \cdot \cos\frac{x}{4} \cdot \cos\frac{x}{8} \cdots }[/math]

Substituting x = π/2 in this formula yields:[17]

[math]\displaystyle{ \frac{2}{\pi} = \cos\frac{\pi}{4} \cdot \cos\frac{\pi}{8} \cdot \cos\frac{\pi}{16} \cdots }[/math]

Then, expressing each term of the product on the right as a function of earlier terms using the half-angle formula:

[math]\displaystyle{ \cos\frac{x}{2} = \sqrt\frac{1+\cos x}{2} }[/math]

gives Viète's formula.[9]

It is also possible to derive from Viète's formula a related formula for π that still involves nested square roots of two, but uses only one multiplication:[18]

[math]\displaystyle{ \pi = \lim_{k\to\infty} 2^{k} \underbrace{\sqrt{2-\sqrt{2+\sqrt{2+\sqrt{2+\sqrt{2+\cdots+\sqrt{2}}}}}}}_{k\text{ square roots}}, }[/math]

which can be rewritten compactly as

[math]\displaystyle{ \begin{align} \pi &= \lim_{k\to\infty}2^k\sqrt{2-a_k} \\[5px] a_1&=0 \\ a_k&=\sqrt{2+a_{k-1}}. \end{align} }[/math]

Many formulae for π and other constants such as the golden ratio are now known, similar to Viète's in their use of either nested radicals or infinite products of trigonometric functions.[8][18][19][20][21][22][23][24]

Derivation

A sequence of regular polygons with numbers of sides equal to powers of two, inscribed in a circle. The ratios between areas or perimeters of consecutive polygons in the sequence give the terms of Viète's formula.

Viète obtained his formula by comparing the areas of regular polygons with 2n and 2n + 1 sides inscribed in a circle.[1][2] The first term in the product, 2/2, is the ratio of areas of a square and an octagon, the second term is the ratio of areas of an octagon and a hexadecagon, etc. Thus, the product telescopes to give the ratio of areas of a square (the initial polygon in the sequence) to a circle (the limiting case of a 2n-gon). Alternatively, the terms in the product may be instead interpreted as ratios of perimeters of the same sequence of polygons, starting with the ratio of perimeters of a digon (the diameter of the circle, counted twice) and a square, the ratio of perimeters of a square and an octagon, etc.[25]

Another derivation is possible based on trigonometric identities and Euler's formula. Repeatedly applying the double-angle formula

[math]\displaystyle{ \sin x = 2\sin\frac{x}{2}\cos\frac{x}{2}, }[/math]

leads to a proof by mathematical induction that, for all positive integers n,

[math]\displaystyle{ \sin x = 2^n \sin\frac{x}{2^n}\left(\prod_{i=1}^n \cos\frac{x}{2^i}\right). }[/math]

The term 2n sin x/2n goes to x in the limit as n goes to infinity, from which Euler's formula follows. Viète's formula may be obtained from this formula by the substitution x = π/2.[9][13]

References

  1. 1.0 1.1 1.2 1.3 1.4 A History of π (2nd ed.). Boulder, Colorado: The Golem Press. 1971. pp. 94–95. ISBN 978-0-88029-418-8. 
  2. 2.0 2.1 2.2 Trigonometric Delights. Princeton, New Jersey: Princeton University Press. 2011. pp. 50, 140. ISBN 978-1-4008-4282-7. 
  3. 3.0 3.1 Eymard, Pierre; Lafon, Jean Pierre (2004). "2.1 Viète's infinite product". The Number pi. Providence, Rhode Island: American Mathematical Society. pp. 44–46. ISBN 978-0-8218-3246-2. https://books.google.com/books?id=qZcCSskdtwcC&pg=PA44. 
  4. 4.0 4.1 4.2 4.3 Kreminski, Rick (2008). "π to thousands of digits from Vieta's formula". Mathematics Magazine 81 (3): 201–207. doi:10.1080/0025570X.2008.11953549. 
  5. 5.0 5.1 Cullerne, J. P.; Goekjian, M. C. Dunn (December 2011). "Teaching wave propagation and the emergence of Viète's formula". Physics Education 47 (1): 87–91. doi:10.1088/0031-9120/47/1/87. 
  6. Beckmann 1971, p. 67.
  7. De Smith, Michael J. (2006). Maths for the Mystified: An Exploration of the History of Mathematics and Its Relationship to Modern-day Science and Computing. Leicester: Matador. p. 165. ISBN 978-1905237-81-4. https://books.google.com/books?id=HlOVuklrQPAC&pg=PA165. 
  8. 8.0 8.1 Moreno, Samuel G.; García-Caballero, Esther M. (2013). "On Viète-like formulas". Journal of Approximation Theory 174: 90–112. doi:10.1016/j.jat.2013.06.006. 
  9. 9.0 9.1 9.2 Morrison, Kent E. (1995). "Cosine products, Fourier transforms, and random sums". The American Mathematical Monthly 102 (8): 716–724. doi:10.2307/2974641. 
  10. Oldham, Keith B.; Myland, Jan C.; Spanier, Jerome (2010). An Atlas of Functions: with Equator, the Atlas Function Calculator. New York: Springer. p. 15. doi:10.1007/978-0-387-48807-3. ISBN 978-0-387-48807-3. https://books.google.com/books?id=UrSnNeJW10YC&pg=PA15. 
  11. Very similar infinite trigonometric series for [math]\displaystyle{ \pi }[/math] appeared earlier in Indian mathematics, in the work of Madhava of Sangamagrama (c. 1340 – 1425), but were not known in Europe until much later. See: "7.3.1 Mādhava on the circumference and arcs of the circle". Mathematics in India. Princeton, New Jersey: Princeton University Press. 2009. pp. 221–234. ISBN 978-0-691-12067-6. https://books.google.com/books?id=DHvThPNp9yMC&pg=PA221. 
  12. 12.0 12.1 12.2 Sidoli, Nathan, ed. "The life of Pi: From Archimedes to ENIAC and beyond". Berlin & Heidelberg: Springer. pp. 531–561. doi:10.1007/978-3-642-36736-6_24. ISBN 978-3-642-36735-9. https://www.carma.newcastle.edu.au/jon/pi-2010.pdf. 
  13. 13.0 13.1 "Chapter 1: From Vieta to the notion of statistical independence". Statistical Independence in Probability, Analysis and Number Theory. Carus Mathematical Monographs. 12. New York: John Wiley & Sons for the Mathematical Association of America. 1959. pp. 1–12. 
  14. {{cite journal | journal = Historisch-litterarische Abteilung der Zeitschrift für Mathematik und Physik | jfm = 23.0263.02 | language = de | pages = 139–140 | title = die Conver rrührenden eigentümlichen Produktentwicklung
  15. "A simple geometric method of estimating the error in using Vieta's product for π". International Journal of Mathematical Education in Science and Technology 38 (1): 136–142. 2007. doi:10.1080/00207390601002799. 
  16. "De variis modis circuli quadraturam numeris proxime exprimendi" (in la). Commentarii Academiae Scientiarum Petropolitanae 9: 222–236. 1738. https://scholarlycommons.pacific.edu/euler-works/74/.  Translated into English by Thomas W. Polaski. See final formula. The same formula is also in "Variae observationes circa angulos in progressione geometrica progredientes" (in la). Opuscula Analytica 1: 345–352. 1783. https://scholarlycommons.pacific.edu/euler-works/561/.  Translated into English by Jordan Bell, arXiv:1009.1439. See the formula in numbered paragraph 3.
  17. Wilson, Robin J. (2018). Euler's pioneering equation: the most beautiful theorem in mathematics (First ed.). Oxford, United Kingdom. pp. 57–58. ISBN 9780198794929. https://swab.zlibcdn.com/dtoken/d543f85cc6cb3fd6ea37e161f89dea19/Euler%E2%80%99s%20Pioneering%20Equation%20The%20most%20beautiful%20theorem%20in%20mathematics%20%28Robin%20Wilson%29%20%28z-lib.org%29.pdf. 
  18. 18.0 18.1 Servi, L. D. (2003). "Nested square roots of 2". The American Mathematical Monthly 110 (4): 326–330. doi:10.2307/3647881. 
  19. Nyblom, M. A. (2012). "Some closed-form evaluations of infinite products involving nested radicals". Rocky Mountain Journal of Mathematics 42 (2): 751–758. doi:10.1216/RMJ-2012-42-2-751. 
  20. Levin, Aaron (2006). "A geometric interpretation of an infinite product for the lemniscate constant". The American Mathematical Monthly 113 (6): 510–520. doi:10.2307/27641976. 
  21. Levin, Aaron (2005). "A new class of infinite products generalizing Viète's product formula for π". The Ramanujan Journal 10 (3): 305–324. doi:10.1007/s11139-005-4852-z. 
  22. "Vieta-like products of nested radicals with Fibonacci and Lucas numbers". Fibonacci Quarterly 45 (3): 202–204. 2007. 
  23. Stolarsky, Kenneth B. (1980). "Mapping properties, growth, and uniqueness of Vieta (infinite cosine) products". Pacific Journal of Mathematics 89 (1): 209–227. doi:10.2140/pjm.1980.89.209. https://projecteuclid.org/euclid.pjm/1102779384. 
  24. Allen, Edward J. (1985). "Continued radicals". The Mathematical Gazette 69 (450): 261–263. doi:10.2307/3617569. 
  25. Rummler, Hansklaus (1993). "Squaring the circle with holes". The American Mathematical Monthly 100 (9): 858–860. doi:10.2307/2324662. 

External links