# Use of Calculus in Hindu Mathematics

## Kripa Shankar Shukla

### Argara, Hussaingunj Crossing Lucknow 226 001

#### 1. Differential Calculus

A Controversy
Attention was first drawn to the occurrence of the differential formula
δ (sin θ) = cos θ δ θ
in Bhaskara II's (1150) Siddhanta Siromani by Pandit BapuDeva Sastri in 1858. The Pandit published a summarised translation of the passages which involve the use of the above formula. His summary was defective in so far as it did not bring into prominence the idea of the infinitesimal increment which underlies Bhaskara's analysis. Without making clear to his readers, the full significance of Bhaskara's result, the Pandit made the mistake of asserting — what was plain to him — that Bhaskara was fully acquainted with the principles of the differential calculus.

The Pandit was adversely criticised by Spotiswoode, who without consulting the original on which the Pandit based his conclusions, remarked (1) that Bapu Deva Sastri had overstated his case in saying that Bhaskaraciirya was fully acquainted with the principles of the differential calculus, (2) that there was no allusion to the most essential feature of the differential calculus, viz. the infinitesimal magnitudes of the intervals of time and space therein employed, and (3) that the approximative character of the result was not realized.

Since the above controversy took place no serious investigation of the subject seems to have been made by any scholar. In order that the reader may be better able to judge the merit of the Hindu claim to the invention of the differential calculus, it is desirable that the problems which required the use of the above differential formula be stated first.

Problems in Astronomy
The calculation of eclipses is one of the most important problems of Astronomy. In ancient days this problem was probably more important than it is now, because the exact time and duration of the eclipses could not be foretold on account of lack of the necessary mathematical equipment on the part of the astronomer. In India, the Hindus observed fast and performed various other religious rites on the occasion of eclipses. Thus their calculation was a matter of national importance. It afforded the Hindu astronomer a means of demonstrating the accuracy of his science and his own ability to the public who patronised him. The problem of the calculation of conjunction of planets and occupation of stars was equally important both from scientific as well as religious view points.

