Brahmagupta's interpolation formula

Brahmagupta's interpolation formula is a second-order polynomial interpolation formula developed by the Indian mathematician and astronomer Brahmagupta (598–668 CE) in the early 7th century CE. The Sanskrit couplet describing the formula can be found in the supplementary part of Khandakadyaka a work of Brahmagupta completed in 665 CE.[1] The same couplet appears in Brahmagupta's earlier Dhyana-graha-adhikara, which was probably written "near the beginning of the second quarter of the 7th century CE, if not earlier."[1] Brahmagupta was one of the first to describe and use an interpolation formula using second-order differences.[2][3]

Brahmagupta's interpolation formula is equivalent to modern-day second-order Newton–Stirling interpolation formula.


  x    x2    ...     xr    xr+1    ...     xn 
  Differences   D1   ...   Dr  Dr+1 ...   Dn

Mathematicians prior to Brahmagupta used a simple linear interpolation formula. The linear interpolation formula to compute f(a) is

where .

For the computation of f(a), Brahmagupta replaces Dr with another expression which gives more accurate values and which amounts to using a second-order interpolation formula.

  1. ^ a b Gupta, R. C. "Second-order interpolation in Indian mathematics upto the fifteenth century". Indian Journal of History of Science. 4 (1 & 2): 86–98.
  2. ^ Van Brummelen, Glen (2009). The mathematics of the heavens and the earth: the early history of trigonometry. Princeton University Press. p. 329. ISBN 9780691129730. (p.111)
  3. ^ Meijering, Erik (March 2002). "A Chronology of Interpolation From Ancient Astronomy to Modern Signal and Image Processing". Proceedings of the IEEE. 90 (3): 319–321. doi:10.1109/5.993400.

© MMXXIII Rich X Search. We shall prevail. All rights reserved. Rich X Search