# Infinitesimal Calculus

INFINITESIMAL CALCULUS, i. The infinitesimal calculus is the body of rules and processes by means of which continuously varying magnitudes are dealt with in mathematical analysis. The name " infinitesimal " has been applied to the calculus because most of the leading results were first obtained by means of arguments about " infinitely small " quantities; the " infinitely small " or " infinitesimal " quantities were vaguely conceived as being neither zero nor finite but in some intermediate, nascent or evanescent, state. There was no necessity for this confused conception, and it came to be understood that it can be dispensed with; but the calculus was not developed by its first founders in accordance with logical principles from precisely defined notions, and it gained adherents rather through the impressiveness and variety of the results that could be obtained by using it than through the cogency of the arguments by which it was established. A similar statement might be made in regard to other theories included in mathematical analysis, such, for instance, as the theory of infinite series. Many, perhaps all, of the mathematical and physical theories which have survived have had a similar history a history which may be divided roughly into two periods: a period of construction, in which results are obtained from partially formed notions, and a period of criticism, in which the fundamental notions become progressively more and more precise, and are shown to be adequate bases for the constructions previously built upon them. These periods usually overlap. Critics of new theories are never lacking. On the other hand, as E. W. Hobson has well said, " pertinent criticism of fundamentals almost invariably gives rise to new construction." In the history of the infinitesimal calculus the Variable Quantities.

17th and 18th centuries were mainly a period of construction, the 19th century mainly a period of criticism.

I. Nature of the Calculus.

2. The guise in which variable quantities presented themselves to the mathematicians of the 17th century was that of the lengths of variable lines. This method of representing variable quantities dates from the 14th century, < riatl when it was employed by Nicole Oresme, who studied npnseatand afterwards taught at the College de Navarre in ff Paris from 1348 to 1361. He represented one of two variable quantities, e.g. the time that has elapsed since some epoch, by a length, called the "longitude," measured along a particular line; and he represented the other of the two quantities, e.g. the temperature at the instant, by a length, called the " latitude," measured at right angles to this line. He recognized that the variation of the temperature with the time was represented by the line, straight or curved, which joined the ends of all the lines of " latitude." Oresme's longitude and latitude were what we should now call the abscissa and ordinate. The same method was used later by many writers, among whom Johannes Kepler and Galileo Galilei may be mentioned. In Galileo's investigation of the motion of falling bodies (1638) the abscissa OA represents the time during which a body has been falling, and the ordinate AB represents the velocity acquired during that time (see fig. i). The velocity being proportional to the time, the " curve " obtained is a straight line OB, and Galileo showed FIG. i.

that the distance through which the body has fallen is represented by the area of the triangle OAB.

The most prominent problems in regard to a curve were the problem of finding the points at which the ordinate is a maximum or a minimum, the problem of drawing a tangent to The pnbthe curve at an assigned point, and the problem of lem* of determining the area of the curve. The relation of Maxima the problem of maxima and minima to the problem /H ln i ma of tangents was understood in the sense that maxima Tangents, or minima arise when a certain equation has equal "td Quadroots, and, when this is the case, the curves by which ntare *' the problem is to be solved touch each other. The reduction of problems of maxima and minima to problems of contact was known to Pappus. The problem of finding the area of a curve was usually presented in a particular form in which it is called the " problem of quadratures." It was sought to determine the area contained between the curve, the axis of abscissae and two ordinates, of which one was regarded as fixed and the other as variable. Galileo's investigation may serve as an example. In that example the fixed ordinate vanishes. From this investigation it may be seen that before the invention of the infinitesimal calculus the introduction of a curve into discussions of the course of any phenomenon, and the problem of quadratures for that curve, were not exclusively of geometrical import; the purpose for which the area of a curve was sought was often to find something which is not an area for instance, a length, or a volume or a centre of gravity.

3. The Greek geometers made little progress with the problem of tangents, but they devised methods for investigating the problem of quadratures. One of these methods was Qleek afterwards called the " method of exhaustions," and method*. the principle on which it is based was laid down in the lemma prefixed to the 12th book of Euclid's Elements as follows: " If from the greater of two magnitudes there be taken more than its half, and from the remainder more than its half, and so on, there will at length remain a magnitude less than the smaller of the proposed magnitudes." The method adopted by Archimedes was more general. It may be described as the enclosure of the magnitude to be evaluated between two others which can be brought by a definite process to differ from each other by less than any assigned magnitude. A simple example of its application is the 6th proposition of Archimedes' treatise On the Sphere and Cylinder, in which it is proved that the area contained between a regular polygon inscribed in a circle and a similar polygon circumscribed to the same circle can be made less than any assigned area by increasing thenumberof sides of the polygon. The methods of Euclid and Archimedes were specimens of rigorous limiting processes (see FUNCTION). The new problems presented by the analytical geometry and natural philosophy of the 17th century led to new limiting processes.

4. In the problem of tangents the new process may be described as follows. Let P, P' be two points of a curve (see fig. 2). Let x, y be the coordinates of P, and x+Ax, y+Ay those erea- o{ p/ The sym b o l Ax means " the difference of two Vat/oo. x < g an j tnere ; s a ijkg mean i n g f or the symbol Ay. The fraction Ay/Ax is the trigonometrical tangent of the angle which the secant PP' makes with the axis of x. Now let A* be continually diminished towards zero, so that P' continually approaches P. If the curve has a tangent at P the secant P P' approaches _ a limiting position (see 33 below). When * this is the case the fraction Ay/Ax tends to a limit, and this limit is the trigonometrical tangent of the angle which the tangent at P to the curve makes with the axis of .r. The limit is denoted by Z--" p IG 2 If the equation of the curve is of the form y=f(x) where/ is a functional symbol (see FUNCTION), then Ay /(r+Ax)-/(x)

Ax~ Ax ' and d x -< l "-A*-0 Ax The limit expressed by the right-hand member of this defining equation is often written /'(*).

and is called the " derived function " of /(x), sometimes the " derivative " or " derivate " of /(x). \\%en the function f(x) is a rational integral function, the division by Ax can be performed, and the limit is found by substituting zero for Ax in the quotient. For example, if /(x) =x ! , we have Ax and /'(*)= 2*.

The process of forming the derived function of a given function is called differentiation. The fraction Ay /Ax is called the " quotient of differences," and its limit dyjdx is called the " differential coefficient of y with respect to x." The rules for forming differential coefficients constitute the differential calculus.

The problem of tangents is solved at one stroke by the formation of the differential coefficient; and the problem of maxima and minima is solved, apart from the discrimination of maxima from minima and some further refinements, by equating the differential coefficient to zero (see MAXIMA AND MINIMA).

5. The problem of quadratures leads to a type of limiting process which may be described as follows: Let y=f(x) be the equation of a curve, and let AC and BD be the ordinates of the points C and D (see fig. 3). Let a, b be the abscissae of these points. Let the segment AB be divided into a number of segments by means of intermediate points such as M, and let MN be one such segment. Let PM and QN be those ordinates of the curve which have M and N as their feet. On MN as base describe two rectangles, of which the heights are the greatest and least values of y which correspond to points on the arc PQ of the curve. In fig. 3 these are the rectangles RM.SN. Let the sum of the areas of such rectangles as RM be formed, and likewise the sum of the areas of such rectangles as SN. When the number of the points such as M is increased without limit, and the lengths of all the FIG. 3.

segments such as MN are diminished without limit, these two sums of areas tend to limits. When they tend to the same limit the curvilinear figure ACDB has an area, and the limit is the measure of this area (see 33 below). The limit in question is the same whatever law may be adopted for inserting the points such as M between A and B, and for diminishing the lengths of the segments such as MN. Further, if P' is any point on the arc PQ, and P'M' is the ordinate of P', we may construct a rectangle of which the height is P'M' and the base is MN.and the limit of the sum of the areas of all such rectangles is the area of the figure as before. If x is the abscissa of P, x+Ax that of Q, x' that of P', the limit in question might be written where the letters o, b written below and above the sign of summation S indicate the extreme values of x. This limit is called " the definite integral of /(x) between the limits a and b," and the notation for it is The germs of this method of formulating the problem of quadratures are found in the writings of Archimedes. The method leads to a definition of a definite integral, but the direct application of it to the evaluation of integrals is in general difficult. Any process for evaluating a definite integral is a process of integration, and the rules for evaluating integrals constitute the integral calculus.

6. The chief of these rules is obtained by regarding the extreme ordinate BD as variable. Let now denote the abscissa of B.

The area A of the figure ACDB is represented by the ft Tneo.

integral I f(x)dx, and it is a function of . Let BD of In be displaced to B'D' so that becomes +A (see s ' oa> fig. 4). The area of the figure ACD'B' is represented by the integral | f(x)dx, and the increment AA 'of the area is given by Theorem ver- the formula which represents the area BDD'B'. between those two rectangles, having as a common base the segment BB , and as heights the greatest and least ordinates of points on the arc DD' of the curve. Let these heights be H and h. Then AA is intermediate between HA and AAJI, and the quotient of differences AA/A is intermediate between H and h. If the function /(x) is continuous at B (see FUNCTION), This area is intermediate FIG. 4.

then, as A is diminished without limit, H and h tend to BD, or /, as a limit, and we have The introduction of the process of differentiation, together with the theorem here proved, placed the solution of the problem of quadratures on a new basis. It appears that we can always find the area A if we know a function r (x) which has /(x) as its differential coefficient. If /(x) is continuous between a and b, we can prove that When we recognize a function F(x) which has the property expressed by the equation ^? =/(*)

we are said to integrate the function /(x), and F(x) is called the indefinite integral of /(x) with respect to x, and is written //(*)<**.

Differeatials - 7. In the process of 4 the increment Ay is not in general equal to the product of the increment Ax and the derived f unction /'(x). In general we can write down an equation of the form Ay=/'(x)Ax+R, in which R is different from zero when Ax is different from zero; and then we have not only lim -Az_0 R=0 ' but also We may separate Ay into two parts: the part /'(x)Ax and the part R. The partf'(x) Ax alone is useful for forming the differential coefficient, and it is convenient to give it a name. It is called the differential of /(x), and is written df(x), or ay when y is written for /(x). When this notation is adopted dx is written instead of Ax, and is called the " differential of x," so that we have df(x)=f'(x)dx.

Thus the differential of an independent variable such as x is a finite difference; in other words it is any number we please. The differential of a dependent variable such as y, or of a function of the independent variable x, is the product of the differential of x and the differential coefficient or derived function. It is important to observe that the differential coefficient is not to be defined as the ratio of differentials, but the ratio of differentials is to be defined as the previously introduced differential coefficient. The differentials are either finite differences, or are so much of certain finite differences as are useful for forming differential coefficients.

Again let F(x) be the indefinite integral of a continuous function f(x), so that we have When the points M of the process explained in 5 are inserted between the points whose abscissae are a and b, we may take them to be n I in number, so that the segment AB is divided into n segments. Let Xi, x 2 , . . .Xn-i be the abscissae of the points in order. The integral is the limit of the sum /(a) (*,-a)+/(*i) (*-*!) + ...+/(*) (x r+ i-x r )

every term of which is a differential of the (ormf(x)dx. Further the integral is equal to the sum of differences -F(*j)}+ ... -HF(* r+1 )-F(;v)} "tai Artifice for this sum is F(6) F(a). Now the difference F(av+i) F(x r ) is not equal to the differential 1(x r ) (x r+ ix r ), but the sum of the differences is equal to the limit of the sum of these differentials. The differential may be regarded as so much of the difference as is required to form the integral. From this point of view a differential is called a differential element of an integral, and the integral is the limit of the sum of differential elements. In like manner the differential element ydx of the area of a curve ( 5) is not the area of the portion contained between two ordinates, however near together, but is so much of this area as need be retained for the purpose of finding the area of the curve by the limiting process described.

8. The notation of the infinitesimal calculus is intimately bound up with the notions of differentials and sums of elements. The letter " d " is the initial letter of the word differentia (difference) 1 and the symbol "/ " is a conventionally written " S," the initial letter of the word summa (sum or whole). The notation was introduced by Leibnitz (see 25-27, below).

9. The fundamental artifice of the calculus is the artifice of forming differentials without first forming differential coefficients. From an equation containing x and y we can deduce a new equation, containing also Ax and Ay, by substituting x+Ax for x anc ' y+Ay f r y- M there is a differential coefficient of y with respect to x, then Ay can be expressed in the form <#>.Ax+R, where lim. A . c=0 (R/Aa:)=o l as in 7 above. The artifice consists in rejecting ab initio all terms of the equation which belong to R. We do not form R at all, but only <f>.Ax, or <t>. dx, which is the differential dy. In the same way, in all applications of the integral calculus to geometry or mechanics we form the element of an integral in the same way as the element of area y. dx is formed. In fig. 3 of 5 the element of area y. dx is the area of the rectangle RM. The actual area of the curvilinear figure PQNM is greater than the area of this rectangle by the area of the curvilinear figure PQR ; but the excess is less than the area of the rectangle PRQS, which is measured by the product of the numerical measures of MN and QR, and we have MN.QR lim. MN _ MN Thus the artifice by which differential elements of integrals are formed is in principle the same as that by which differentials are formed without first forming differential coefficients.

10. This principle is usually expressed by introducing the notion of orders of small quantities. If x, y are two variable numbers which are Ord rs ol connect ed together by any relation, and if when x tends to zero y also tends to zero, the fraction y/x may tend to a finite limit. In this case x and y are said to be " of the same order." When this is not the case we may have small quantities.