In problems of the above nature it is essential to determine the true instantaneous motion of a planet or star at any particular instant. This instantaneous motion was called by the Hindu astronomers tat-kalika-gati. The formula giving the tat-kalika-gati (instantaneous motion) is given by Aryabhata and Brahmagupta in the following form:
u'- v' = v' - v ± e (sin w' - sin w)                   (i)
where u, v, w are the true longitude, mean longitude, mean anomaly respectively at any particular time and u', v', w' the values of the respective quantities at a subsequent instant; and e is the eccentricity or the sine of the greatest equation of the orbit. The tat-kalika-gati is the difference u'-u between the true longitudes at the two positions under consideration. Aryabhata and Brahmagupta used the sine table to find the value of (sin w' — sin w). The sine table used by them was tabulated at intervals of 3° 45' and thus was entirely unsuited for the purpose. To get the values of sines of angles, not occurring in the table, recourse was taken to interpolation formulae, which were incorrect because the law of variation of the difference was not known.

A Differential Formula
Manjula (932) was the first Hindu astronomer to state that the difference of the sines,             sin w' - sin w = (w' - w) cos w,
where (w' - w) is small.

He says:
"True motion in minutes is equal to the cosine (of the mean anomaly) multiplied by the difference (of the mean anomalies) and divided by the cheda, added or subtracted contrarily (to the mean motion)."
Thus according to Manjula formula (i) becomes
u' - u = v' - v ± e(w' - w) cos w,                  (ii)
which, in the language of the differential calculus, may be written as
δu = δv ± e cos θ δ θ.

We cannot say exactly what was the method employed by Manjula to obtain formula (ii). The formula occurs also in the works of Aryabhata II (950). Bhaskara II (1150), and later writers. Bhaskara II indicates the method of obtaining the differential of sine θ. His method is probably the same as that employed by his predecessors.

Proof of the Differential Formula
Let a point P (See Fig. 1) move on a circle. Let its position at two successive intervals be denoted by P and Q. Now, if P and Q are taken very near each other, the direction of motion in the interval PQ is the same as that of the tangent at P. Let PT be measured along the tangent at P equal to the arc PQ. Then PT would be the motion of the point P if its velocity at P had not changed direction.

Discussing the motion of planets, Bhaskaracarya says: "The difference between the longitudes of a planet found at any time on a certain day and at the same time on the following day is called its (sphuta) gati (true rate of motion) for that interval of time."

"This is indeed rough motion (sthulagati). I now describe the fine (suksma) instantaneous (tat-kalika) motion. The tatkalika-gati (instantaneous motion) of a planet is the motion which it would have, had its velocity during any given interval of time remained uniform."

During the course of the above statement, Bhaskara II observes that the tat-kalika-gati is suksma ("fine" as opposed lo rough), and for that the interval must be taken to be very small, so that the motion would be very small. This small interval of time has been said to be equivalenttoa ksana which according to the Hindus is an infinitesimal interval of time (immeasurably small). It will be apparent from the above that Bhaskara did really employ the notion of the infinitesimal in his definition of Tat-kalika-gati.

But in actual practice, the intervals that are considered are not infinitesimal. How are we, then, to apply the notion of tat-kalika-gati to actual problems ? The answer to the above question is given by Bhaskara II as follows:

In equation (i) we have to consider the sine-difference (sin w' - sin w). Let an arc of 90° be divided into n parts each equal to A, and let us consider the sine differences R(sin A - sin O), R (sin 2A - sin A), R (sin 3A - sin 2A)... These differences are termed Bhogya khanda. Bhaskara-II says: "These are not equal to each other but gradually decrease, and consequently while the increase of the arc is uniform, the increment of the sine varies — on account of deflection of the arc."

In the figure given above, let the arc PQ = A. Then
R (sin BOQ - sin BOP) = QN - PM = Qn
which is the Bhogya Khanda. Bhaskara introduces the notion of Tat-kalika Bhogya Khanda (instantaneous sine difference) in order to find the variation of the sine at P. According to him if the arc BP instead of being deflected towards Q, be increased in the direction of the tangent, so that PT = PQ = A, then TS - PM = Tr is the Tat-kalika Bhogya Khanda of the sine PT, i.e. the "instantaneous sine difference". By having recourse to this artifice Bhaskara II avoids the use of the infinitesimal in his analysis. It should be borne in mind that the "instantaneous sine difference" for a finite arc PQ, is a purely artificial quantity created with a special end in view, and is different from the actual "sine difference" R (sin BOQ - sin BOP).

Now from the similar triangles PTr and PMO, we at once derive the proportion
R : PT : : R cos w : Tr ................             (iii)
Tr = PT cos w.
But     Tr = R(sin w' - sin w) and PT = R(w' - w)
(sin w' - sin w) = (w' - w) cos w.
Thus the Tat-kalika Bhogya Khanda (the instantaneous sine difference) in modern notation is
δ (sin θ) = cos θ δ θ.
This formula has been used by Bhaskara to calculate the ayana-valana ("angle of position").

If the above were the only result occurring in Bhaskara II's work, one would be justified in not accepting the conclusions of Pandit Bapu Deva Sastri. There is however other evidence in Bhaskara II's work to show that he did actually know the principles of the differential calculus. This evidence consists partly in the occurrence of the two most important results of the differential calculus:
(i) He has shown that when a variable attains the maximum value its differential vanishes.
(ii) He shows that when a planet is either in apogee or in perigee the equation of the centre vanishes, hence he concludes that for some intermediate position the increment of the equation of centre (i.e. the differential) also vanishes.
The second of the above results is the celebrated Rolle's Theorem, the mean value theorem of the differential calculus.

Remarks
The use of a formula involving differentials in the works of ancient Hindu mathematicians has been established beyond the possibility of any doubt. That the motions of instantaneous variation and that of motion entered into the Hindu idea of differentials as found in works of Manjula, Aryabhata II and Bhaskara II is apparent from the epithet Tat-kalika (instantaneous) gati (motion) to denote these differentials. The main contribution of Bhaskara II to the theory of these differentials, which were already worked out by his predecessors, seems to be his proof of the formula by the rule of proportion without actually using the infinitesimal or varying quantities. He has, however, made it quite clear that the differentials give true results only when very small variations are concerned.

Nilakarttha's Result
Nilakantha (c. 1500) in his commentary on the Aryabhatiya has given proofs, on the theory of proportion (similar triangles) of the following results.
(1) The sine-difference sin (θ + δ θ) - sin θ varies as the cosine and decreases as θ increases.
(2) The cosine-difference cos (θ + δ θ) - cos θ varies as the sine negatively and numerically increases as θ increases.
He has obtained the following formulae:
$(a) sin(\theta + \delta \theta)-sin\theta = 2sin \frac{\delta \theta}{2} cos(\theta + \frac{\delta \theta}{2} ). \\ (b) cos(\theta + \delta \theta)-cos\theta = -2sin \frac{\delta \theta}{2} sin(\theta + \frac{\delta \theta}{2}).$

The above results are true for all values of δθ whether big or small. There is nothing new in the above results. They are simply expressions as products of sine and cosine differences.

But what is important in Nilakantha's work is his study of the second differences. These are studied geometrically by the help of the property of the circle and of similar triangles. Denoting by Δ2 (sin θ) and Δ2 (cos θ), the second differences of these functions, Nilakarttha's results may be stated as follows:
(1) The difference of the sine-difference varys as the sine negatively and increases (numerically) with the angle.
(2) The difference of the cosine-difference varys as the cosine negatively and decreases (numerically) with the angle.
For Δ2 (sin θ), Nilakantha has obtained the following formula
$\Delta_2(sin \theta) = -sin \theta \left[sin (\frac{\Delta \theta}{2})\right]^2$
Besides the above, Nilakantha, has made use of a result involving the differential of an inverse sine function. This result, expressed in modern notation, is
$\delta (sin^{-1} e sin w) = \frac{e cos w}{\sqrt{1-e^2 sin^2 w}}\delta w$

In the writings of Acyuta (1550-1621 A.D.) we find use of the differential of a quotient also
$\delta \left[\frac{e sin w}{1 \pm e cos w}\right] = \frac{e cos w \pm (e sin w)^2/(1 \pm e cos w)}{1 \pm e cos w}\delta w$

#### 2. Method of infinitesimal-integration

Surface of the Sphere
For calculating the area of the surface of a sphere Bhaskara II (1150) describes two methods which are almost the same as we usually employ now for the same purpose.

First Method: "Make a spherical ball of clay or of wood. On it take a (vertical) circumference circle and divide this into 21600 parts. Mark a point on the top of it. With that point as the centre and with the radius equal to the 96th part of the circumference, i.e. to 225', describe a circle. Again with same point as the centre with twice that arc as radius describe another circle; with thrice that another circle; and so on up to 24 times. Thus there will be 24 circles in all. The radii of these circles will be the jya 225' (= R sin 225'), etc. From them the lengths of the circles can be determined by proportion. Now the length of the extreme circle is 21600' and its radius is 3438'. Multiplying the Rsines (of 225', 450', etc.) by 21600 and dividing by 3438, we shall obtain the lengths of the circles. Between two and two of these circles there lies annular strips and there are altogether 24 such. They will be many more in case of many R-sines being taken into consideration (bahujya-pakse-vahuni syuh). In each annulus considering the larger circle at the lower end as the base and the smaller circle at the top as the face and 225'as the altitude (of the trapezium), find its area by means of the rule 'half the sum of the base and the face multiplied by the altitude etc. Similarly the areas of all the annular figures severally can be found. The sum of all these areas is equal to the area of the surface of half of the sphere. So twice that is the area of the surface of the whole sphere. And that is equal to the product of the diameter and the circumference."

In other words, if Tn denotes the nth jya (or Rsine), Cn the circumference of the corresponding circle, An the area of the nth annulus and S the area of the surface of the sphere, then we shall have
$C_n = \frac{21600}{3438}T_n \\ A_n = \frac{C_n + C_{n+1}}{2}x225 \\ = \frac{225 x 21600}{2 x 3438}(T_n + T_{n+1})$ Therefore,
$\frac{1}{2}S = \Sigma A_n \\ = \frac{225 x 21600}{2 x 3438}\Sigma(T_n + T_{n+1})$
the summation being taken so as to include all the Rsines in a quadrant of the circle. Since there are ordinarily 24 Rsines in a Hindu trigonometrical table, we have
$\frac{1}{2}S = \frac{225 x 21600}{3438}\Sigma(T_1 + T_2 + ... + T_{23} + \frac{1}{2}T_{24}) \\ = \frac{21600 x 225 x 52513}{3438} \\ = 21600 x 3436.7 \ldots$
Hence approximately
S = 21600 x 2 x 3437
Area of the surface = circumference x diameter.

Second Method: "Suppose the (horizontal) circumference-circle on the surface of the sphere to be divided into parts as many as four times the number of Rsines (in a quadrant). As the surface of an emblic myrobalan is seen divided into vapras (i.e. lunes) by lines passing through its face (or top) and bottom, so the surface of the sphere should be divided into lunes by vertical circles as many as the parts of the above mentioned (horizontal) circumference-circle. Then the area of each lune should be determined by (breaking it up into) parts. And this area of a lune is equal to the sum of all the Rsines diminished by half the radius and divided by the semi-radius. Since that is again equal to the diameter of the sphere, so it has been said that the area of the surface of a sphere is equal to the product of its circumference and diameter."

The method has been further elucidated by him in his gloss thus:

"As many as are Rsines in the table of any particular work selected, take four times that number, and suppose the (horizontal) circumference-circle on the sphere is to be divided into, as many parts. Like the natural lines seen on the surface of a round emblic myrobalan passing through its face and base and thus dividing it into lunes, draw circles on the surface of the given sphere, passing through its top and bottom and thereby dividing it into lunes as many as the number of parts into which the (horizontal) circumference-circle is divided. Next the area of each lune has to be determined. It can be done thus: For instance in the Dhivrddhida, there are 24 Rsines. So suppose the (horizontal) circumference-circle measures 96 cubits. On drawing the vertical circles through every cubit, there will be as many lunes. Then the upper half of any one lune ondrawing the transverse arcs at distances of every cubit, will be divided into portions equal to the number of Rsines, that is, 24. The lengths of these transverse lines will be obtained by dividing the Rsines severally by the radius. Of these the lowest line measures one cubit; but the upper and upper ones are a little smaller and smaller according to the Rsines. But the altitude is all along one cubit in length. Now by finding the area of each portion in accordance with the rule, "half the sum of the top and the base multiplied by the altitude etc: they should be added together. This sum gives the area of half a lune; twice that is the area of a lune. For the determination of that the rule is, "the sum of all the Rsines minus half the radius etc." Now the sum of all the Rsines, 225 etc., is 54233. This diminished by the semi-radius becomes 52514. Dividing the result by the semi-radius we get the area of each lune as 30;33. Now 30;33 is equal to the diameter of a circle whose circumference measures 96. And as the number of lunes is equal to the number of portions of the circumference it is consequently proved that the area of the surface of a sphere is equal to the product of its circumference and diameter".

If ln denotes the length of the nth transverse arc, we have
$l_n = \frac{T_nxl}{R}$
Therefore,
$\textrm{area of a lune } = 2 x \Sigma \frac{1}{2}(l_n + l_{n+1})x1 \\ = 2 \Sigma \frac{1}{2R}(T_n + T_{n+1})$
the summation being taken so as to include all the Rsines. Hence
$\textrm{area of a lune } = 2 x \frac{1}{R}(T_1 + T_2 + \ldots + T_{23}+\frac{1}{2}T_{24}) \\ = \frac{1}{R/2} \left(T_1 + T_2 + \ldots + T_{24}-\frac{R}{2} \right) \\ = 30;32,94 \ldots$
$\textrm{Now } 96 x \frac{1250}{3927} = 30; 33,46 \ldots$
Hence the area of a lune is numerically equal to the diameter of the sphere. As the number of limes is equal to the number of parts of the circumference of the sphere, we get

Area of the surface = circumference x diameter.

Volume of the Sphere
To find the volume of a sphere Bhaskara II states the following method:
"Consider on the surface of the sphere pyramidal excavations, each of a base of an unit area having unit sides and of a depth equal to the radius, as many as the number of units of area in the surface. The apices of these pyramids meet at the centre of the sphere. The sum of the volumes of the pyramids is equal to the volume of the sphere. So it is proved (that the volume of a sphere is equal to the sixth part of the product of the surface area and diameter).

The above results are the nearest approach to the method of the integral calculus in Hindu Mathematics. It will be observed that the modern idea of the "limit of a sum" is not present. This idea, however, is of comparatively recent origin so that credit must be given to Bhaskara II for having used the same method as that of the integral calculus, although in a crude form.

Feedback   Contact email: