Cavalieri's quadrature formula

From HandWiki
Short description: Mathematical term in calculus
Cavalieri's quadrature formula computes the area under the cubic curve, together with other higher powers.

In calculus, Cavalieri's quadrature formula, named for 17th-century Italian mathematician Bonaventura Cavalieri, is the integral

[math]\displaystyle{ \int_0^a x^n\,dx = \tfrac{1}{n+1}\, a^{n+1} \qquad n \geq 0, }[/math]

and generalizations thereof. This is the definite integral form; the indefinite integral form is:

[math]\displaystyle{ \int x^n\,dx = \tfrac{1}{n+1}\, x^{n+1} + C \qquad n \neq -1. }[/math]

There are additional forms, listed below. Together with the linearity of the integral, this formula allows one to compute the integrals of all polynomials.

The term "quadrature" is a traditional term for area; the integral is geometrically interpreted as the area under the curve y = xn. Traditionally important cases are y = x2, the quadrature of the parabola, known in antiquity, and y = 1/x, the quadrature of the hyperbola, whose value is a logarithm.

Forms

Negative n

For negative values of n (negative powers of x), there is a singularity at x = 0, and thus the definite integral is based at 1, rather than 0, yielding:

[math]\displaystyle{ \int_1^a x^n\,dx = \frac{1}{n+1} (a^{n+1} - 1) \qquad n \neq -1. }[/math]

Further, for negative fractional (non-integer) values of n, the power xn is not well-defined, hence the indefinite integral is only defined for positive x. However, for n a negative integer the power xn is defined for all non-zero x, and the indefinite integrals and definite integrals are defined, and can be computed via a symmetry argument, replacing x by −x, and basing the negative definite integral at −1.

Over the complex numbers the definite integral (for negative values of n and x) can be defined via contour integration, but then depends on choice of path, specifically winding number – the geometric issue is that the function defines a covering space with a singularity at 0.

n = −1

There is also the exceptional case n = −1, yielding a logarithm instead of a power of x:

[math]\displaystyle{ \int_1^a \frac{1}{x}\,dx = \ln a, }[/math]
[math]\displaystyle{ \int \frac{1}{x}\,dx = \ln x + C, \qquad x \gt 0 }[/math]

(where "ln" means the natural logarithm, i.e. the logarithm to the base e = 2.71828...).

The improper integral is often extended to negative values of x via the conventional choice:

[math]\displaystyle{ \int \frac{1}{x}\,dx = \ln |x| + C, \qquad x \neq 0. }[/math]

Note the use of the absolute value in the indefinite integral; this is to provide a unified form for the integral, and means that the integral of this odd function is an even function, though the logarithm is only defined for positive inputs, and in fact, different constant values of C can be chosen on either side of 0, since these do not change the derivative. The more general form is thus:[1]

[math]\displaystyle{ \int\frac{1}{x} \, dx= \begin{cases} \ln |x| + C^- & x \lt 0 \\ \ln |x| + C^+ & x \gt 0 \end{cases} }[/math]

Over the complex numbers there is not a global antiderivative for 1/x, due this function defining a non-trivial covering space; this form is special to the real numbers.

Note that the definite integral starting from 1 is not defined for negative values of a, since it passes through a singularity, though since 1/x is an odd function, one can base the definite integral for negative powers at −1. If one is willing to use improper integrals and compute the Cauchy principal value, one obtains [math]\displaystyle{ \int_{-c}^c \frac{1}{x}\,dx = 0, }[/math] which can also be argued by symmetry (since the logarithm is odd), so [math]\displaystyle{ \int_{-1}^1 \frac{1}{x}\,dx = 0, }[/math] so it makes no difference if the definite integral is based at 1 or −1. As with the indefinite integral, this is special to the real numbers, and does not extend over the complex numbers.

Alternative forms

The integral can also be written with indexes shifted, which simplify the result and make the relation to n-dimensional differentiation and the n-cube clearer:

[math]\displaystyle{ \int_0^a x^{n-1}\,dx = \tfrac{1}{n} x^n \qquad n \geq 1. }[/math]
[math]\displaystyle{ \int x^{n-1}\,dx = \tfrac{1}{n} x^n + C \qquad n \neq 0. }[/math]

More generally, these formulae may be given as:

[math]\displaystyle{ \int (ax + b)^n dx= \frac{(ax + b)^{n+1}}{a(n + 1)} + C \qquad\mbox{(for } n\neq -1\mbox{)}\,\! }[/math]
[math]\displaystyle{ \int\frac{1}{ax + b} dx= \frac{1}{a}\ln\left|ax + b\right| + C }[/math]
More generally:
[math]\displaystyle{ \int\frac{1}{ax + b} \, dx= \begin{cases} \frac{1}{a}\ln\left|ax + b\right| + C^- & x \lt -b/a \\ \frac{1}{a}\ln\left|ax + b\right| + C^+ & x \gt -b/a \end{cases} }[/math]

Proof

The modern proof is to use an antiderivative: the derivative of xn is shown to be nxn−1 – for non-negative integers. This is shown from the binomial formula and the definition of the derivative – and thus by the fundamental theorem of calculus the antiderivative is the integral. This method fails for [math]\displaystyle{ \int \frac{1}{x}\,dx, }[/math] as the candidate antiderivative is [math]\displaystyle{ \frac{1}{0} \cdot x^0 }[/math], which is undefined due to division by zero. The logarithm function, which is the actual antiderivative of 1/x, must be introduced and examined separately.

The derivative [math]\displaystyle{ (x^n)'=nx^{n-1} }[/math] can be geometrized as the infinitesimal change in volume of the n-cube, which is the area of n faces, each of dimension n − 1.
Integrating this picture – stacking the faces – geometrizes the fundamental theorem of calculus, yielding a decomposition of the n-cube into n pyramids, which is a geometric proof of Cavalieri's quadrature formula.

For positive integers, this proof can be geometrized:[2] if one considers the quantity xn as the volume of the n-cube (the hypercube in n dimensions), then the derivative is the change in the volume as the side length is changed – this is xn−1, which can be interpreted as the area of n faces, each of dimension n − 1 (fixing one vertex at the origin, these are the n faces not touching the vertex), corresponding to the cube increasing in size by growing in the direction of these faces – in the 3-dimensional case, adding 3 infinitesimally thin squares, one to each of these faces. Conversely, geometrizing the fundamental theorem of calculus, stacking up these infinitesimal (n − 1) cubes yields a (hyper)-pyramid, and n of these pyramids form the n-cube, which yields the formula. Further, there is an n-fold cyclic symmetry of the n-cube around the diagonal cycling these pyramids (for which a pyramid is a fundamental domain). In the case of the cube (3-cube), this is how the volume of a pyramid was originally rigorously established: the cube has 3-fold symmetry, with fundamental domain a pyramids, dividing the cube into 3 pyramids, corresponding to the fact that the volume of a pyramid is one third of the base times the height. This illustrates geometrically the equivalence between the quadrature of the parabola and the volume of a pyramid, which were computed classically by different means.

Alternative proofs exist – for example, Fermat computed the area via an algebraic trick of dividing the domain into certain intervals of unequal length;[3] alternatively, one can prove this by recognizing a symmetry of the graph y = xn under inhomogeneous dilation (by d in the x direction and dn in the y direction, algebraicizing the n dimensions of the y direction),[4] or deriving the formula for all integer values by expanding the result for n = −1 and comparing coefficients.[5]

History

Archimedes computed the area of parabolic segments in his The Quadrature of the Parabola.

A detailed discussion of the history, with original sources, is given in (Laubenbacher Pengelley); see also history of calculus and history of integration.

The case of the parabola was proven in antiquity by the ancient Greek mathematician Archimedes in his The Quadrature of the Parabola (3rd century BC), via the method of exhaustion. Of note is that Archimedes computed the area inside a parabola – a so-called "parabolic segment" – rather than the area under the graph y = x2, which is instead the perspective of Cartesian geometry. These are equivalent computations, but reflect a difference in perspective. The Ancient Greeks, among others, also computed the volume of a pyramid or cone, which is mathematically equivalent.

In the 11th century, the Islamic mathematician Ibn al-Haytham (known as Alhazen in Europe) computed the integrals of cubics and quartics (degree three and four) via mathematical induction, in his Book of Optics.[6]

The case of higher integers was computed by Cavalieri for n up to 9, using his method of indivisibles (Cavalieri's principle).[7] He interpreted these as higher integrals as computing higher-dimensional volumes, though only informally, as higher-dimensional objects were as yet unfamiliar.[8] This method of quadrature was then extended by Italian mathematician Evangelista Torricelli to other curves such as the cycloid, then the formula was generalized to fractional and negative powers by English mathematician John Wallis, in his Arithmetica Infinitorum (1656), which also standardized the notion and notation of rational powers – though Wallis incorrectly interpreted the exceptional case n = −1 (quadrature of the hyperbola) – before finally being put on rigorous ground with the development of integral calculus.

Prior to Wallis's formalization of fractional and negative powers, which allowed explicit functions [math]\displaystyle{ y=x^{p/q}, }[/math] these curves were handled implicitly, via the equations [math]\displaystyle{ x^p=ky^q }[/math] and [math]\displaystyle{ x^py^q=k }[/math] (p and q always positive integers) and referred to respectively as higher parabolae and higher hyperbolae (or "higher parabolas" and "higher hyperbolas"). Pierre de Fermat also computed these areas (except for the exceptional case of −1) by an algebraic trick – he computed the quadrature of the higher hyperbolae via dividing the line into equal intervals, and then computed the quadrature of the higher parabolae by using a division into unequal intervals, presumably by inverting the divisions he used for hyperbolae.[9] However, as in the rest of his work, Fermat's techniques were more ad hoc tricks than systematic treatments, and he is not considered to have played a significant part in the subsequent development of calculus.

Of note is that Cavalieri only compared areas to areas and volumes to volumes – these always having dimensions, while the notion of considering an area as consisting of units of area (relative to a standard unit), hence being unitless, appears to have originated with Wallis;[10][11] Wallis studied fractional and negative powers, and the alternative to treating the computed values as unitless numbers was to interpret fractional and negative dimensions.

The exceptional case of −1 (the standard hyperbola) was first successfully treated by Grégoire de Saint-Vincent in his Opus geometricum quadrature circuli et sectionum coni (1647), though a formal treatment had to wait for the development of the natural logarithm, which was accomplished by Nicholas Mercator in his Logarithmotechnia (1668).

References

  1. "Reader Survey: log|x| + C", Tom Leinster, The n-category Café, March 19, 2012
  2. (Barth 2004), (Carter Champanerkar)
  3. See Rickey.
  4. (Wildberger 2002)
  5. (Bradley 2003)
  6. Victor J. Katz (1995), "Ideas of Calculus in Islam and India", Mathematics Magazine 68 (3): 163–174 [165–9 & 173–4]
  7. (Struik 1986)
  8. (Laubenbacher Pengelley) – see Informal pedagogical synopsis of the Analysis chapter for brief form
  9. See Rickey reference for discussion and further references.
  10. Ball, 281
  11. Britannica, 171

History

Proofs

  • Wildberger, N. J. (2002). "A new proof of Cavalieri's quadrature formula". The American Mathematical Monthly 109 (9): 843–845. doi:10.2307/3072373. 
  • Bradley, David M. (May 2003). "Remark on Cavalieri's quadrature formula". The American Mathematical Monthly 110 (5): 437. Bibcode2005math......5059B.  Appears in print at end of:
    • Sondow, Jonathan (May 2003). "Zeros of the Alternating Zeta Function on the Line R(S) = 1". The American Mathematical Monthly 110 (5): 435-437. 
  • Barth, N. R. (2004). "Computing Cavalieri's quadrature formula by a symmetry of the n-cube". The American Mathematical Monthly 111 (9): 811–813. doi:10.2307/4145193. 
  • Carter, J. Scott; Champanerkar, Abhijit (2006). "A geometric method to compute some elementary integrals". arXiv:math/0608722.
  • Malik, M.A. (1984) "A Note on Cavalieri Integration", Mathematics Magazine 57(3): 154–6 doi:10.2307/2689662
  • V. Frederick Rickey (2011) Fermat's Integration of Powers", in Historical Notes for Calculus Teachers

External links