either r-o^ = - lim., In the former case y is said to be " of a lower order " than x; in the latter case y is said to be " of a higher order " than x. In accordance with this notion we may say that the fundamental artifice of the infinitesimal calculus consists in the rejection of small quantities of an unnecessarily high order. This artifice is now merely an incident in the conduct of a limiting process, but in the 17th century, when limiting processes other than the Greek methods for quadratures were new, the introduction of the artifice was a great advance.

n. By the aid of this artifice, or directly by carrying out Rules of t ' le appropriate limiting processes, we may obtain the Dlfferen- ru ' es &V which differential coefficients are formed. These tiatioo rules may be^classified as " formal rules " and " particular results." The formal rules may be stated as follows:

(i.) The differential coefficient of a constant is zero.

(ii.) For a sum u+v+ . . . +z, where u, v, . . .are functions of x, d(u+v+ ... +g)_du dv dz dx '~dx + dx + ' (iii.) For a product uv (iv.) For a quotient u/v d(u - (v.) For a function of a function, that is to say, for a function y expressed in terms of a variable z, which is itself expressed as a function of x, dy _dy_ dz 3* dz ' Hx In addition to these formal rules we have particular results as to the differentiation of simple functions. The most important results are written down in the following table:

nx"- 1 for all values of n loga* ar 1 loga* a 1 logo sin x cos x cos x sin x sin" 1 * (I -*)-! tan" 1 * (l+* 2 )-> Each of the formal rules, and each of the particular results in the table, is a theorem of the differential calculus. All functions (or rather expressions) which can be made up from those in the table by a finite number of operations of addition, subtraction, multiplication or division can be differentiated by the formal rules. All such functions are called explicit functions. In addition to these we have implicit functions, or such as are determined by an equation containing two variables when the equation cannot be solved so as to exhibit the one variable expressed in terms of the other. We have also functions of several variables. Further, since the derived function of a given function is itself a function, we may seek to differentiate it, and thus there arise the second and higher differential coefficients. We postpone for the present the problems of differential calculus which arise from these considerations. Again, we may have explicit functions which are expressed as the results of limiting operations, or by the limits of the results obtained by performing an infinite number of algebraic operations upon the simple functions. For the problem of differentiating such functions reference may be made to FUNCTION.

12. The processes of the integral calculus consist largely in transformations of the functions to be integrated into such indefinite forms that they can be recognized as differential co- integrals. efficients of functions which have previously been differentiated. Corresponding to the results in the table of II we have those in the following table:

/(*)

Jf(x)dx x n+l n + i for all values of n except i log e x e<" a~ l e" COS X sin x sin x cos x (a 2 -*")- sin-'j i a*+x 2 I . . x a tan "a The formal rules of 1 1 give us means for the transformation of integrals into recognizable forms. For example, the rule (ii.) for a sum leads to the result that the integral of a sum of a finite number of terms is the sum of the integrals of the several terms. The rule (iii.) for a product leads to the method of integration by parts. The rule (v.) for a function of a function leads to the method of substitution (see 48 below).

II. History.

13. The new limiting processes which were introduced in the development of the higher analysis were in the first instance Kepler's related to problems of the integral calculus. Johannes methods Kepler in his Astronomia nova . . . de motibus stellae ' ' texn ' Martis (1609) stated his laws of planetary motion, to the effect that the orbits of the planets are ellipses with the Sun at a focus, and that the radii vectores drawn from the Sun to the planets describe equal areas in equal times. From these statements it is to be concluded that Kepler could measure the areas of focal sectors of an ellipse. When he made out these laws there was no method of evaluating areas except the Greek methods. These methods would have sufficed for the purpose, but Kepler invented Jiis own method. He regarded the area as measured by the " sum of the radii " drawn from the focus, and he verified his laws of planetary motion by actually measuring a large number of radii of the orbit, spaced according to a rule, and adding their lengths.

He had observed that the focal radius vector SP (fig. 5) is equal to the perpendicular SZ drawn from S to the tangent at p to the auxiliary ^ circle, and he had further established the theorem which we should now express in the form the differential element of the area ASp as S/> turns about S, is equal to the product of SZ and the differential ad<t>, where a is the radius of the auxiliary circle, and <t> is the angle ACp, that is !A the eccentric angle of P on the ellipse. The area ASP bears to the area ASp the ratio of the minor to the major axis, a result known to FIG. 5. Archimedes. Thus Kepler's radii are spaced according to the rule that the eccentric angles of their ends are equidifferent, and his " sum of radii " is proportional to the expression which we should now write I (a -foe cos <t>)d4>.

where e is the eccentricity. Kepler evaluated the sum as proportional to <t>+e sin <t>.

Kepler soon afterwards occupied himself with the volumes of solids. The vintage of the year 1612 was extraordinarily abundant, and the question of the cubic content of wine casks was brought under his notice. This fact accounts for the title of his work, Nova stereometria doliorum; accessit stereometriae Archimedeae supplementum (1615). In this treatise he regarded solid bodies as being made up, as it were (veluli), of " infinitely " many " infinitely " small cones or " infinitely " thin disks, and he used the notion of summing the areas of the disks in the way he had previously used the notion of summing the focal radii of an ellipse.

14. In connexion with the early history of the calculus it must not be forgotten that the method by which logarithms Lo ar- were mvented (!6i4) was effectively a method of ithms. infinitesimals. Natural logarithms were not invented as the indices of a certain base, and the notation e for the base was first introduced by Euler more than a century after the invention. Logarithms were introduced as numbers which increase in arithmetic progression when other related numbers increase in geometric progression. The two sets of numbers were supposed to increase together, one at a uniform rate, the other at a variable rate, and the increments were regarded for purposes of calculation as very small and as accruing discontinuously.

15. Kepler's methods of integration, for such they must be called, were the origin of Bonaventura Cavalieri's theory of Cava- the summation of indivisibles. The notion of a Herfs continuum, such as the area within a closed curve, visible*. as bem S made U P f indivisible parts, " atoms " of area, if the expression may be allowed, is traceable to the speculations of early Greek philosophers; and although the nature of continuity was better understood by Aristotle and many other ancient writers yet the unsound atomic conception was revived in the 13th century and has not yet been finally uprooted. It is possible to contend that Cavalieri did not himself hold the unsound doctrine, but his writing on this point is rather obscure. In his treatise Geomelria indivisibilibus conlinuorum nova quadam ralione promola (1635) he regarded a plane figure as generated by a line moving so as to be always parallel to a fixed line, and a solid figure as generated by a plane moving so as to be always parallel to a fixed plane; and he compared the areas of two plane figures, or the volumes of two solids, by determining the ratios of the sums of all the indivisibles of which they are supposed to be made up, these indivisibles being segments of parallel lines equally spaced in the case of plane figures, and areas marked out upon parallel planes equally spaced in the case of solids. By this method Cavalieri was able to effect numerous integrations relating to the areas of portions of conic sections and the volumes generated by the revolution of these portions about various axes. At a later date, and partly in answer to an attack made upon him by Paul Guldin, Cavalieri published a treatise entitled Exerdtationes geometricae sex (1647), in which he adapted his method to the determination of centres of gravity, in particular for solids of variable density.

Among the results which he obtained is that which we should now write i ">+! integral).

He regarded the problem thus solved as that of determining the sum of the mth powers of all the lines drawn across a parallelogram parallel to one of its sides.

At this period scientific investigators communicated their results to one another through one or more intermediate persons. Such intermediaries were Pierre de Carcavy and Pater Marin Mersenne; and among the writers thus in communication were Bonaventura Cavalieri, Cavalier!. Christiaan Huygens, Galileo Galilei, Giles Personnier de Roberval, Pierre de Fermat, Evangelista Torricelli, and a little later Blaise Pascal; but the letters of Carcavy or Mersenne would probably come into the hands of any man who was likely to be interested in the matters discussed. It often happened that, when some new method was invented, or some new result obtained, the method or result was quickly known to a wide circle, although it might not be printed until after the lapse of a long time. When Cavalieri was printing his two treatises there was much discussion of the problem of quadratures. Roberval (1634) regarded an area as made up of " infinitely " many " infinitely " narrow strips, each of which may be considered to be a rectangle, and he had similar ideas in regard to lengths and volumes. He knew how to approximate to the quantity which we express by ( xTdx by the process of forming the sum o m + i m +2 m + . .. (n-i)" n m+l and he claimed to beabletoprove that thissumtendstoi/(m+i), as n increases for all positive integral values of m. The method of integrating x m by forming this sum was found also Permit's by Fermat (1636), who stated expressly that he method of arrived at it by generalizing a method employed by Iatf * r '- Archimedes (for the cases m = i and m = 2) in his books * on Conoids and Spheroids and on Spirals (see T. L. Heath, The Works of Archimedes, Cambridge, 1897). Fermat extended the result to the case where m is fractional (1644), and to the case where m is negative. This latte- extension and the proofs were given in his memoir, Proporlionis geometricae in quadrandis parabolis et hyperbolis us s, which appears to have received a final form before 1659, although not published until 1679. Fermat did not use fractional or negative indices, but he regarded his problems as the quadratures of parabolas and hyperbolas of various orders. His method was to divide the interval of integration into parts by means of intermediate points the abscissae of which are in geometric progression. In the process of s above, the points M must be chosen according to this rule. This restrictive condition being understood, we may say that Fermat's formulation of the problem of quadratures is the same as our definition of a definite integral.

The result that the problem of quadratures could be solved for any curve whose equation could be expressed in the form y = or in the form where none of the indices is equal to i, was used by John Wallis in his Arithmetica infinitorum (1655) as well as by Fermat (1659). The case in which m=i was that of the ordinary rectangular hyperbola; and Gregory of St Vincent in his Opus geometricum quadraturae circuli et sectionum coni (1647) had proved by the method of exhaustions that the area contained between the curve, one asymptote, and two ordinates parallel to the other asymptote, increases in arithmetic progression as the distance between the ordinates (the one nearer to the centre being kept fixed) increases in geometric progression. Fermat described his method of integration as a logarithmic method, and thus it is clear that the relation between the quadrature of the hyperbola and logarithms was understood although it was not expressed analytically. It was not very long before the relation was used for the calculation of logarithms by Nicolaus Mercator in his Logarilhmotechnia (1668). He began by writing the equation of the curve in the form y=i/(i+x), expanded this expression in powers of x by the method of division, and integrated it term by term in accordance with the well-understood rule for finding the quadrature of a curve given by such an equation as that written at the foot of p. 325.

By the middle of the 17th century many mathematicians could perform integrations. Very many particular results had Integra- been obtained, and applications of them had been tion before made to the quadrature of the circle and other conic the integral sections, and to various problems concerning the Calculus. i en g t hs O f curves, the areas they enclose, the volumes and superficial areas of solids, and centres of gravity. A systematic account of the methods then in use was given, along with much that was original on his part, by Blaise Pascal in his Leltres de Amos Deltonmlle sur qudques-unes de ses inventions en geotnetrie (1659).

16. The problem of maxima and minima and the problem of tangents had also by the same time been effectively solved. Fermat' a Oresme in the 14th century knew that at a point where methods of the ordinate of a curve is a maximum or a minimum DiHereo- its variation from point to point of the curve is slowest ; nation. ln( j K e pi er j n the Stereometria doliorum remarked that at the places where the ordinate passes from a smaller value to the greatest value and then again to a smaller value, its variation becomes insensible. Fermat in 1629 was in possession of a method which he then communicated to one Despagnet of Bordeaux, and which he referred to in a letter to Roberval of 1636. He communicated it to Rene Descartes early in 1638 on receiving a copy of Descartes's Gecmetrie (1637), and with it he sent to Descartes an account of his methods for solving the problem of tangents and for determining centres of gravity.

Fermat's method for maxima and minima is essentially our method. Expressed in a more modern notation, what he did was to begin by connecting the ordinate y and the abscissa x of a point of a curve by an equation which holds at all points of the curve, then to subtract the value of y in terms of * from the value obtained by substituting x + E for x, then to divide the difference by E, to put E=o in the quotient, and to equate the quotient to zero. Thus he differentiated with respect to x and equated the differential coefficient to zero.

FIG. 6. Fermat's method for solving the problem of tangents may be explained as follows:

Let (x, y) be the coordinates of a point P of a curve, (x', y'), those of a neighbouring point P' on the tangent at P, and let MM' = E (fig. 6). From the similarity of the triangles P'TM', PTM we have y':A-E=y:A, where A denotes the subtangent TM. The point P' being near the curve, we may substitute in the equation of the curve x E for x and (yA yE)/A for y. The equation of the curve is approximately satisfied. If it is taken to be satisfied exactly, the result is an equation of the form <j>(x, y, A, E)=o, the left-hand member of which is divisible by E. Omitting the factor E, and putting E=o in the remaining factor, we have an equation which gives A. In this problem of tangents also Fermat found the required result by a process equivalent to differentiation.

Fermat gave several examples of the application of his method; among them was one in which he showed that he could differentiate very complicated irrational functions. For such functions his method was to begin by obtaining a rational equation. In rationalizing equations Fermat, in other writings, used the device of introducing new variables, but he did not use this device to simplify the process of differentiation. Some of his results were published by Pierre Herigone in his Supplementum cursus mathemalici (1642). His communication to Descartes was not published in full until after his death (Fermat, Opera varia, 1679). Methods similar to Fermat's were devised by Rene de Sluse (1652) for tangents, and by Johannes Hudde (1658) for maxima and minima. Other methods for the solution of the problem of tangents were devised by Roberval and Torricelli, and published almost simultaneously in 1644. These methods were founded upon the composition of motions, the theory of which had been taught by Galileo (1638), and, less completely, by Roberval (1636). Roberval and Torricelli could construct the tangents of many curves, but they did not arrive at Fermat's artifice. This artifice is that which we have noted in 10 as the fundamental artifice of the infinitesimal calculus.

17. Among the comparatively few mathematicians who before 1665 could perform differentiations was Isaac Barrow. In his book entitled Lectiones opticae et geometricae, written apparently in 1663, 1664, and published in 1669, 1670, he gave a method of tangents like that of Roberval and Torricelli, compounding two velocities in the directions of the axes of x and y to obtain a resultant along the tangent to a curve. In an appendix to this book he gave another method which differs from Fermat's in the introduction of a differential equivalent to our dy as well as dx. Two neighbouring ordinates PM and QN of a curve (fig. 7) are regarded as containing an indefinitely small (indefinite parvum) arc, and PR is drawn parallel to the axis of x. The tangent PT at P is regarded as identical with the secant PQ, and the Barrow'* Differential Triangle.

FIG. 7.

position of the tangent is determined by the similarity of the triangles PTM, PQR. The increments QR, PR of the ordinate and abscissa are denoted by a and e; and the ratio of a to e is determined by substituting x+e for x and y+a for y in the equation of the curve, rejecting all terms which are of order higher than the first in a and e, and omitting the terms which do not contain a or e. This process is equivalent to differentiation. Barrow appears to have invented it himself, but to have put it into his book at Newton's request. The triangle PQR is sometimes called " Barrow's differential triangle."

The reciprocal relation between differentiation and integration ( 6) was first observed explicitly by Barrow in the book cited above. If the quadrature of a curve y=f(x) is known, so that the g,,,,,,,,,, area up to the ordinate x is given by F(x), the curve /nvers/on . y = F(x) can be drawn, and Barrow showed that the taeoretttf subtangent of this curve is measured by the ratio of its ordinate to the ordinate of the original curve. The curve y = F(x) is often called the " quadratrix " of the original curve; and the result has been called " Barrow's inversion-theorem." He did not use it as we do for the determination of quadratures, or indefinite integrals, but for the solution of problems of the kind which were then called " inverse problems of tangents." In these problems it was sought to determine a curve from some property of its tangent, e.g. the property that the subtangent is proportional to the square of the abscissa. Such problems are now classed under " differential equations." When Barrow wrote, quadratures were familiar and differentiation unfamiliar, just as hyperbolas were trusted while logarithms were strange. The functional notation was not invented till long afterwards (see FUNCTION), and the want of it is felt in reading all the mathematics of the i?th century.

18. The great secret which afterwards came to be called the " infinitesimal calculus " was almost discovered by Fermat, and still more nearly by Barrow. Barrow went farther than Fermat in the theory of differentiation, though not in the practice, for he compared two increments; he went farther in the theory of integration, for he obtained the inversiontheorem. The great discovery seems to consist partly in the recognition of the fact that differentiation, known to be a useful process, could always be performed, at least for the functions then known, and partly in the recognition T'rf7- 0/ of tne * act tnat l ^ e i nvers i n - tneorem could be coveiy applied to problems of quadrature. By these steps called the the problem of tangents could be solved once for all, intini- an( j t jj e operation of integration, as we call it, could be rendered systematic. A further step was necessary in order that the discovery, once made, should become accessible to mathematicians in general; and this step was the introduction of a suitable notation. The definite abandonment of the old tentative methods of integration in favour of the method in which this operation is regarded as thw inverse of differentiation was especially the work of Isaac Newton; the precise formulation of simple rules for the process of differentiation in each special case, and the introduction of the notation which has proved to be the best, were especially the work of Gottfried Wilhelm Leibnitz. This statement remains true although Newton invented a systematic notation, and practised differentiation by rules equivalent to those of Leibnitz, before Leibnitz had begun to work upon the subject, and Leibnitz effected integrations by the method of recognizing differential coefficients before he had had any opportunity of becoming acquainted with Newton's methods.

19. Newton was Barrow's pupil, and he knew to start with in 1664 all that Barrow knew, and that was practically all that was known about the subject at that time. His a- or igi na l thinking on the subject dates from the year tions. of the great plague (1665-1666), and it issued in the invention of the " Calculus of Fluxions," the principles and methods of which were developed by him in three tracts entitled De analyst per aequationes numero terminorum infinitas, Methodus fluxionum et serierum infinitarum, and De quadratura cunarum. None of these was published until long after they were written. The Analysis per aequationes was composed in 1666, but not printed until 1711, when it was published by William Jones. The Methodus fluxionum was composed in 1671 but not printed till 1736, nine years after Newton's death, when an English translation was published by John Colson. In Horsley's edition of Newton's works it bears the title Gcometria analytica. The Quadratura appears to have been composed in 1676, but was first printed in 1704 as an appendix to Newton's Oplicks.

20. The tract De Analyst per aequationes . . . was sent by Newton to Barrow, who sent it to John Collins with a request that , it might be made known. One way of making it known method of WO V' have been to print it in the Philosophical TransSeries actions of the Royal Society, but this course was not adopted. Collins made a copy of the tract and sent it to Lord Brouncker, but neither of them brought it before the Royal Society. The tract contains a general proof of Barrow's inversion-theorem which is the same in principle as that in 6 above. In this proof and elsewhere in the tract a notation is introduced for the momentary increment (momentum) of the abscissa or area of a curve; this " moment " is evidently meant to represent a moment of time, the abscissa representing time, and it is effectively the same as our differential element the thing that Fermat had denoted by E, and Barrow by e, in the case of the abscissa. Newton denoted the moment of the abscissa by o, that of the area z by ov. He used the letter v for the ordinate y, thus suggesting that his curve is a velocitytime graph such as Galileo had used. Newton gave the formula for the area of a curve v = x m (m*-l) in the form z = m+1 /(m + l). In the proof he transformed this formula to the form z n =c"x", where n and p are positive integers, substituted x+o for x and Z+OT for z, and expanded by the binomial theorem for a positive integral exponent, thus obtaining the relation from which he deduced the relation by omitting the equal terms z" and c"x p and dividing the remaining terms by o, tacitly putting = after division. This relation is the same as v = x m . Newton pointed out that, conversely, from the relation v-x m the relation z = x m+l /(m + l) follows. He applied his formula to the quadrature of curves whose ordinates can be expressed as the sum of a finite number of terms of the form ax* 1 -, and gave examples of its application to curves in which the ordinate is expressed by an infinite series, using for this purpose the binomial theorem for negative and fractional exponents, that is to say, the expansion of (i+x)" in an infinite series of powers of x. This theorem he had discovered ; but he did not in this tract state it in a general form or give any proof of it. He pointed out, however, how it may be used for the solution of equations by means of infinite series. He observed also that all questions concerning lengths of curves, volumes enclosed by surfaces, and centres of gravity, can be formulated as problems of quadratures, and can thus be solved either in finite terms or by means of infinite series. In the Quadratura (1676) the method of integration which is founded upon the inversiontheorem was carried out systematically. Among other results there given is the quadrature of curves expressed by equations of the Form y = x n (a+bx m ) p ; this has passed into text-books under the title " integration of binomial differentials " (see 49). Newton announced the result in letters to Collins and Oldenburg of 1676.

21. In the Methodus fluxionum (1671) Newton introduced his characteristic notation. He regarded variable quantities as generated by the motion of a point, or line, or plane, and called Ne ^ aa , s the generated quantity a " fluent " and its rate of genera- , ethod of tion a " fluxion." The fluxion of a fluent x is represented Fluxlons by x, and its moment, or " infinitely " small increment accruing in an " infinitely " short time, is represented by xo. The problems of the calculus are stated to be (i.) to find the velocity at any time when the distance traversed is given; (ii.) to find the distance traversed when the yelocity is given. The first of these leads to differentiation. In any rational equation containing * and y the expressions x-j-xo and y+yo_ are to be substituted for x and y, the resulting equation is to be divided by o, and afterwards o is to be omitted. In the case of irrational functions, or rational functions which are not integral, new variables are introduced in such a way as to make the equations contain rational integral terms only. Thus Newton's rules of differentiation would be in our notation the rules (i.), (ii.), (v.) of II, together with the particular result which we write ^p = mx m ~ l , (m integral).

a result which Newton obtained by expanding (x+xo) m by the binomial theorem. The second problem is the problem of integration, and Newton's method for solving it was the method of series founded upon the particular result which we write x m dx-- Newton added applications of his methods to maxima and minima, tangents and curvature. In a letter to Collins of date 1672 Newton stated that he had certain methods, and he described certain results which he had found by using them. These methods and results are those which are to be found in the Methodus fluxionum ; but the letter makes no mention of fluxions and fluents or of the characteristic notation. The rule for tangents is said in the letter to be analogous to de Sluse's, but to be applicable to equations that contain irrational terms.

22. Newton gave the fluxional notation also in the tract De Quadratura cunarum (1676), and he there added to it notation for the higher differential coefficients and for indefinite . ... integrals, as we call them. Just as.r, y, z, . . . are fluents of which x, y, t, . . . are the fluxions, so x, y, z, . . . can t '" ''[ be treated as fluents of which the fluxions may be denoted ., ' . by *, y', 2, . . . In like manner the fluxions of these may be denoted by X, y, Z, . . and so on. Again x, y, z, . . . may be regarded as fluxions of which the fluents may be denoted by x, y, i,. . . , and these again as fluxions of other quantities denoted by x,y,z, . . . and so on. No use was made of the notation x , x, , . . in the course of the tract. The first publication of the fluxional notation was made by Wallis in the second edition of his Algebra (1693) in the form of extracts from communications made to him by Newton in 1692. In this account of the method the symbols o, x, X, . . . occur, but not the symbols x, x, . . . Wallis's treatise also contains Newton's formulation of the problems of the calculus in the words Data aequatione fluentes quotcumque quantitates involvente fluxiones imienire et vice versa (" an equation containing any number of fluent quantities being given, to find their fluxions and vice versa "). In the Philosophiae naturalis principia mathematica (1687), commonly called the " Principia," the words " fluxion " and " moment " occur in a lemma in the second book; but the notation which is characteristic of the calculus of fluxions is nowhere used.

23. It is difficult to account for the fragmentary manner of publication of the Fluxional Calculus and for the long delays which took place. At the time (1671) when Newton Ketanled composed the Methodus fluxionum he contemplated Publicsbringing out an edition of Gerhard Kinckhuysen's ttonoithe treatise on algebra and prefixing his tract to this treatise. In the same year his " Theory of Light and Colours " was published in the Philosophical Transactions, and the opposition which it excited led to the abandonment of the project with regard to fluxions. In 1680 Collins sought the assistance of the Royal Society for the publication of the tract, and this was granted in 1682. Yet it remained unpublished. The reason is unknown; but it is known that about 1679, 1680, Newton took up again the studies in natural philosophy which he had intermitted for several years, and that in 1684 he wrote the tract De motu which was in some sense a first draft of the Principia, and it may be conjectured that the fluxions were held over until the Principia should be finished. There is also reason to think that Newton had become dissatisfied with the arguments about infinitesimals on which his calculus was based. In the preface to the De quadrature, curvarum (1704), in which he describes this tract as something which he once wrote (" olim scripsi ") he says that there is no necessity to introduce into the method of fluxions any argument about infinitely small quantities; and in the Principia (1687) he adopted instead of the method of fluxions a new method, that of " Prime and Ultimate Ratios." By the aid of this method it is possible, as Newton knew, and as was afterwards seen by others, to found the calculus of fluxions on an irreproachable method of limits. For the purpose of explaining his discoveries in dynamics and astronomy Newton used the method of limits only, without the notation of fluxions, and he presented all his results and demonstrations in a geometrical form. There is no doubt that he arrived at most of his theorems in the first instance by using the method of fluxions. Further evidence of Newton's dissatisfaction with arguments about infinitely small quantities is furnished by his tract Methodus differentialis, published in 1711 by William Jones, in which he laid the foundations of the " Calculus of Finite Differences."

26. In the Ada erudilorum of 1684 Leibnitz published a short memoir entitled Nova metkodus pro maximis et minimis, itemque tangentibus, quae nee fractas nee irralionales Leibnitz's quantitates moratur, et singulare pro illis calculi genus. DifferIn this memoir the differential dx of a variable x, """' considered as the abscissa of a point of a curve, is said Cajcalus - to be an arbitrary quantity, and the differential dy of a related variable y, considered as the ordinate of the point, is defined as a quantity which has to dx the ratio of the ordinate to the subtangent, and rules are given for operating with differentials. These are the rules for forming the differential of a constant, a sum (or difference), a product, a quotient, a power (or root). They are equivalent to our rules (i.)-(iv.) of 1 1 and the particular result The rule for a function of a function is not stated explicitly but is illustrated by examples in which new variables are introduced, in much the same way as in Newton's Methodus fluxionum. In connexion with the problem of maxima and minima, it is noted that the differential of y is positive or negative according as y increases or decreases when x increases, and the discrimination of maxima from minima depends upon the sign of ddy, the differential of dy. In connexion with the problem of tangents the differentials are said to be proportional to the momentary increments of the abscissa and ordinate. A tangent is denned as a line joining two " infinitely " near points of a curve, and the " infinitely " small distances (e.g., the distance between the feet of the ordinates of such points) are said to be expressible by means of the differentials (e.g.,dx). The method is illustrated by a few examples, and one example is given of its application to " inverse problems of tangents." Barrow's inversion-theorem and its application to quadratures are not mentioned. No proofs are given, but it is stated that they can be obtained easily by any one versed in such matters. The new methods in regard to differentiation which were contained in this memoir were the use of the second differential for the discrimination of maxima and minima, and the introduction of new variables for the purpose of differentiating complicated expressions. A greater novelty was the use of a letter (d), not as a symbol for a number or magnitude, but as a symbol of operation. None of these novelties account for the far-reaching effect which this memoir has had upon the development of mathematical analysis. This effect was a consequence of the simplicity and directness with which the rules of differentiation were stated. Whatever indistinctness might be felt to attach to the symbols, the processes for solving problems of tangents and of maxima and minima were reduced once for all to a definite routine.

27. This memoir was followed in 1686 by a second, entitled De Geometria recondita et analysi indivisibilium atque infinitorum, Deveio * n wn ' crl Leibnitz described the method of using his meat new differential calculus for the problem of quadratures. of the This was the first publication of the notation fydx. Calculus. -p ne new me thod was called calculus summalorius. The brothers Jacob (James) and Johann (John) Bernoulli were able by 1690 to begin to make substantial contributions to the development of the new calculus, and Leibnitz adopted their word " integral " in 1695, they at the same time adopting his symbol "/." In 1696 the marquis de 1'Hospital published the first treatise on the differential calculus with the title Analyse des infiniment petits pour I' intelligence des lignes courbes. The few references to fluxions in Newton's Principia (1687) must have been quite unintelligible to the mathematicians of the time, and the publication of the fluxional notation and calculus by Wallis in 1693 was too late to be effective. Fluxions had been supplanted before they were introduced.

The differential calculus and the integral calculus were rapidly developed in the writings of Leibnitz and the Bernoullis. Leibnitz (1695) was the first to differentiate a logarithm and an exponential, and John Bernoulli was the first to recognize the property possessed by an exponential (a 1 ) of becoming infinitely great in comparison with any power (*") when x is increased indefinitely. Roger Cotes (1722) was the first to differentiate a trigonometrical function. A great development of infinitesimal methods took place through the founding in 1696-1697 of the " Calculus of Variations " by the brothers Bernoulli.

" The differential method is one and the same with the method of fluxions, excepting the name and mode of notation; Mr Leibnitz calling those quantities differences which Mr Newton calls moments or fluxions, and marking them with the letter d, a mark not used by Mr Newton. And therefore we take the proper question to be, not who invented this or that method, but who was the first inventor of the method; and we believe that those who have reputed Mr Leibnitz the first inventor, knew little or nothing of his correspondence with Mr Collins and Mr Oldenburg long before; nor of Mr Newton's having that method above fifteen years before Mr. Leibnitz began to publish it in the Acta Eruditorum of Leipzig. For which reasons we reckon Mr Newton the first inventor, and are of opinion that Mr Keill, in asserting the same, has been no ways injurious to Mr Leibnitz."

The report with the letters and other documents was printed (1712) under the title Commercium Epistolicum D. Johannis Collins et aliorum de analysi promota, jussu Societatis Regiae in lucent editum, not at first for publication. An account of the contents of the Commercium Epistolicum was printed in the Philosophical Transactions for 1715. A second edition of the Commercium Epistolicum was published in 1722. The dispute was continued for many years after the death of Leibnitz in 1716. To translate the words of Moritz Cantor, it " redounded to the discredit of all concerned."

29. One lamentable consequence of the dispute was a severance of British methods from continental ones. In Great Britain it became a point of honour to use fluxions and other Newtonian methods, while on the continent the notation of Leibnitz was universally adopted. This tiaeatai severance did not at first prevent a great advance in Schools of mathematics in Great Britain. So long as attention ^ was directed to problems in which there is but one independent variable (the time, or the abscissa of a point of a curve), and all the other variables depend upon this one, the fluxional notation could be used as well as the differential and integral notation, though perhaps not quite so easily. Up to about the middle of the 18th century important discoveries continued to be made by the use of the method of fluxions. It was the introduction of partial differentiation by Leonhard Euler (1734) and Alexis Claude Clairaut (1739), and the developments which followed upon the systematic use of partial differential coefficients, which led to Great Britain being left behind; and it was not until after the reintroduction of continental methods into England by Sir John Herschel, George Peacock and Charles Babbage in 1815 that British mathematics began to flourish again. The exclusion of continental mathematics from Great Britain was not accompanied by any exclusion of British mathematics from the continent. The discoveries of Brook Taylor and Colon Maclaurin were absorbed into the rapidly growing continental analysis, and the more precise conceptions reached through a critical scrutiny of the true nature of Newton's fluxions and moments stimulated a like scrutiny of the basis of the method of differentials.

31. Defective modes of expression did not hinder constructive work. It was the great merit of Leibnitz's symbolism that a mathematician who used it knew what was to be done in order to formulate any problem analytically, even though he might not be absolutely clear as to the basls of proper interpretation of the symbols, or able to render a satisfactory account of them. While new and varied results were promptly obtained by using them, a long time elapsed before the theory of them was placed on a sound basis. Even after Cauchy had formulated his theory much remained to be done, both in the rapidly growing department of complex variables, and in the regions opened up by the theory of expansions in trigonometric series. In both directions it was seen that rigorous demonstration demanded greater precision in regard to fundamental notions, and the requirement of precision led to a gradual shifting of the basis of analysis from geometrical intuition to arithmetical law. A sketch of the outcome of this movement the " arithmetization of analysis," as it has been called will be found in FUNCTION. Its general tendency has been to show that many theories and processes, at first accepted as of general validity, are liable to exceptions, and much of the work of the analysts of the latter half of the 19th century was directed to discovering the most general conditions in which particular processes, frequently but not universally applicable, can be used without scruple.

III. Outlines of the Infinitesimal Calculus.

32. The general notions of functionality, limits and continuity are explained in the article FUNCTION. Illustrations of the more immediate ways in which these notions present themselves in the development of the differential and integral calculus will be useful in what follows.

33. Let y be given as a function of x, or, more generally, let x and y be given as functions of a variable t. The first of these cases is included in the second by putting x = t. If certain conditions are satisfied the aggregate of the points determined by the functional relations form a curve. The first condition is that the aggregate of the values of / to which values of x and y correspond must be continuous, or, in other words, that these values must consist of all real numbers, or of all those real numbers which lie between assigned extreme numbers. When this condition is satisfied the points are " ordered," and their order is determined by the order of the numbers t, supposed to be arranged in order of increasing or decreasing magnitude; also there are two senses of description of the curve, according as t is taken to increase or to diminish. The second condition is that the aggregate of the points which are determined by the functional relations must be " continuous." This condition means that, if any point P determined by a value of / is taken, and any distance 6, however small, is chosen, it is possible to find two points Q, Q' of the aggregate which are such that (i.) P is between Q and Q', (ii.) if R, R are any points between Q and Q' the distance RR' is less than 8. The meaning of the word " between " in this statement is fixed by the ordering of the points. Sometimes additional conditions are imposed upon the functional relations before they are regarded as defining a curve. An aggregate of points which satisfies the two conditions stated above is sometimes called a " Jordan curve." It by no means follows that every curve of this kind has a tangent. In order that the curve may have a tangent T ats . at P it is necessary metrical limits.

that, if any angle o, however small, is specified, a distance 6 can be found such that when P is between Q and Q', and PQ and PQ' are less than , the angle RPR' is less than o for all pairs of points R, R' which are between P and Q, or between P and Q' (fig. 8). When this condition is satisfied y is a function of * which has a differential coefficient. The only way of FIG. 8.

finding out whether this condition is satisfied or not is to attempt to form the differential coefficient. If the quotient of differences Ay/A* has a limit when Ax tends to zero, y is a differentiable function of x, and the limit in question is the differential coefficient. The derived function, or differential coefficient, of a function f(x) is always defined by the formula Rules for the formation of differential coefficients in particular cases have been given in 11 above. The definition of a differential coefficient, and the rules of differentiation, are quite independent of any geometrical interpretation, such as that concerning tangents to a curve, and the tangent to a curve is properly defined by means of the differential coefficient of a function, not the differential coefficient by means of the tangent.

It may happen that the limit employed in defining the differential coefficient has one value when h approaches zero through positive Proms- va ' ues * an d a different value when h approaches zero J through negative values. The two limits are then called Regressive tne P r g ress i ve " a d ' regressive " differential coDlfferen- efficients. In applications to dynamics, when * denotes tlal Co- a coordinate and / the time, dx/dt denotes a velocity. If efficients. tne ve 'city is changed suddenly the progressive differ- ential coefficient measures the velocity just after the change, and the regressive differential coefficient measures the velocity just before the change. Variable velocities are properly defined by means of differential coefficients.

AH geometrical limits may be specified in terms similar to those employed in specifying the tangent to a curve; in difficult cases Areas they must be so specified. Geometrical intuition may fail to answer the question of the existence or non-existence of the appropriate limits. In the last resort the definitions of many quantities of geometrical import must be analytical, not geometrical. As illustrations of this statement we may take the definitions of the areas and lengths of curves. We may not assume that every curve has an area or a length. To find out whether a curve has an area or not, we must ascertain whether the limit expressed by fydx exists. _ When the limit exists the curve has an area. The definition of the integral is quite independent of any geometrical interpretation. The length of a curve again is denned by means of a limiting process. Let P, Q be two points of a curve, and Ri, R 2 , . . . R^_i a set of intermediate points of the curve, supposed to be described in the sense in which Q comes after P. The points R are supposed to be reached successively in the order of the suffixes when the curve is described in this sense. We form a sum of lengths of chords PR 1 +R,R 2 + . . . +R^Q.

If this sum has a limit when the number of the points R is increased indefinitely and the lengths of all the chords are diminished indeLenrths fi mte 'y. tn ' s limit is the length of the arc PQ. The limit of Curves. ' s tn . e same whatever law may be adopted for inserting the intermediate points R and diminishing the lengths of the chords. It appears from this statement that the differential element of the arc of a curve is the length of the chord joining two neighbouring points. In accordance with the fundamental artifice for forming differentials ( 9, 10), the differential element of arc ds may be expressed by the formula ds = V [(<**)*+ (<ty)'l, of which the right-hand member is really the measure of the distance between two neighbouring points on the tangent. The square root must be taken to be positive. We may describe this differential element as being so much of the actual arc between two neighbouring points as need be retained for the purpose of forming the integral expression for an arc. This is a description, not a definition, because the length of the short arc itself is only definable by means of the integral expression. Similar considerations to those used in defining the areas of plane figures and the lengths of plane curves are applicable to the formation of expressions for differential elements of volume or of the areas of curved surfaces.

34. In regard to differential coefficients it is an important theorem that, if the derived f unction f'(x) vanishes at all points of an interval, Constant* tne function /(x) is constant in the interval. It follows of late- tnat ' 'f two functions have the same derived function gration. !- nev can on ' v differ by a constant. Conversely, indefinite integrals are indeterminate to the extent of an additive constant.

35- .The differential coefficient dy/dx, or the derived function / (*), is itself a function of x, and its differential coefficient is denoted Higher tyjf'W or ^P**' In the second of these notations Dlffereu- J. Is r . e K ar ded as the symbol of an operation, that of tlal Co- differentiation with respect to x, and the index 2 means efficients. tnat tne operation is repeated. In like manner we may express the results of n successive differentiations by /"> (x) or by dy/dx. When the second differential coefficient exists, or the first is differentiable, we have the relation /' (*) = n~ ,._/(+*)-/(*)+/(-) (L)

The limit expressed by the right-hand member of this equation may exist in cases in which f(x) does not exist or is not differentiable. The result that, when the limit here expressed can be shown to vanish at all points of an interval, then/(x) must be a linear function of x in the interval, is important.

The relation (i.) is a particular case of the more general relation As in the case of relation (i.) the limit expressed by the right-hand member may exist although some or all of the derived functions f'(x), f*(x), . . .f- n ~^(x) do not exist.

Corresponding to the rule iii. of 1 1 we have the rule for forming the wth differential coefficient of a product in the form d n (uv) d n v i du d n ~ l v_>n(n i) d**u d n ~^ z v, \d n u where the coefficients are those of the expansion of (i +x) in powers of x (n being a positive integer). The rule is due to Leibnitz, (1695).

Differentials of higher orders may be introduced in the same way as the differential of the first order. In general when y=f(x), the wth differential d n y is defined by the equation in which dx is the (arbitrary) differential of x.

When d/dx is regarded as a single symbol of operation the symbol f...dx represents the inverse operation. If the former is denoted by D, the latter may be denoted by Dy l . D" means that _ the operation D is to be performed n times in succession ; J 60 ' 8 D~* that the operation of forming the indefinite integral ^ "*"*' is to be performed n times in succession. Leibnitz's course of thought ( 24) naturally led him to inquire after an interpretation of D" where n is not an integer. For an account of the researches to which this inquiry gave rise, reference may be made to the article by A. Voss in Ency. d. math. Wiss. Bd. ii. A, 2 (Leipzig, 1889). The matter is referred to as " fractional "or " generalized" differentiation.

36. After the formation of differential coefficients the most important theorem of the differential calculus is the theorem of intermediate value (" theorem of mean _. value," " theorem of finite increments," " Rolle's theorem," are other names for it). This theorem may be explained as follows: Let A, B be two points of a curve y=f(x) (fig. 9). Then there is a point P between A and B at which the tangent is parallel to the secant AB. This theorem is expressed analytically in the statement that if f'(x) is continuous between a and b, there is a value Xi of x between o and b which has the property expressed by the equation FIG. 9.

The value Xi can be expressed in the form 0+6(60) where 6 is a number between o and I.

A slightly more general theorem was gjven by Cauchy (1823) to the effect that, if f'(x) and F'(*) are continuous between x = a and x = b, then there is a number 6 between o and i which has the property expressed by the equation F(6)-F(q)_F'|q+9(&-q)} /(*) ~/(o) f'\a+6(b~a)\' The theorem expressed by the relation (i.) was first noted by Rolle (1690) for the case where /(x) is a rational integral function which vanishes when x = a and also when x = b. The general theorem was given by Lagrange (1797). Its fundamental importance was first recognized by Cauchy (1823). It may be observed here that the theorem of integral calculus expressed by the equation follows at once from the definition of an integral and the theorem of intermediate value.

The theorem of intermediate value may be generalized in the statement that, if /(*) and all its differential coefficients up to the nth inclusive are continuous in the interval between x = a and* = 6, then there is a number 9 between o and I which has the property expressed by the equation 17. This theorem provides a means for computing the values of a function at points near to an assigned point when the value of the function and its differential coefficients at the assigned _ . point are known. The function is expressed by a termin- J ated series, and, when the remainder tends to zero as n ""' increases, it may be transformed into an infinite series. The theorem was first given by Brook Taylor in his Methodus Incrementorum (1717) as a corollary to a theorem concerning finite differences. Taylor gave the expression for f(x+z) in terms of f(x), /'(*),... as an infinite series proceeding by powers of z. His notation was that appropriate to the method of fluxions which he used. This rule for expressing a function as an infinite series is known as Taylor's theorem. The relation (i.), in which the remainder after n terms is put in evidence, was first obtained by Lagrange (1797). Another form of the remainder was given by Cauchy (1823) viz., The conditions of validity of Taylor's expansion in an infinite series have been investigated very completely oy A. Pringsheim (Math. Ann. Bd. xliv., 1894). It is not sufficient that the function and all its differential coefficients should be finite at x = a ; there must be a neighbourhood of a within which Cauchy's form of the remainder tends to zero as n increases (cf. FUNCTION).

An example of the necessity of this condition is afforded by the function /(*) which is given by the equation The sum of the series power series.

is the same as that of the series It is easy to prove that this is less than e~ l when x lies between o and i, and also that/(jc) is greater than e' 1 when x = i/V3- Hence the sum of the series (i.) is not equal to the sum of the series (ii.).

The particular case of Taylor's theorem in which = is often called Maclaurin's theorem, because it was first explicitly stated by Colin Maclaurin in his Treatise of Fluxions (1742). Maclaurin like Taylor worked exclusively with the fluxional calculus.

Examples of expansions in series had been known for some time. The series for log (i+x) was obtained by Nicolaus Mercator (1668) by expanding (i + x)~ l by the method of algebraic division, and integrating the series term by term. He regarded his result as a " quadrature of the hyperbola." Newton (1669) obtained the expansion of sin" 1 * by expanding (i-x 2 )~i by the binomial theorem and Integrating the series term by term. James Gregory (1671) gave the series for tan" 1 *:. Newton also obtained the series for sin x, cos x, and e x by reversion of series (1669). The symbol e for the base of the Napierian logarithms was introduced by Euler (1739). All these series can be obtained at once by Taylor's theorem. James Gregory found also the first few terms of the series for tan x and sec x; the terms of these series may be found successively by Taylor's theorem, but the numerical coefficient of the general term cannot be obtained in this way.

Taylor's theorem for the expansion of a function in a power series was the basis of Lagrange's theory of functions, and it is fundamental also in the theory of analytic functions of a complex variable as developed later by Karl Weierstrass. It has also numerous applications to problems of maxima and minima and to analytical geometry. These matters are treated in the appropriate articles.

The forms of the coefficients in the series for tan x and sec x can be expressed most simply in terms of a set of numbers introduced by James Bernoulli in his treatise on probability entitled Ars Conjectandi (1713). These numbers BI, B 2 , . . . called Bernoulli's numbers, are the coefficients so denoted in the formula e'-l and they are connected with the sums of powers of the reciprocals of the natural numbers by equations of the type The function 2 2\ has been called Bernoulli's function of the with order by J. L. Raabe (Crelle's J. /. Math. Bd. xlii., 1851). Bernoulli's numbers and functions are of especial importance in the calculus of finite differences (see the article by D. Seliwanoff in Ency. d. math. Wiss. Bd. i., E., 1901).

When * is given in terms of y by means of a power series of the form .x= y (C<,+C iy +Crf+...) (Co=t=o)=y/ (y), say, there arises the problem of expressing y as a power series in x. This problem is that of reversion of series. It can be shown that provided the absolute value of x is not too great, To this problem is reducible that of expanding y in powers of * when x and y are connected by an equation of the form y = a+xf(y), for which problem Lagrange (1770) obtained the formula For the history of the problem and the generalizations of Lagrange's result reference may be made to O. Stolz, Grundzuge d. Diff. u. Int. Rechnung, T. 2 (Leipzig, 1896).

38. An important application of the theorem of intermediate value and its generalization can be made to the problem of evaluating certain limits. If two functions <(*) and t(x) both . ,, vanish at x=a, the fraction <*>(*) AK*) may have a finite imit at a. This limit is described as the limit of an indeterminate form." Such indeterminate forms were considered first by de 1'Hospital (1696) to whom the problem of evaluating the limit presented itself in the form of tracing the curve y=#(*)AK*) near the ordinate x=a, when the curves y =<K*0 and y = <l/(x) both cross the axis of x at the same point as this ordinate. In fig. 10 PA and QA represent short arcs of the curves \$,' \$, chosen so that P and Q have the same abscissa. The value of the ordinate of the corresponding point R of the compound curve is given by the ratio of the ordinates PM, QM. _ De 1'Hospital treated PM and QM // M -r as " infinitesimal," so that the equations PM :AM ='(o) and QM :AM =t'(a) could FIG. 10.

be assumed to hold, and he arrived at the result that the " true value " of <t>(a)l\j/(a) is <j>' '(a) l\j/' '(a) . It can be proved rigorously that, if '/''(x) does not vanish at x = a, while <t>(a) =o and (f-(a) =o, then It can be proved further if that <t> m (x) and t n (x) are the differential coefficients of lowest order of <j>(x) and <j/ (x) which do not vanish at x = a, and if m = n, then lim.,.

If m>n the limit is zero; but if m<n the function represented by the quotient <t>(x)lt(x) " becomes infinite " at x = a. If the value of the function at x =a is not assigned by the definition of the function, the function does not exist at x = a, and the meaning of the statement that it " becomes infinite " is that it has no finite limit. The statement does not mean that the function has a value which we call infinity. There is no such value (see FUNCTION).

Such indeterminate forms as that described above are said to be of the form o/o. Other indeterminate forms are presented in the form ox oo, or i, or oo/oo, or oo-oo. The most notable of the forms i is lim.i=o(i+^) 1 / 1 , which is e. The case in which <j>(x) and tf/(x) both tend to become infinite at x = a is reducible to the case in which both the functions tend to become infinite when x is increased indefinitely. If <t>'(x) and ^'(x) have determinate finite limits when x is increased indefinitely, while <t>(x) and ^(x) are determinately (positively or negatively) infinite, we have the result expressed by the equation lim..

lim..

.*'<*)

For the meaning of the statement that <t> (x) and \l/ (x) are determinately infinite reference may be made to the article FUNCTION. The evaluation of forms of the type =0/00 leads to a scale of increasing " infinities," each being infinite in comparison with the preceding. Such a scale is log x,...x, x"*,...x n ,...e x ,...x x ; each of the limits expressed by such forms as \\rn.x-x <t>(x}l\l/(x), where <t>(x) precedes ^/(x) in the scale, is zero. The construction of such scales, along with the problem of constructing a complete scale, was discussed in numerous writings by Paul du Bois-Reymond (see in particular, Math. Ann. Bd. xi., 1877). For the general Croblem of indeterminate forms reference may be made to the article y A. Pringsheim in Ency. d. math. Wiss. Bd. ii., A. i_(i89g). Forms of the type 0/0 presented themselves to early writers on analytical geometry in connexion with the determination of the tangents at a double point of a curve; forms of the type oo/oo presented themselves in like manner in connexion with the determination of asymptotes of curves. The evaluation of limits has innumerable applications in all parts of analysis. Cauchy's Analyse algbbrique (1821) was an epoch-making treatise on limits.

If a function <(*) becomes infinite at x = a, and another function iHx) also becomes infinite at x=a in such a way that <t>(x)/<l/(x) has a finite limit C, we say that <j>(x) and ^(x) become " infinite of the same order." We may write <t>(x) =C<l/(x)+<tn(x), where Hm.i=o0iW/<A(*) =o, and thus <#>i(x) is of a lower order than* (A;); it may be finite or infinite at x = a. If it is finite, we describe C<l>(x)

xiv. 1 8 as the " infinite part " of <#>(*). The resolution of a function which becomes infinite into an infinite part and a finite part can often be effected by taking the infinite part to be infinite of the same order as one of the functions in the scale written above, or in some more comprehensive scale. This resolution is the inverse of the process of evaluating an indeterminate form of the type oo oo .

For example lim. a ._ {(e 1 l)" 1 x' 1 } is finite and equal to =, and the function (e* i)" 1 x"" 1 can be expanded in a power series in x.

39. The nature of a function of two or more variables, and the meaning, to be attached to continuity and limits in respect of such functions, have been explained under FUNCTION. The Functions theorems of differential calculus which relate to such severa/ faggfaM are ( n general the same whether the number variables. Q j var i a ki es ; s two or anv greater number, and it will generally be convenient to state the theorems for two variables.

40. Let u or f (x f , y) denote a function of two variables * and y. If we regard y as constant, u or/ becomes a function of one variable x, and we may seek to differentiate it with respect to *. Partial jj t ^ e f unc tion of x is differentiable, the differential differen- coefficient which is formed in this way is called the nation. p art j a i differential coefficient " of u or/ with respect to x, and is denoted by ^ or g*. The symbol " d " was appropriated for partial differentiation by C. G. J. Jacobi (1841). It had before been written indifferently with " d as a symbol of differentiation.

Euler had written " (+] " for the partial differential coefficient of / with respect to x. Sometimes it is desirable to put in evidence the variable which is treated as constant, and then the partial differential coefficient is written " (-, )" or "(a ) " This course is often \ax/ \ox; a adopted by writers on Thermodynamics. Sometimes the symbols d or d are dropped, and the partial differential coefficient is denoted by MI orfx. As a definition of the partial differential coefficient we have the formula In the same way we may form the partial differential coefficient with respect to y by treating x as a constant.

The introduction of partial differential coefficients enables us to solve at once for a surface a problem analogous to the problem of tangents for a curve ; and it also enables us to take the first step in the solution of the problem of maxima and minima for a function of several variables. If the equation of a surface is expressed in the form z=f(x, y), the direction cosines of the normal to the surface d/ d/ at any point are in the ratios -f- : -: = I. If / is a maximum or a minimum at (x, y), then dfjdx and df/dy vanish at that point.

In applications of the differential calculus to mathematical physics we are in general concerned with functions of three variables x, y, z, which represent the coordinates of a point; and then considerable importance attaches to partial differential coefficients which are formed by a particular rule. Let F(x, y, z) be the function, P a point (x, y, z), P' a neighbouring point (x+Ax, y+Ay, z +Az) , and let AJ be the length of PP'. The value of F(x, y, z) at P may be denoted shortly by F(P). A limit of the same nature as a partial differential coefficient is expressed by the formula =F(P) ' in which As is diminished indefinitely by bringing P' up to P, and P' is supposed to approach P along a straight line, for example, the tangent to a curve or the normal to a surface. The limit in question is denoted by dF/dh, in which it is understood that h indicates a direction, that of PP'. If /, m, n are the direction cosines of the limiting direction of the line PP', supposed drawn from P to P', then dF 9F , dF . 9F The operation of forming dF/dh is called " differentiation with respect to an axis " or " vector differentiation."

41. The most important theorem in regard to partial differential coefficients is the theorem of the total differential. We may write down the equation Total +f(a,b+k)-f(a,b).

Dltteren- If /* is a continuous function of * when x lies between a Ual. and a-\-h and y = b-\-k, and if further/, is a continuous function of y when y lies between b and d-\-k, there exist values of 9 and 17 which lie between o and I and have the properties expressed by the equations f(a+h, b+k)-f(a, b+k)=hf,(a+0h, b+K), f(a, b+k) -f(a, b) = kf u (a, b+yk).

Further, f*(a+eh, b+k) and / (a, b+rik) tend to the limits f, (a, b) and f t (a, b) when h and k tend to zero, provided the differential coefficients/i, / are continuous at the point (a, b). Hence in this case the above equation can be written f(a+h, b+k)-f(a, b) = hf,(a, b)+kf y (a,.b)+R, where TD T)

Hm-A-o, fc-oF =0 and lim 'A-o, fc-oj = - In accordance with the notation of differentials this equation gives Just as in the case of functions of one variable, dx and ay are arbitrary finite differences, and df is not the difference of two values of/, but is so much of this difference as need be retained for the purpose of forming differential coefficients.

The theorem of the total differential is immediately applicable to the differentiation of implicit functions. When y is a function of x which is given by an equation of the form/Gt, y) =o, and it is either impossible or inconvenient to solve this equation so as to express y as an explicit function of x, the differential coefficient dy/dx can be formed without solving the equation. We have at once =_ dx dx I dy' This rule was known, in all essentials, to Fermat and de Sluse before the invention of the algorithm, of the differential calculus.

An important theorem, first proved by Euler, is immediately deducible from the theorem of the total differential. If /(*, y) is a homogeneous function of degree n then The theorem is applicable to functions of any number of variables and is generally known as Euler 's theorem of homogeneous functions. 42. Many problems in which partial differential coefficients occur are simplified by the introduction of certain determinants called " Jacobians " or " functional determinants." They were introduced into Analysis by C. G. J. Jacobi Jacoblaas - (J. f. Math., Crelle, Bd. 22, 1841, p.>3io,). The Jacobian of i, U), . . . Un with respect to xi, xi, . . . x n is the determinant dui dxT dxT dUn dx, du\ in which the constituents of the rth row are the n partial differential coefficients of UT with respect to the n variables x. This determinant is expressed shortly by d\X\, Xtt . . ., Xn)' Jacobians possess many properties analogous to those of ordinary differential coefficients, for example, the following:

9 (MI. 2 u n )^d(x t , d(xi, Xi,..., xj * d(u>, ,...,X n )^d(u^ /N !/ d(yi< yt, . y) d(xi, * 2l . . . flx n ) d(xi, x t *)' If functions ( lt ^, . . . ) of n variables (*i, %,..., x n ) are not independent, but are connected by a relation /(i, ut, . . . ) =o, then d(Ui, U, "")^ - d(Xl, Xl,..., Xn)

and, conversely, when this condition is satisfied identically the functions i, j, . . . , Un are not independent. H3 43. Partial differential coefficients of the second and higher orders can be formed in the same way as those of the first order. For example, when there are two variables x, y, the first partial derivatives df/dx and df/dy are functions of x and ' y, which we may seek to differentiate partially with ' '%" ge ? respect to * or y. The most important theorem in re- alfferenlation to partial differential coefficients of orders higher WaWons than the first is the theorem that the values of such coefficients do not depend upon the order in which the differentiations are performed. For example, we have the equation ~dx (dy) = dy (dx) W This theorem is not true without limitation. The conditions for its validity have been investigated very completely by H. A. Schwarz (see his Ges. math. Abhandlungen, Bd. 2, Berlin, 1890, p.. 275). It is a sufficient, though not a necessary, condition that all the differential coefficients concerned should be continuous functions of x, y. In consequence of the relation (i.) the differential coefficients expressed in the two members of this relation are written 'JPf_ d*f dxdy r dydx' The differential coefficient dx'dydf in which p+q+r = n, is formed by differentiating p times with respect to x, q times with respect to y, r times with respect to z, the differentiations being performed in any order. Abbreviated notations are sometimes used in such forms as Differentials of higher orders are introduced by the defining equation in which the expression (dx-^-+dy^-\ is developed by the binomial <\ theorem in the same way "as if dx^- and dy-^- were numbers, and (d\ r / d \ tl "" r d n f die) \dv) f ' S re P' ace d by Q-pg n-r' When there are more than two variables the multinomial theorem must be used instead of the binomial theorem.

The problem of forming the second and higher differential coefficients of implicit functions can be solved at once by means of partial differential coefficients. For example, if / (x, y)=o is the equation defining y as a function of x, we have \dy I \dy dx* dx dy dxdy\dx d The differential expression Xdx+Ydy, in which both X and Y are functions of the two variables x and y, is a total differential if there exists a function / of x and y which is such that When this is the case we have the relation = dX/dy. (ii.)

Conversely, when this equation is satisfied there exists a function / which is such that df-JUx+Ydy.

The expression Xdx+~Ydy in which X and Y are connected by the relation (ii.) is often described as a " perfect differential." The theory of the perfect differential can be extended to functions of re variables, and in this case there are in(n-i) such relations as (ii.).

In the case of a function of two variables x,y an abbreviated notation is often adopted for differential coefficients. The function being denoted by z, we write . , dz dz d 2 z 3 2 z 5 2 z p, q, r, s, t for , , Partial differential coefficients of the second order are important in geometry as expressing the curvature of surfaces. When a surface is given by an equation of the formz=/(*, y), the lines of curvature are determined by the equation -{(i+p*)s-pqr}(dx)* = o, and the principal radii of curvature are the values of R which satisfy the equation /I/]. The problem of change of variables was first considered by Brook Taylor in his Melhodus incrementorum. In the case conChan of s 'dered by Taylor y is expressed as a function of z, and z . 17 as a function of x, and it is desired to express |the differential coefficients of y with respect to x without eliminating z. The result can be obtained at once by the rules for differentiating a product and a function of a function. We have dy_dy dz_ dx~ dz 'd~x' d*y_dy d?z ~~''' y /<fc\ * \dx)

_ dx 3 ~dz y_ ldz\ z 3 \dx)

The introduction of partial differential coefficients enables us to deal with more general cases of change of variables than that considered above. If u, v are new variables, and x, y are connected with them by equations of the type *=/i , y =/!, (')

while y is either an explicit or an implicit function of x, we have the problem of expressing the differential coefficients of various orders of y with respect to x in terms of the differential coefficients of v with respect to u. We have dy- t s Jin.<*h<*L\ / /2fij.26* dx~ \du^dv du) I \du^~dv du by the rule of the total differential. In the same way, by means of differentials of higher orders, we may express d?yldx*, and so on.

Equations such as (i.) may be interpreted as effecting a transformation by which a point (u, v) is made to correspond to a point (x, y). The whole theory of transformations, and of functions, or differential expressions, which remain invariant under groups of transformations, has been studied exhaustively by Sophus Lie (see, in particular, his Theorie der Transformationsgruppen, Leipzig, 1888-1893). (See also DIFFERENTIAL EQUATIONS and GROUPS).

A more general problem of change of variables is presented when it is desired to express the partial differential coefficients of a function V with respect to x, y, . . . in terms of those with respect to u, v, . . ., where u, v, . . . are connected with x, y, . . . by any functional relations. When there are two variables x, y, and u, v are given functions of #, y, we have dV_ = dVdu, svdv_, dx du dx~^dv dx _ = _ dy du dy dv dy and the differential coefficients of higher orders are to be formed by repeated applications of the rule for differentiating a product and the rules of the type d__du__d_, d _d__ dx~dxdu dxdv When *, y are given functions of ,,... we have, instead of the above, such equations as du dx du dy du' and dV/dx, dV/dy can be found by solving these equations, provided the Jacobian d(x,y)/d(u,v) is not zero. The generalization of this method for the case of more than two variables need not detain us.

In cases like that here considered it is sometimes more convenient not to regard the equations connecting *, y with u, v as effecting a point transformation, but to consider the loci w = const., n=const. as two " families " of curves. Then in any region of the plane of (x, y) in which the Jacobian d(x, y)jd(u, v) does not vanish or become infinite, any point (x, y) is uniquely determined by the values of and v which belong to the curves of the two families that pass through the point. Such variables as u,v are then described as "(curvilinear coordinates " of the point. This method is applicable to any number of variables. When the loci = const., . . . intersect each other at right angles, the variables are " orthogonal " curvilinear coordinates. Three-dimensional systems of such coordinates have important applications in mathematical physics. Reference may be made to G. Lam6, Lemons sur les coordonnees curvilignes (Paris, 1859), and to G. Darboux, Legons sur les coordonnees curvilignes et sys&mes orthogonaux (Paris, 1898).

When such a coordinate as is connected with x and y by a functional relation of the form f(x,y,u)=o the curves = const. are a family of curves, and this family may be such that no two curves of the family have a common point. When this is not the case the points in which a curve f(x,y,u)=o is intersected by a curve /(x, y, u+&u) =o tend to limiting positions as AM is diminished indefinitely. The locus of these limiting positions is the " envelope " of the family, and in general it touches all the curves of the family. It is easy to see that, if u, v are the parameters of two families of curves which have envelopes, the Jacobian d(x,y)/d(u,v) vanishes at all points on these envelopes. It is easy to see also that at any point where the reciprocal Jacobian d(u,v)/d(x,y) vanishes, a curve of the family u touches a curve of the family v.

If three variables x, y,z are connected by a functional relation f(x, y,z)=o, one of them, z say, may be regarded as an implicit function of the other two, and the partial differential coefficients of z with respect to x and y can be formed by the rule of the total differential. We have =_ .=- dx dx/ dz' dy dyl dz' and there is no difficulty in proceeding to express the higher differential coefficients. There arises the problem of expressing the partial differential coefficients of x with respect to y and z in terms of those of z with respect to x and y. The problem is known as that of " changing the dependent variable." It is solved by applying the rule of the total differential. Similar considerations are applicable to all cases in which n variables are connected by fewer than n equations.

45. Taylor's theorem can be extended to functions of several variables. In the case of two variables the general for- Extension mula, with a remainder after n terms, can be written O f Taylor's most simply in the form theorem.

f(a+h, 6+fe) =/(a, b)+df(a, in which and finff _L/)/. h-i-fib~\ I I Ji l fr L\ dx' dyi J x-a+gh, v -b+ft The last expression is the remainder after n terms, and in it 9 denotes some particular number between o and I. The results for three or more variables can be written in the same form. The extension of Taylor's theorem was given by Lagrange (1797); the form written above is due to Cauchy (1823). For the validity of the theorem in this form it is necessary that all the differential coefficients up to the nth should be continuous in a region bounded by x = a=h, y = bk. When all the differential coefficients, no matter how high the order, are continuous in such a region, the theorem leads to an expansion of the function in a multiple power series. Such expansions are just as important in analysis, geometry and mechanics as expansions of functions of one variable. Among the problems which are solved by means of such expansions are the problem of maxima and minima for functions of more than one variable (see MAXIMA and MINIMA).

46. In treatises on the differential calculus much space is usually Plane devoted to the differential geometry of curves and curves. surfaces. A few remarks and results relating to the differential geometry of plane curves are set down here.

(i.) If ^ denotes the angle which the radius vector drawn from the origin makes with the tangent to a curve at a point whose polar coordinates are r, 8 and if p denotes the perpendicular from the origin to the tangent, then cos ^ = dr/ds, sin ^ = rd8/ds=p/r, ' where ds denotes the element of arc. The curve may be determined by an equation connecting p with r.

(ii.) The locus of the foot of the perpendicular let fall from the origin upon the tangent to a curve at a point is called the pedal of the curve with respect to the origin. The angle \f/ for the pedal is the same as the angle \f/ for the curve. Hence the (p, r) equation of the pedal can be deduced. If the pedal is regarded as the primary curve, the curve of which it is the pedal is the " negative pedal " of the primary. We may have pedals of pedals and so on, also negative pedals of negative pedals and so on. Negative pedals are usually determined as envelopes.

(iii.) If </> denotes the angle which the tangent at any point makes with a fixed line, we have (iv.) The " average curvature " of the arc As of a curve between two points is measured by the quotient where the upright lines denote, as usual, that the absolute value of the included expression is to be taken, and <t> is the angle which the tangent makes with a fixed line, so that A< is the angle between the tangents (or normals) at the points. As one of the points moves up to coincidence with the other this average curvature tends to a limit which is the " curvature " of the curve at the point. It is denoted by 1 4m |sr| Sometimes the upright lines are omitted and a rule of signs is given : Let the arc s of the curve be measured from some point along the curve in a chosen sense, and let the normal be drawn towards that side to which the curve is concave; if the normal is directed towards the left of an observer looking along the tangent in the chosen sense of description the curvature is reckoned positive, in the contrary case negative. The differential d<t> is often called the " angle of contingence." In the 14th century the size of the angle between a curve and its tangent seems to have been seriously debated, and the name " angle of contingence " was then given to the supposed angle.

(v.) The curvature of a curve at a point is the same as that of a certain circle which touches the curve at the point, and the "radius of curvature" p is the radius of this circle. We have - = p \ds The centre of the circle is called the " centre of curvature "; it is the limiting position of the point of intersection of the normal at the point and the normal at a neighbouring point, when the second point moves up to coincidence with the first. If a circle is described to intersect the curve at the point P and at two other points, and one of these two points is moved up to coincidence with P, the circle touches the curve at the point P and meets it in another point ; the centre of the circle is then on the normal. As the third point now moves up to coincidence with P, the centre of the circle moves to the centre of curvature. The circle is then said to " osculate " the curve, or to have_ " contact of the second order " with it at P.

(vi.) The following are formulae for the radius of curvature :

a v "') T e P oin ' s at which the curvature vanishes are " points of inflection. If P is a point of inflection and Q a neighbouring point, then, as Q moves up to coincidence with P, the distance from P to the point of intersection of the normals at P and Q becomes greater than any distance that can be assigned. The equation which gives the abscissae of the points in which a straight line meets the curve being expressed in the form f(x) =o, the function f(x) has a factor (x xo) 3 , where xt> is the abscissa of the point of inflection P, and the line is the tangent at P. When the factor (x Xo) occurs (n + I ) times in /(x), the curve is said to have " contact of the nth order " with the line. There is an obvious modification when the line is parallel to the axis of y.

(viii.) The locus of the centres of curvature, or envelope of the normals, of a curve is called the " evolute." A curve which has a given curve as evolute is called an " involute " of the given curve. All the involutes are " parallel " curves, that is to say, they are such that one is derived from another by marking off a constant distance along the normal. The involutes are " orthogonal trajectories " of the tangents to the common evolute.

(ix.) The equation of an algebraic curve of the nth degree can be expressed in the form WO+KI + WJ+ . . . +*= o, where u<> is a constant, and u r is a. homogeneous rational integral function of x, y of the rth degree. When the origin is on the curve, vanishes, and tti=o represents the tangent at the origin. If HI also vanishes, the origin is a double point and u% = o represents the tangents at the origin. If 2 has distinct factors, or is of the form a(y pix)(y pix), the value of y on either branch of the curve can be expressed (for points sufficiently near the origin) in a power series, which is either . . . , or where qi, . . . and qi, . . . are determined without ambiguity. If pi and pi are real the two branches have radii of curvature pi, pj determined by the formulae -<+*> When pi and /> 2 are imaginary the origin is the real point of intersection of two imaginary branches. In the real figure of the curve it is an isolated point. If 2 is a square, a(y px)*, the origin is a cusp, and in general there is not a series for y in integral powers of x, which is valid in the neighbourhood of the origin. The further investigation of cusps and multiple points belongs rather to analytical geometry and the theory of algebraic functions than to differential calculus.

(x.) When the equation of a curve is given in the form Wo+i + . . . +Hn-i+ttn = o where the notation is the same as that in (ix.), the factors of n determine the directions of the asymptotes. If these factors are all real and distinct, there is an asymptote corresponding to each factor. If u.^LiLj . . . L n , where Li, . . . are linear in x, y, we may resolve Mn_i/ into partial fractions according to the formula n-l Ai . A 2 , and then Li+Ai=o, L 2 +A2 = o, . . . are the equations of the asymptotes. When a real factor of is repeated we may have two parallel asymptotes or we may have a " parabolic asymptote." Sometimes the parallel asymptotes coincide, as in the curve x*(x'+y 2 a?')=a 4 , where x = o is the only real asymptote. The whole theory of asymptotes belongs properly to analytical geometry and the theory of algebraic functions.

47. The formal definition of an integral, the theorem of the existence of the integral for certain classes of functions, a list of classes of " integrable " functions, extensions of the notion integral of integration to functions which become infinite or in- calculus determinate, and to cases in which the limits of integration become infinite, the definitions of multiple integrals, and the possibility of defining functions by means of definite integrals all these matters have been considered in FUNCTION. The definition of integration has been explained in 5 above, and the results of some of the simplest integrations have been given in 12. A few theorems relating to integrations have been noted in 34, 35, 36 above.

48. The chief methods for the evaluation of indefinite methods of integrals are the method of integration by parts, and the / n<e rar / on introduction of new variables.

From the equation d(uv) =udv+vdu we deduce the equation dv , C du , dx=uv-J Vfadx, ) u dx a or, as it may be written uwdx = This is the rule of " integration by parts." As an example we have 1 1/-*! /*-?-/?<<* When we introduce a new variable z in place of x, by means of an equation giving x in terms of z, we express f(x) in terms of z. Let <t>(z) denote the function of z into which f(x) is transformed. Then from the equation j dx , we deduce the equation ) gf As an example, in the integral /V (i -x 2 )<2x put x = sin z; the integral becomes Integration la terms of elementary f unctions.

49. The indefinite integrals of certain classes of functions can be expressed by means of a finite number of operations of addition or multiplication in terms of the so-called " elementary " functions. The elementary functions are rational algebraic functions, implicit algebraic functions, exponentials and logarithms, trigonometrical and inverse circular functions. The following are among the classes of functions whose integrals involve the elementary functions only: (i.) all rational functions; (ii.) all irrational functions of the form f(x, y), where / denotes a rational algebraic function of x and y, and y is connected with x by an algebraic equation of the second degree; (iii.) all rational functions of sin x and cos x; (iv.) all rational functions of e?; (v.) all rational integral functions of the variables x, e" x , & x , ... sin mx, cos mx, sin nx, cos nx, ... in which a, 6, ... and m, n, . . . are any constants. The integration of a rational function is generally effected by resolving the function into partial fractions, the function being first expressed as the quotient of two rational integral functions. Corresponding to any simple root of the denominator there is a logarithmic term in the integral. If any of the roots of the denominator are repeated there are rational algebraic terms in the integral. The operation of resolving a fraction into partial fractions requires a knowledge of the roots of the denominator, but the algebraic part of the integral can always be found without obtaining all the roots of the denominator. Reference may be made to C. Hermite, Cours d' analyse, Paris, 1873. The integration of other functions, which can be integrated in terms of the elementary functions, can usually be effected by transforming the functions into rational functions, (possibly after preliminary integrations by parts. In the case of rational functions of x and a radical of the form V (ax*+bx+c) the radical can be reduced by a linear substitution tojone of the forms V (a 1 *>), V (x 2 a 2 ), V (x 2 +a 2 ). The substitutions x = a sin 8, x = a sec 8, x = a tan are then effective in the three cases. By these substitutions the subject of integration becomes a rational function of sin 8 and cos 0, and it can be reduced to a rational function of t by the substitution tan \S = t. There are many other substitutions by which such integrals can be determined. Sometimes we may have information as to the functional character of the integral without being able to determine it. For example, when the subject of integration is of theform (ax t +bx 3 +cx 2 +dx+e)-l the integral cannot be expressed explicitly in terms of elementary functions. Such integrals lead to new functions (see FUNCTION).

Methods of reduction and substitution for the evaluation of indefinite integrals occupy a considerable space in text-books of the integral calculus. In regard to the functional character of the integral reference may be made to G. H. Hardy's tract, The Integration of Functions of a Single Variable (Cambridge, 1905), and to the memoirs there quoted. A few results are added here (i.) f(x t +a)-Ux='\og(x+(if+aW.

Ca " be evaluated ^ the substitution x-p = ilz, andj (x _^ v( g. +ate+e) can be deduced by differentiating (re i) times with respect to p.

(Hi ) f , ( "*+, I P d f, . , . can be reduced by the sub- J (ox 2 +2/3x+7)V (ax 2 +20X+c)

stitution y* = (ax*+2bx+c)!(ax*+2px+y) to the form A /V(xf-^) +B /V(/-x 2 )

where A and B are constants, and Xi and X 2 are the two values of X for which (a-Xa)x 2 +2(6-X/3)x+c-XT is a perfect square (see A. G. Greenhill, A Chapter in the Integral Calculus, London, 1888). (iv.) (x m (ax n +V)fdx, in which m, n, p are rational, can be reduced, by putting ax" = bt, to depend upon /(i +t)"dt. If p is an integer and q a fraction r/s, we put t = *. If q is an integer and p = r/s we put i+t = u'. If p+q is an integer and p = r/s we put l+t = tu'. These integrals, called " binomial integrals," were investigated by Newton (De quadratura cuniarum).

(vii.) fe al s (viii.) / sin" 1 x cos" x dx can be reduced by differentiating a function of the form sin p x cos' x; Hence d sin x _ e '&' dx cos 9 x '' dx cos" x (n I ) cos (ix.)

l l-IXnaninteger).

/'AlT /*!JT (x.) I sin z +1 x</x= I cos 2 " +1 x(fx = 4---2 i>(nanintcger).

'o Jo 3 -5- \2n-r i)

(xi.) J ; ; . ecosx )n can be reduced by one of the substitutions _ e+cosx , _ e+cosx i+cosx' i+ecosx' of which the first or the second is to be employed according as e<or>i.

50. Among the integrals of transcendental functions New transwhich lead to new transcendental functions we may notice ccndeats.

~" -' x t <Ldz, called the " logarithmic integral," and denoted by " Li x," also the integrals called the " sine integral " and the " cosine integral," and denoted by " Si x " and " Ci x, also the integral called the " error-function integral," and denoted by " Erf x." All these functions have been tabulated (seeTABLES, MATHEMATICAL). 51. New functions can be introduced also by means of the definite integrals of functions of two or more variables with re- Eulerlaa spect to one of the variables, the limits of integration integrals, being fixed. Prominent among such functions are the Beta and Gamma functions expressed by the equations B(/, ro) = When n is a positive integer r(re + i)=n ! . The Beta function (or " Eulerian integral of the first kind ") is expressible in terms of Gamma functions (or " Eulerian integrals of the second kind ") by the formula The Gamma function satisfies the difference equation r(* + l)=xr(x), and also the equation r(x).r(l -*)=*/ sin (XT), with the particular result The number is called " Euler's constant," and is equal to the limit lim. [(i + l + l+...+i) -logn]; its value to 15 decimal places is 0^577 215 664 901 532. The function log T(i +x) can be expanded in the series log r(i+*)-J i -J log where and the series for log r(i+x) converges when x lies between I and i.

52. Definite integrals can sometimes be evaluated when the limits of integration are some particular numbers, although Dellalte the corresponding indefinite integrals cannot be found. For example, we have the result P(i-* 2 )-i \ogxdx= -i?rlog2, Jo although the indefinite integral of (i 3C 2 )~J log x cannot be found. Numbers of definite integrals are expressible in terms of the transcendental functions mentioned in 50 or in terms of Gamma functions. For the calculation of definite integrals we have the following methods :

(i.) Differentiation with respect to a parameter. (ii.) Integration with respect to a parameter.

(iii.) Expansion in infinite series and integration term by term.

(iv.) Contour integration.

The first three methods involve an interchange of the order of two limiting operations, and they are valid only when the functions satisfy certain conditions of continuity, or, in case the limits of integration are infinite, when the functions tend to zero at infinite distances in a sufficiently high order (see FUNCTION). The method of contour integration involves the introduction of complex variables (see FUNCTION : Complex Variables). A few results are added (ii.)j ?- 3^ dx=ir(cotair cotbir), (o<oor6<i), (iv.) *".cos zx.er^dx = - JC-'VTT, (vii.) I log(i 2 J o d* = oor2irlogaaccordingasa<or>i, y o cos ax-cos ftx / \ r (xn.)

./ (xiii.)f" ./ /"" (xiv.) I .Jo | j Jo 53. The meaning of integration of a function of n variables through a domain of the same number of dimensions is explained in the article FUNCTION. In the case of two variables x, y we Multiple integrate a function f(x,y) over an area; in the case of Integrals. tnree variables x, y, z we integrate a function f(x, y, z) through a volume. The integral of a f unction f(x, y) over an area in the plane of (x, y) is denoted by ///(*, y)dxdy.

The notation refers to a method of evaluating the integral. We may suppose the area divided into a very large number of very small rectangles by lines parallel to the axes. Then we multiply the value of / at any point within a rectangle by the measure of the area of the rectangle, sum for all the rectangles, and pass to a limit by increasing the number of rectangles indefinitely and diminishing all their sides indefinitely. The process is usually effected by summing first for all the rectangles which lie in a strip between two lines parallel to one axis, say the axis of y, and afterwards for all the strips. This process is equivalent to integrating /(*, y) with respect to y, keeping x constant, and taking certain functions of x as the limits of integration for y, and then integrating the result with respect to x between constant limits. The integral obtained in this way may be written in such a form as and is called a " repeated integral." The identification of a surface integral, such as Jff(x, y)dxdy, with a repeated integral cannot always be made, but implies that the function satisfies certain conditions of continuity. In thej same way volume integrals are usually evaluated by regarding them as repeated integrals, and a volume integral is written in the form S!JS(x,y,z}dxdydz.

Integrals such as surface and volume integrals are usually called " multiple integrals." Thus we have " double " integrals, " triple " integrals, and so on. In contradistinction to multiple integrals the ordinary integral of a function of one variable with respect to that variable is called a " simple integral.

A more general type of surface integral may be defined by taking an arbitrary surface, with or without an edge. We suppose in the Surface place that the surface is closed, or has no edge. We Integrals.

may mark a large number of points on the surface, and draw the tangent planes at all these points. These tangent planes form a polyhedron having a large number of faces, one to each marked point ; and we may choose the marked points so that all the linear dimensions of any face are less than some arbitrarily chosen length. We may devise a rule for increasing the number of marked points indefinitely and decreasing the lengths of all the edges of the polyhedra indefinitely. If the sum of the areas of the faces tends to a limit, this limit is the area of the surface. If we multiply the value of a function /at a point of the surface by the measure of the area of the corresponding face of thejpolyhedron, sum for all the faces, and pass to a limit as before, the result is a surface integral, and is written ///as- The extension to the case of an open surface bounded by an edge presents no difficulty. A line integral taken "' along a curve is defined in a similar way, and is written ffds where ds is the element of arc of the curve ( 33). The direction cosines of the tangent of a curve are dx/ds, dy/ds, dz/ds, and line integrals usually present themselves in the form "tegrais.

S s .

In like manner surface integrals usually present themselves in the where /, m, n are the direction cosines of the normal to the surface drawn in a specified sense.

The area of a bounded portion of the plane of (x, y) may be expressed either as V(xdy-ydx), or as JJdxdy, the former integral being a line integral taken round the boundary of the portion, and the latter a surface integral taken over the area within this boundary. In forming the line integral the boundary is supposed to be described in the positive sense, so that the included area is on the left hand.

530. We have two theorems of transformation connect- Theorems ing volume integrals with surface integrals and surface of Green integrals with line integrals. The first theorem, called and " Green's theorem," is expressed by the equation Stokes.

JIf (8+i? + S where the volume integral on the'left is taken through the volume within a closed surface S, and the surface integral on the right is taken over S, and I, m, n denote the direction cosines of the normal to S drawn outwards. There is a corresponding theorem for a closed curve in two dimensions, viz., the sense of description of s being the positive sense. This theorem is a particular case of a more general theorem called " Stokes's theorem." Let s denote the edge of an open surface S, and let S be covered with a network of curves so that the meshes of the network are nearly plane', then we can choose a sense of description of the edge of any mesh, and a corresponding sense for the normal to S at any point within the mesh, so that these senses are related like the directions of rotation and translation in a right-handed screw. This convention fixes the sense of the normal (I, m, n) at any point on S when the sense of description of s is chosen. If the axes of x, y, z are a right-handed system, we have Stokes's theorem in the form where the integral on the left is taken round the curve j in the chosen sense. When the axes are left-handed, we may either reverse the sense of /, m, n and maintain the formula, or retain the sense of /, TO, n and change the sign of the right-hand member of the equation. For the validity of the theorems of Green and Stokes it is in general necessary that the functions involved should satisfy certain conditions of continuity. For example, in Green's theorem the differential coefficients d/dx, dy/dy, df/dz must be continuous within S. Further, there are restrictions upon the nature of the curves or surfaces involved. For example, Green's theorem, as here stated, applies only to simply-connected regions of space. The correction for multiply-connected regions is important in several physical theories.

54. The process of changing the variables in a multiple integral, such as a surface or volume integral, is divisible into two stages. It is necessary in the first place to determine the differential element expressed by the product of the differentials of the first set of variables in terms of the differentials of the second set of variables. It is necessary in the second place to determine the limits of integration which must be employed when the integral in terms of the new variables is evaluated as a repeated integral. The first part of the problem is solved at once by the introduction of the Jacobian. If the variables of one set are denoted by x it xt ..... *, and those of the other set by MI, 2 , . . ., n , we have the relation r, ."*". ! * " In regard to the second stage of the process the limits of integration must be determined by the rule that the integration with respect to the second set of variables is to be taken through the same domain as the integration with respect to the first set.

For example, when we have to integrate a f unction f(x, y) over the area within a circle given by x i +y' = a i , and we introduce polar coordinates so that x = r cos 8, y = r sin 8, we find that r is the value of the Jacobian, and that all points within or on the circle are given by a^ r^o, 2ir>0^ o, and we have If we have to integrate over the area of a rectangle a^x^o, 6^ y^o, and we transform to polar coordinates, the integral becomes the sum of two integrals, as follows:

:an-16/a , fa sec 9..

f(rcos\$, rsmO)rdr J t f(rcos8,rsin8)rdr. tan-14/a J 55. A few additional results in relation to line integrals and multiple integrals are set down here.

(i.) Any simple integral can be regarded as a line-integral taken . to along a portion of the axis of x. When a change of . . . variables is made, the limits of integration with respect to the new variable must be such that the domain of '" . . integration is the same as before. This condition may munipie require the replacing of the original integral by the sum ' of two or more simple integrals.

(ii.) The line integral of a perfect differential of a one-valued function, taken along any closed curve, is zero.

(iii.) The area within any plane closed curve can be expressed by either of the formulae flr*dO orfipds, where r, 8 are polar coordinates, and p is the perpendicular drawn from a fixed point to the tangent. The integrals are to be understood as line integrals taken along the curve. When the same integrals are taken between limits which correspond to two points of the curve, in the sense of line integrals along the arc between the points, they represent the area bounded by the arc and the terminal radii vectores.

(iv.) The volume enclosed by a surface which is generated by the revolution of a curve about the axis of x is expressed by the formula and the area of the surface is expressed by the formula 2vfyds, where ds is the differential element of arc of the curve. When the former integral is taken between assigned limits it represents the volume contained between the surface and two planes which cut the axis of x at right angles. The latter integral is to be understood as a line integral taken along the curve, and it represents the area of the portion of the curved surface which is contained between two planes at right angles to the axis of x.

(v.) When we use curvilinear coordinates {, i\ which are conjugate functions of x, y, that is to say are such that the Jacobian d(Gr., -n)ld(x, y) can be expressed in the form ,3*.

and in a number of equivalent forms. The area of any portion of the plane is represented by the double integral where J denotes the above Jacobian, and the integration is taken through a suitable domain. When the boundary consists of portions of curves for which = const., or ij = const., the above is generally the simplest way of evaluating it.

(vi.) The problem of " rectifying " a plane curve, or finding its length, is solved by evaluating the integral or, in polar coordinates, by evaluating the integral In both cases the integrals are line integrals taken along the curve.

(vii.) When we use curvilinear coordinates , ij as in (v.) above, the length of any portion of a curve {= const, is given by the integral taken between appropriate limits for r/. There is a similar formula for the arc of a curve 17 = const.

(viii.) The area of a surface z=f(x, y) can be expressed by the formula When the coordinates of the points of a surface are expressed as functions of two parameters u, v, the area is expressed by the formula f f\~ ( d(y, z) ) 2 ( d(z, x) ) 2 , ( d(x, y) ) t ~\ t , J J L / d(j< t>) I ' J d(u v) \ ' I d (u v) \ dudv. When the surface is referred to three-dimensional polar coordinates r, 6, \$ given by the equations x=r sin cos <, y = r sin0 sm<t>, z = r cos0, and the equation of the surface is of the form r=/(0,<#>), the area is expressed by the formula The surface integral of a function of (8, <j>) over the surface of a Sphere r = const, can be expressed in the form In every case the domain of integration must be chosen so as to include the whole surface.

(ix.) In three-dimensional polar coordinates the Jacobian d(x,y, d(r, 8, < The volume integral of a function F (r, 6, <p) through the volume of a sphere r = a is = r 2 sin9.

(x.) Integrations of rational functions through the volume of an ellipsoid x f /a?+y 1 /b'*+z'/c'* = i are often effected by means of a general theorem due to Lejeune Dirichlet (1839), which is as follows: when the domain of integration is that given by the inequality where the a's and o's are positive, the value of the integral //. . . *!!-!. X2-2- 1 . . .

If, however, the object aimed at is an integration through the volume of an ellipsoid it is simpler to reduce the domain of integration to that within a sphere of radius unity by the transformation x = a, y = OTt, 2 = cf, and then to perform the integration through the Sphere by transforming to polar coordinates as in (ix).

56. Methods of approximate integration began to be devised very early. Kepler's practical measurement of the focal sectors A of ellipses (1609) was an approximate integration, as also ' was the method for the quadrature of the hyperbola given by James Gregory in the appendix to his Exercitationes ' geometricae (1668). In Newton's Methodus differential ' (i 7 1 1 ) the subject was taken up systematically. Newton's te * ral object was to effect the approximate quadrature of a given curve by making a curve of the type pass through the vertices of (+i) equidistant ordinates of the given curve, and by taking the area of the new curve so determined as an approximation to the area of the given curve. In 1743 Thomas Simpson in his Mathematical Dissertations published a very convenient rule, obtained by taking the vertices of three consecutive equidistant ordinates to be points on the same parabola. The distance between the extreme ordinates corresponding to the abscissae x = a and* = 6 is divided into2n equal segments by ordinates yi,yi, . . .ytn-i, and the extreme ordinates are denoted by yo, y2n. The vertices of the ordinates yo, yi, y2 lie on a parabola with its axis parallel to the axis of y, so do the vertices of the ordinates y 2 , y 3 , yt, and so on. The area is expressed approximately by the formula which is known as Simpson's rule. Since all simple integrals can be represented as areas such rules are applicable to approximate integration in general. For the recent developments reference may be made to the article by A. Voss in Ency. d. Math. Wiss., Bd. II., A. 2 (1899), and to a monograph by B. P. Moors, Valeur approximative d'une integrate definie (Paris, 1905).

Many instruments have been devised for registering mechanically the areas of closed curves and the values of integrals. The best known are perhaps the " planimeter " of J. Amsler (1854) and the " integraph " of Abdank-Abakanowicz (1882).

BIBLIOGRAPHY. For historical questions relating to the subject th chief authority is M. Cantor, Geschichte d. Mathematik (3 Bde., Leipzig, 1894-1901). For particular matters, or special periods, the following may be mentioned: H. G. Zeuthen, Geschichte d. Math, im Altertum u. Mittelalter (Copenhagen, 1896) and Gesch. d. Math, im XVI. u. XVII. Jahrhundert (Leipzig, 1903); S. Horsley, Isaaci Newtoni opera quae exstant omnia (5 vols., London, 1779-1785); C. I. Gerhardt, Leibnizens math. Schriften (7 Bde., Leipzig, 1849- 1863); Joh. Bernoulli, Opera omnia (4 Bde., Lausanne and Geneva, 1742). Other writings of importance in the history of the subject are cited in the course of the article. A list of some of the more important treatises on the differential and integral calculus is appended. The list has no pretensions to completeness ; in particular, most of the recent books in which the subject is presented in an elementary way for beginners or engineers are omitted. L. Euler, Institutions! calculi differentialis (Petrop., 1755) and Institutiones calculi integralis (3 Bde., Petrop., 1768-1770); J. L. Lagrange, Lemons sur le calcul des fonctions (Paris, 1806, (Euvres, t. x.), and Theorie des fonctions analytiques (Paris, 1797, 2nded., 1813, (Euvres, t. ix.) ; S. F. Lacroix, Traite de calcul diff. et de calcul int. _(3 tt., Paris, 1808-1819). There have been numerous later editions; a translation by Herschel, Peacock and Babbage of an abbreviated edition of Lacroix's treatise was published at Cambridge in 1816. G. Peacock, Examples of the Differential and Integral Calculus (Cambridge, 1820); A. L. Cauchy, Resume des lemons . . . sur le calcul infinitesimals (Paris, 1823), and Lemons sur le calcul differential (Paris, 1829; (Euvres, ser. 2, t. iv.) ; F. Minding, Handbuchd. Diff.-u. Int.-Rechnung (Berlin, 1836) ; F. Moigno, Legons sur le calcul diff. (4 tt., Paris, 1840-1861) ; A. de Morgan, Diff. and Int. Calc. (London, 1842); D. Gregory, Examples on the Diff. and Int. Calc. (2 vols., Cambridge, 1841-1846) ; I. Todhunter, Treatise on the Diff. Calc. and Treatise on the Int. Calc. (London, 1852), numerous later editions ; B. Price, Treatise on the Infinitesimal Calculus (2 vols., Oxford, 1854), numerous later editions; D. Bierens de Haan, Tables d'integrales definies (Amsterdam, 1858) ; M. Stegemann, Grundriss d. Diff.- u. Int.-Rechnung (2 Bde., Hanover, 1862) numerous later editions; J. Bertrand, Traite de calc. diff. et int. (2 tt., Paris, 1864-1870); J. A. Serret, Cours de calc. diff. et int. (2 tt., Paris, 1868, 2nd ed., 1880, German edition by Harnack, Leipzig, 1884-1886, later German editions by Bohlmann, 1896, and Scheffers,| 1906,! incomplete); B.Williamson, Treatise on the Diff. Calc. (Dublin, 1872), and Treatise on the Int. Calc. (Dublin, 1874) numerous later editions of both ; also the article " Infinitesimal Calculus " in the 9th ed. of the Ency. Brit.; C. Hermite, Cours d' analyse (Paris, 1873); O. Schlomilch, Compendium d. hoheren Analysis (2 Bde., Leipzig, 1874) numerous later editions; J. Thomae, Einleitung in d. Theorie d. bestimmten Integrate (Halle, 1875); R. Lipschitz, Lehrbuch d. Analysis (2 Bde., Bonn, 1877, 1880); A. Harnack, Elemente d. Diff.-u. Int.-Rechnung (Leipzig, 1882, Eng. trans, by Cathcart, London, 1891); M. Pasch, Einleitung in d. Diff.-u. Int.-Rechnung (Leipzig, 1882) ; Genocchi and Peano, Calcolo differencials (Turin, 1884, German edition by Bohlmann and Schepp, Leipzig, 1898, 1899); H. Laurent, Traite d'analyse (7 tt., Paris, 1885-1891); J. Edwards, Elementary Treatise on the Diff. Calc. (London, 1886), several later editions; A. G. Greenhill, Diff. and Int. Calc. (London, 1886, 2nd ed., 1891); E. Picard, Traite d'analyse (3 tt., Paris, 1891-1896); O. Stolz, Grundziige d. Diff.- u. Int.-Rechnung (3 Bde., Leipzig, 1893-1899); C. Jordan, Cours d'analyse (3 tt., Paris, 1893-1896); L. Kronecker, Vorlesungen u. d. Theorie d. einfachen u. vielfachen Integrale (Leipzig, 1894); J. Perry, The Calculus for Engineers (London, 1897); H. Lamb, An Elementary Course of Infinitesimal Calculus (Cambridge, 1897) ; G. A. Gibson, An Elementary Treatise on the Calculus (London, 1901) ; E. Goursat, Cours d'analyse mathematique (2 tt., Paris, 1902- 1905) ; C.-J. de la Vallee Poussin, Cours d'analyse infinMsimale (2 tt., Louvainand Paris, 1903-1906); A. E. H. Love, Elements of the Diff. and Int. Calc. (Cambridge, 1909) ; W. H. Young, The Fundamental Theorems of the Diff. Calc. (Cambridge, 1910). A r6sum(5 of the infinitesimal calculus is given in the articles " Diff.-u. Int- Rechnung " by A. Voss, and " Bestimmte Integrale " by G. Brunei in Ency. d. math. Wiss. (Bde. ii. A. 2, and ii. A. 3, Leipzig, 1899, 1900). Many questions of principle are discussed exhaustively by E. W. Hobson, The Theory of Functions of a Real Variable (Cambridge, 1907). (A. E. H. L.) Note - this article incorporates content from Encyclopaedia Britannica, Eleventh Edition, (1910-1911)  