Bhaskara I's sine approximation formula

In mathematics, Bhaskara I's sine approximation formula is a rational expression in one variable for the computation of the approximate values of the trigonometric sines discovered by Bhaskara I (c. 600 – c. 680), a seventh-century Indian mathematician.[1] This formula is given in his treatise titled Mahabhaskariya. It is not known how Bhaskara I arrived at his approximation formula. However, several historians of mathematics have put forward different theories as to the method Bhaskara might have used to arrive at his formula. The formula is elegant, simple and enables one to compute reasonably accurate values of trigonometric sines without using any geometry whatsoever.[2]

The approximation formula

The formula is given in verses 17 – 19, Chapter VII, Mahabhaskariya of Bhaskara I. A translation of the verses is given below:[3]

(The reference "Rsine-differences 225" is an allusion to Aryabhata's sine table.)

In modern mathematical notations, for an angle x in degrees, this formula gives[3]

\sin x^\circ = \frac{4 x (180-x)}{40500 - x(180-x)}

Equivalent forms of the formula

Bhaskara I's sine approximation formula can be expressed using the radian measure of angles as follows.[1]

\sin x = \frac{16x (\pi - x)}{5\pi^2 - 4x (\pi - x)}

For a positive integer n this takes the following form:[4]

\sin \frac{\pi}{n} = \frac{16(n-1)}{5n^2-4n+4}.

Equivalent forms of Bhaskara I's formula have been given by almost all subsequent astronomers and mathematicians of India. For example, Brahmagupta's (598 – 668 CE) Brhma-Sphuta-Siddhanta (verses 23 – 24, Chapter XIV)[3] gives the formula in the following form:

R \sin x^\circ = \frac{R x(180-x)}{10125 - \frac{1}{4}x(180-x)}

Also, Bhaskara II (1114 – 1185 CE) has given this formula in his Lilavati (Kshetra-vyavahara, Soka No.48) in the following form:

2R\sin x^\circ = \frac{4\times 2R \times 2Rx\times (360R - 2Rx)}{\frac{1}{4}\times 5 \times (360R)^2 - 2Rx\times (360R -2Rx)}

Accuracy of the formula

Figure illustrates the level of accuracy of the Bhaskara I's sine approximation formula. The shifted curves 4 x ( 180 - x ) / ( 40500 - x ( 180 - x ) - 0.2 and sin ( x ) + 0.2 look like exact copies of the curve sin ( x ).

The formula is applicable for values of x° in the range from 0 to 180. The formula is remarkably accurate in this range. The graphs of sin ( x ) and the approximation formula are indistinguishable and are nearly identical. One of the accompanying figures gives the graph of the error function, namely the function,

\sin x^\circ - \frac{4 x (180-x)}{40500 - x(180-x)}

in using the formula. It shows that the maximum absolute error in using the formula is around 0.0016. From a plot of the percentage value of the absolute error, it is clear that the maximum percentage error is less than 1.8. The approximation formula thus gives sufficiently accurate values of sines for all practical purposes. However it was not sufficient for the more accurate computational requirements of astronomy. The search for more accurate formulas by Indian astronomers eventually led to the discovery the power series expansions of sin x and cos x by Madhava of Sangamagrama (c. 1350 – c. 1425), the founder of the Kerala school of astronomy and mathematics.

Graph of the error in Bhaskara I's sine approximation formula
Graph of the percentage error in Bhaskara I's sine approximation formula

Derivation of the formula

Bhaskara I had not indicated any method by which he arrived at his formula. Historians have speculated on various possibilities. No definitive answers have as yet been obtained. Beyond its historical importance of being a prime example of the mathematical achievements of ancient Indian astronomers, the formula is of significance from a modern perspective also. Mathematicians have attempted to derive the rule using modern concepts and tools. Around half a dozen methods have been suggested, each based on a separate set of premises.[2][3] Most of these derivations use only elementary concepts.

Derivation based on elementary geometry [2][3]

Let the circumference of a circle be measured in degrees and let the radius R of the circle be also measured in degrees. Choosing a fixed diameter AB and an arbitrary point P on the circle and dropping the perpendicular PM to AB, we can compute the area of the triangle APB in two ways. Equating the two expressions for the area one gets (1/2) AB × PM = (1/2) AP × BP. This gives

\frac{1}{PM} = \frac{AB}{AP \times BP}.

Letting x be the length of the arc AP, the length of the arc BP is 180 - x. These arcs are much bigger than the respective chords. Hence one gets

\frac{1}{PM} < \frac{2R}{x(180-x)}.

One now seeks two constants α and β such that

\frac{1}{PM} = \alpha \frac{2R}{x(180-x)} + \beta

It is indeed not possible to obtain such constants. However one may choose values for α and β so that the above expression is valid for two chosen values of the arc length x. Choosing 30° and 90° as these values and solving the resulting equations, one immediately gets Bhaskara I's sine approximation formula.

Derivation starting with a general rational expression

Assuming that x is in radians, one may seek an approximation to sin (x) in the following form:

\sin x = \frac{a+bx+cx^2}{p+qx+rx^2}

The constants a, b, c, p, q and r (only five of them are independent) can be determined by assuming that the formula must be exactly valid when x = 0, π/6, π/2, π, and further assuming that it has to satisfy the property that sin (x) = sin (π - x).[2][3] This procedure produces the formula expressed using radian measure of angles.

An elementary argument[4]

Comparison of graphs of the parabolas
x(180 − x)/8100 and x(180 − x)/9000
with the graph of sin(x) (x in degrees).

The part of the graph of sin(x) in the range from 0° to 180° "looks like" part of a parabola through the points (0, 0) and (180, 0). The general such parabola is

k x ( 180 - x ).\,

The parabola that also passes through (90, 1) (which is the point corresponding to the value sin(90°) = 1) is

\frac{x ( 180 - x )}{ 90 \times 90} = \frac{x ( 180 - x )}{ 8100}.

The parabola which also passes through (30, 1/2) (which is the point corresponding to the value sin(30°) = 1/2) is

\frac{x ( 180 - x )}{2 \times 30 \times 150} = \frac{x(180-x)}{9000}.

These expressions suggest a varying denominator which takes the value 90 × 90 when x = 90 and the value 2 × 30 × 150 when x = 30. That this expression should also be symmetrical about the line ' x = 90' rules out the possibility of choosing a linear expression in x. Computations involving x(180 − x) might immediately suggest that the expression could be of the form

8100a + bx ( 180 - x ).\,

A little experimentation (or by setting up and solving two linear equations in a and b) will yield the values a = 5/4, b = −1/4. These give Bhaskara I's sine approximation formula.

See also

References

  1. ↑ 1.0 1.1 J J O'Connor and E F Robertson (November 2000). "Bhaskara I". School of Mathematics and Statistics University of St Andrews, Scotland. Archived from the original on 23 March 2010. Retrieved 22 April 2010.
  2. ↑ 2.0 2.1 2.2 2.3 Glen Van Brummelen (2009). The mathematics of the heavens and the earth: the early history of trigonometry. Princeton University Press. ISBN 978-0-691-12973-0. (p.104)
  3. ↑ 3.0 3.1 3.2 3.3 3.4 3.5 R.C. Gupta (1967). "Bhaskara I' approximation to sine". Indian Journal of HIstory of Science 2 (2). Retrieved 20 April 2010.
  4. ↑ 4.0 4.1 George Gheverghese Joseph (2009). A passage to infinity : Medieval Indian mathematics from Kerala and its impact. New Delhi: SAGE Publications India Pvt. Ltd. ISBN 978-81-321-0168-0. (p.60)

Further references

  1. R.C..Gupta, On derivation of Bhaskara I's formula for the sine, Ganita Bharati 8 (1-4) (1986), 39-41.
  2. T. Hayashi, A note on Bhaskara I's rational approximation to sine, Historia Sci. No. 42 (1991), 45-48.