Madhava series

In mathematics, a Madhava series is one of the three Taylor series expansions for the sine, cosine, and arctangent functions discovered in 14th or 15th century Kerala by the mathematician and astronomer Madhava of Sangamagrama (c. 1350 – c. 1425) or his followers in the Kerala school of astronomy and mathematics.[1] Using modern notation, these series are:

All three series were later independently discovered in 17th century Europe. The series for sine and cosine were rediscovered by Isaac Newton in 1669,[2] and the series for arctangent was rediscovered by James Gregory in 1671 and Gottfried Leibniz in 1673,[3] and is conventionally called Gregory's series. The specific value can be used to calculate the circle constant π, and the arctangent series for 1 is conventionally called Leibniz's series.

In recognition of Madhava's priority, in recent literature these series are sometimes called the Madhava–Newton series,[4] Madhava–Gregory series,[5] or Madhava–Leibniz series[6] (among other combinations).[7]

No surviving works of Madhava contain explicit statements regarding the expressions which are now referred to as Madhava series. However, in the writing of later Kerala school mathematicians Nilakantha Somayaji and Jyeshthadeva one can find unambiguous attributions of these series to Madhava. These later works also include proofs and commentary which suggest how Madhava may have arrived at the series.

Madhava series in "Madhava's own words"

None of Madhava's works, containing any of the series expressions attributed to him, have survived. These series expressions are found in the writings of the followers of Madhava in the Kerala school. At many places these authors have clearly stated that these are "as told by Madhava". Thus the enunciations of the various series found in Tantrasamgraha and its commentaries can be safely assumed to be in "Madhava's own words". The translations of the relevant verses as given in the Yuktidipika commentary of Tantrasamgraha (also known as Tantrasamgraha-vyakhya) by Sankara Variar (circa. 1500 - 1560 CE) are reproduced below. These are then rendered in current mathematical notations.[8][9]

Madhava's sine series

In Madhava's own words

Madhava's sine series is stated in verses 2.440 and 2.441 in Yukti-dipika commentary (Tantrasamgraha-vyakhya) by Sankara Variar. A translation of the verses follows.

Multiply the arc by the square of the arc, and take the result of repeating that (any number of times). Divide by the squares of the successive even numbers (such that current is multiplied by previous) increased by that number and multiplied by the square of the radius. Place the arc and the successive results so obtained one below the other, and subtract each from the one above. These together give the jiva [sine], as collected together in the verse beginning with "vidvan" etc.

Rendering in modern notations

Let r denote the radius of the circle and s the arc-length.

  • The following numerators are formed first:
  • These are then divided by quantities specified in the verse.
  • Place the arc and the successive results so obtained one below the other, and subtract each from the one above to get jiva:

Transformation to current notation

Let θ be the angle subtended by the arc s at the centre of the circle. Then s = r θ and jiva = r sin θ. Substituting these in the last expression and simplifying we get

which is the infinite power series expansion of the sine function.

Madhava's reformulation for numerical computation

The last line in the verse as collected together in the verse beginning with "vidvan" etc. is a reference to a reformulation of the series introduced by Madhava himself to make it convenient for easy computations for specified values of the arc and the radius. For such a reformulation, Madhava considers a circle one quarter of which measures 5400 minutes (say C minutes) and develops a scheme for the easy computations of the jivas of the various arcs of such a circle. Let R be the radius of a circle one quarter of which measures C. Madhava had already computed the value of π using his series formula for π.[10] Using this value of π, namely 3.1415926535922, the radius R is computed as follows: Then

R = 2 × 5400 / π = 3437.74677078493925 = 3437 arcminutes 44 arcseconds 48 sixtieths of an arcsecond = 3437 44 48.

Madhava's expression for jiva corresponding to any arc s of a circle of radius R is equivalent to the following:

Madhava now computes the following values:

No. Expression Value Value in Katapayadi system
   1       R × (π / 2)3 / 3!       2220   39   40       ni-rvi-ddhā-nga-na-rē-ndra-rung   
   2       R × (π / 2)5 / 5!       273   57   47       sa-rvā-rtha-śī-la-sthi-ro   
   3       R × (π / 2)7 / 7!       16   05   41       ka-vī-śa-ni-ca-ya   
   4       R × (π / 2)9 / 9!       33   06       tu-nna-ba-la   
   5       R × (π / 2)11 / 11!       44       vi-dvān   

The jiva can now be computed using the following scheme:

jiva = s (s / C)3 [ (2220 39 40) (s / C)2 [ (273 57 47) (s / C)2 [ (16 05 41) (s / C)2[ (33 06) (s / C)2 (44 ) ] ] ] ].

This gives an approximation of jiva by its Taylor polynomial of the 11'th order. It involves one division, six multiplications and five subtractions only. Madhava prescribes this numerically efficient computational scheme in the following words (translation of verse 2.437 in Yukti-dipika):

vi-dvān, tu-nna-ba-la, ka-vī-śa-ni-ca-ya, sa-rvā-rtha-śī-la-sthi-ro, ni-rvi-ddhā-nga-na-rē-ndra-rung . Successively multiply these five numbers in order by the square of the arc divided by the quarter of the circumference (5400), and subtract from the next number. (Continue this process with the result so obtained and the next number.) Multiply the final result by the cube of the arc divided by quarter of the circumference and subtract from the arc.

Madhava's cosine series

In Madhava's own words

Madhava's cosine series is stated in verses 2.442 and 2.443 in Yukti-dipika commentary (Tantrasamgraha-vyakhya) by Sankara Variar. A translation of the verses follows.

Multiply the square of the arc by the unit (i.e. the radius) and take the result of repeating that (any number of times). Divide (each of the above numerators) by the square of the successive even numbers decreased by that number and multiplied by the square of the radius. But the first term is (now)(the one which is) divided by twice the radius. Place the successive results so obtained one below the other and subtract each from the one above. These together give the śara as collected together in the verse beginning with stena, stri, etc.

Rendering in modern notations

Let r denote the radius of the circle and s the arc-length.

  • The following numerators are formed first:
  • These are then divided by quantities specified in the verse.
  • Place the arc and the successive results so obtained one below the other, and subtract each from the one above to get śara:

Transformation to current notation

Let θ be the angle subtended by the arc s at the centre of the circle. Then s = and śara = r(1 − cos θ). Substituting these in the last expression and simplifying we get

which gives the infinite power series expansion of the cosine function.

Madhava's reformulation for numerical computation

The last line in the verse as collected together in the verse beginning with stena, stri, etc. is a reference to a reformulation introduced by Madhava himself to make the series convenient for easy computations for specified values of the arc and the radius. As in the case of the sine series, Madhava considers a circle one quarter of which measures 5400 minutes (say C minutes) and develops a scheme for the easy computations of the śaras of the various arcs of such a circle. Let R be the radius of a circle one quarter of which measures C. Then, as in the case of the sine series, Madhava gets R = 3437 44 48.

Madhava's expression for śara corresponding to any arc s of a circle of radius R is equivalent to the following:

Madhava now computes the following values:

No. Expression Value Value in Katapayadi system
   1       R × (π / 2)2 / 2!       4241   09   00       u-na-dha-na-krt-bhu-re-va   
   2       R × (π / 2)4 / 4!       872   03   05       mī-nā-ngo-na-ra-sim-ha   
   3       R × (π / 2)6 / 6!       071   43   24       bha-drā-nga-bha-vyā-sa-na   
   4       R × (π / 2)8 / 8!       03   09   37       su-ga-ndhi-na-ga-nud   
   5       R × (π / 2)10 / 10!       05   12       strī-pi-śu-na   
   6       R × (π / 2)12 / 12!       06       ste-na   

The śara can now be computed using the following scheme:

śara = (s / C)2 [ (4241 09 00) (s / C)2 [ (872 03 05 ) (s / C)2 [ (071 43 24) (s / C)2[ (03 09 37) (s / C)2 [(05 12) − (s / C)2 (06) ] ] ] ] ]

This gives an approximation of śara by its Taylor polynomial of the 12'th order. This also involves one division, six multiplications and five subtractions only. Madhava prescribes this numerically efficient computational scheme in the following words (translation of verse 2.438 in Yukti-dipika):

The six stena, strīpiśuna, sugandhinaganud, bhadrāngabhavyāsana, mīnāngonarasimha, unadhanakrtbhureva. Multiply by the square of the arc divided by the quarter of the circumference and subtract from the next number. (Continue with the result and the next number.) Final result will be utkrama-jya (R versed sign).

Madhava's arctangent series

In Madhava's own words

Madhava's arctangent series is stated in verses 2.206 2.209 in Yukti-dipika commentary (Tantrasamgraha-vyakhya) by Sankara Variar. A translation of the verses is given below.[11] Jyesthadeva has also given a description of this series in Yuktibhasa.[12][13][14]

Now, by just the same argument, the determination of the arc of a desired sine can be (made). That is as follows: The first result is the product of the desired sine and the radius divided by the cosine of the arc. When one has made the square of the sine the multiplier and the square of the cosine the divisor, now a group of results is to be determined from the (previous) results beginning from the first. When these are divided in order by the odd numbers 1, 3, and so forth, and when one has subtracted the sum of the even(-numbered) results from the sum of the odd (ones), that should be the arc. Here the smaller of the sine and cosine is required to be considered as the desired (sine). Otherwise, there would be no termination of results even if repeatedly (computed).

By means of the same argument, the circumference can be computed in another way too. That is as (follows): The first result should by the square root of the square of the diameter multiplied by twelve. From then on, the result should be divided by three (in) each successive (case). When these are divided in order by the odd numbers, beginning with 1, and when one has subtracted the (even) results from the sum of the odd, (that) should be the circumference.

Rendering in modern notations

