Brahmagupta's interpolation formula

Brahmagupata'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 Dhyana-graha-adhikara an earlier work of Brahmagupta but of uncertain date. However internal evidences suggest that Dhyana-graha-adhikara could be dated earlier than Brahmasphuta-siddhanta a work of Brahmagupta composed in 628 CE. "Hence the invention of the second-order interpolation formula by Brahmagupta should be placed near the beginning of the second quarter of the 7th century CE, if not earlier."^[1] Brahmagupta was the first to invent and use an interpolation formula using second-order differences in the history of mathematics.^[2][3]

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

Preliminaries

Given a set of tabulated values of a function f(x) in the table below, let it be required to compute the value of f(a), x_r < a < x_r+1.

x	x1	x2	...	xr	xr+1	xr+2	...	xn
f(xr)	f1	f2	...	fr	fr+1	fr+2	...	fn

Assuming that the successively tabulated values of x are equally spaced with a common spacing of h, Aryabhata had considered the table of first differences of the table of values of a function. Writing

$D_r=f_{r+1}-f_r$

the following table can be formed:

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

$f(a)=f_r+ t D_r$ where $t=\frac{a-x_r}{h}$.

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

Brahmagupta's description of the scheme

In Brahmagupta's terminology the difference D_r is the gatakhanda, meaning past difference or the difference that was crossed over, the difference D_r+1 is the bhogyakhanda which is the difference yet to come. Vikala is the amount in minutes by which the interval has been covered at the point where we want to interpolate. In the present notations it is a − x_r. The new expression which replaces f_r+1 − f_r is called sphuta-bhogyakhanda. The description of sphuta-bhogyakhanda is contained in the following Sanskrit couplet (Dhyana-Graha-Upadesa-Adhyaya, 17; Khandaka Khadyaka, IX, 8):^[1]

This has been translated using Bhattolpala's (10th century CE) commentary as follows:^[1][4]

Multiply the vikala by the half the difference of the gatakhanda and the bhogyakhanda and divide the product by 900. Add the result to half the sum of the gatakhanda and the bhogyakhanda if their half-sum is less than the bhogyakhanda, subtract if greater. (The result in each case is sphuta-bhogyakhanda the correct tabular difference.)

This formula was originally stated for the computation of the values of the sine function for which the common interval in the underlying base table was 900 minutes or 15 degrees. So the reference to 900 is in fact a reference to the common interval h.

In modern notation

Brahmagupta's method computation of shutabhogyakhanda can be formulated in modern notation as follows:

sphuta-bhogyakhanda $\displaystyle = \frac{D_r + D_{r+1}}{2} \pm t\frac{|D_r-D_{r+1}|}{2}.$

The ± sign is to be taken according to whether 1/2(D_r + D_r+1) is less than or greater than D_r+1, or equivalently, according to whether D_r < D_r+1 or D_r > D_r+1. Brahmagupta's expression can be put in the following form:

sphuta-bhogyakhanda $\displaystyle = \frac{D_r + D_{r+1}}{2} + t\frac{D_{r+1}-D_r}{2}.$

This correction factor yields the following approximate value for f(a):

$\begin{align} f(a) & = f_r + t\times\text{sphuta-bhogyakhanda}\\ & = f_r + t \frac{D_r + D_{r+1}}{2} + t^2\frac{D_{r+1} - D_r}{2}. \end{align}$

This is Stirling's interpolation formula truncated at the second-order differences.^[5][6] It is not known how Brahmagupta arrived at his interpolation formula.^[1] It is also interesting to note that Brahmagupta has given a separate formula for the case where the values of the independent variable are not equally spaced.

See also

• Brahmagupta's identity
• Brahmagupta matrix
• Brahmagupta–Fibonacci identity

References

1. 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–342. doi:10.1109/5.993400.
4. Raju, C K (2007). Cultural foundations of mathematics: the nature of mathematical proof and the transmission of the calculus from India to Europe in the 16th c. CE. Pearson Education India. pp. 138–140. ISBN 9788131708712.
5. Milne-Thomson, Louis Melville (2000). The Calculus of Finite Differences. AMS Chelsea Publishing. pp. 67–68. ISBN 9780821821077.
6. Hildebrand, Francis Begnaud (1987). Introduction to numerical analysis. Courier Dover Publications. pp. 138–139. ISBN 9780486653631.