Let s be the arc of the desired sine (jya or jiva) y. Let r be the radius and x be the cosine (kotijya).

  • The first result is .
  • Form the multiplier and divisor .
  • Form the group of results:
  • These are divided in order by the numbers 1, 3, and so forth:
  • Sum of odd-numbered results:
  • Sum of even-numbered results:
  • The arc is now given by

Transformation to current notation

Let θ be the angle subtended by the arc s at the centre of the circle. Then s = rθ, x = kotijya = r cos θ and y = jya = r sin θ. Then y / x = tan θ. Substituting these in the last expression and simplifying we get

  • .

Letting tan θ = q we finally have

Another formula for the circumference of a circle

The second part of the quoted text specifies another formula for the computation of the circumference c of a circle having diameter d. This is as follows.

Since c = π d this can be reformulated as a formula to compute π as follows.

This is obtained by substituting q = (therefore θ = π / 6) in the power series expansion for tan−1 q above.

Comparison of convergence of various infinite series for π


Comparison of the convergence of two Madhava series (the one with 12 in dark blue) 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. (click for detail)

See also

Notes

  1. Gupta 1987; Katz 1995; Roy 2021, Ch. 1. Power Series in Fifteenth-Century Kerala, pp. 1–22
  2. Newton (1669) De analysi per aequationes numero terminorum infinitas was circulated as a manuscript but not published until 1711. For context, see:
    Roy 2021, Ch. 8. De Analysi per Aequationes Infinitas, pp. 165–185.
    Leibniz later included the series for sine and cosine in Leibniz (1676) De quadratura arithmetica circuli ellipseos et hyperbola cujus corollarium est trigonometria sine tabulis, which was only finally published in 1993. However, he had been sent Newton's sine and cosine series by Henry Oldenburg in 1675 and did not claim to have discovered them. See:
    Probst, Siegmund (2015). "Leibniz as reader and second inventor: The cases of Barrow and Mengoli". In Goethe, N.; Beeley, P.; Rabouin, D. (eds.). G.W. Leibniz, Interrelations between Mathematics and Philosophy. Springer. pp. 111–134. doi:10.1007/978-94-017-9664-4_6.
  3. Gregory received a letter from John Collins including Newton's sine and cosine series in late 1670. He discovered the general Taylor series and sent a now-famous letter back to Collins in 1671 including several specific series including the arctangent. See Roy 1990.
    Horvath, Miklos (1983). "On the Leibnizian quadrature of the circle" (PDF). Annales Universitatis Scientiarum Budapestiensis (Sectio Computatorica). 4: 75–83.
  4. For example:
    Plofker, Kim (2005). "Relations between approximations to the sine in Kerala mathematics". In Emch, Gérard G.; Sridharan, R.; Srinivas, M. D. (eds.). Contributions to the History of Indian Mathematics. Gurgaon: Hindustan Book Agency. pp. 135–152. doi:10.1007/978-93-86279-25-5_6.
    Filali, Mahmoud (2012). "Harmonic analysis and applications". Kybernetes. 41: 129–144. doi:10.1108/03684921211213160.
  5. For example: Gupta 1973; Joseph 2011, p. 428;
    Levrie, Paul (2011). "Lost and Found: An Unpublished ζ(2)-Proof". Mathematical Intelligencer. 33: 29–32. doi:10.1007/s00283-010-9179-y.
  6. For example: Gupta 1992;
    Pouvreau, David (2015). "Sur l'accélération de la convergence de la série de Madhava-Leibniz". Quadrature (in French). 97: 17–25.
    Young, Paul Thomas (2022). "From Madhava–Leibniz to Lehmer's Limit". American Mathematical Monthly. 129 (6): 524–538. doi:10.1080/00029890.2022.2051405.
  7. For example,
    Madhava–Gregory–Leibniz series: Benko, David; Molokach, John (2013). "The Basel Problem as a Rearrangement of Series". College Mathematics Journal. 44 (3): 171–176. doi:10.4169/college.math.j.44.3.171.
    Madhava–Leibniz–Gregory series: Danesi, Marcel (2021). "1. Discovery of π and Its Manifestations". Pi (π) in Nature, Art, and Culture. Brill. pp. 1–30. doi:10.1163/9789004433397_002.
    Nilakantha–Gregory series: Campbell, Paul J. (2004). "Borwein, Jonathan, and David Bailey, Mathematics by Experiment". Reviews. Mathematics Magazine. 77 (2): 163. doi:10.1080/0025570X.2004.11953245.
    Gregory–Leibniz–Nilakantha formula: Gawrońska, Natalia; Słota, Damian; Wituła, Roman; Zielonka, Adam (2013). "Some generalizations of Gregory's power series and their applications" (PDF). Journal of Applied Mathematics and Computational Mechanics. 12 (3).
  8. Bag 1976.
  9. Raju 2007, pp. 114–120.
  10. Raju 2007, p. 119.
  11. Raju 2007, p. 231.
  12. O'Connor, John J.; Robertson, Edmund F. (2000). "Madhava of Sangamagramma". MacTutor History of Mathematics archive.
  13. Gupta 1973.
  14. Sarma 1972.

References

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.