The trapezoidal rule works by approximating the region under the graph of the function as a trapezoid and calculating its area. It follows that
The integral can be even better approximated by partitioning the integration interval, applying the trapezoidal rule to each subinterval, and summing the results. In practice, this "chained" (or "composite") trapezoidal rule is usually what is meant by "integrating with the trapezoidal rule". Let be a partition of such that and be the length of the -th subinterval (that is, ), then
The approximation becomes more accurate as the resolution of the partition increases (that is, for larger ). When the partition has a regular spacing, as is often the case, the formula can be simplified for calculation efficiency.
As discussed below, it is also possible to place error bounds on the accuracy of the value of a definite interval estimated using a trapezoidal rule.
When the grid spacing is non-uniform, one can use the formula
For a domain discretized into equally spaced panels, considerable simplification may occur. Let
the approximation to the integral becomes
which requires fewer evaluations of the function to calculate.
The error of the composite trapezoidal rule is the difference between the value of the integral and the numerical result:
There exists a number ξ between a and b, such that
It follows that if the integrand is concave up (and thus has a positive second derivative), then the error is negative and the trapezoidal rule overestimates the true value. This can also be seen from the geometric picture: the trapezoids include all of the area under the curve and extend over it. Similarly, a concave-down function yields an underestimate because area is unaccounted for under the curve, but none is counted above. If the interval of the integral being approximated includes an inflection point, the error is harder to identify.
In general, three techniques are used in the analysis of error:
An asymptotic error estimate for N → ∞ is given by
Further terms in this error estimate are given by the Euler–Maclaurin summation formula.
It is argued that the speed of convergence of the trapezoidal rule reflects and can be used as a definition of classes of smoothness of the functions.
The trapezoidal rule converges rapidly for periodic functions. This is an easy consequence of the Euler-Maclaurin summation formula, which says that if is times continuously differentiable with period
where and is the periodic extension of the th Bernoulli polynomial. Due to the periodicity, the derivatives at the endpoint cancel and we see that the error is .
Although some effort has been made to extending the Euler-Maclaurin summation formula to higher dimensions, the most straightforward proof of the rapid convergence of the trapezoidal rule in higher dimensions is to reduce the problem to that of convergence of Fourier series. This line of reasoning shows that if is periodic on a -dimensional space with continuous derivatives, the speed of convergence is . For very large dimension, the shows that Monte-Carlo integration is most likely a better choice, but for 2 and 3 dimensions, equispaced sampling is efficient. This is exploited in computational solid state physics where equispaced sampling over primitive cells in the reciprocal lattice is known as Monkhorst-Pack integration..
This section needs expansion. You can help by adding to it. (January 2010)
Applicability and alternatives
The trapezoidal rule is one of a family of formulas for numerical integration called Newton–Cotes formulas, of which the midpoint rule is similar to the trapezoid rule. Simpson's rule is another member of the same family, and in general has faster convergence than the trapezoidal rule for functions which are twice continuously differentiable, though not in all specific cases. However, for various classes of rougher functions (ones with weaker smoothness conditions), the trapezoidal rule has faster convergence in general than Simpson's rule.
For non-periodic functions, however, methods with unequally spaced points such as Gaussian quadrature and Clenshaw–Curtis quadrature are generally far more accurate; Clenshaw–Curtis quadrature can be viewed as a change of variables to express arbitrary integrals in terms of periodic integrals, at which point the trapezoidal rule can be applied accurately.
- Ossendrijver, Mathieu (Jan 29, 2016). "Ancient Babylonian astronomers calculated Jupiter's position from the area under a time-velocity graph". Science. 351: 482–484. doi:10.1126/science.aad8085. PMID 26823423.
- Atkinson (1989, equation (5.1.7))
- (Weideman 2002, p. 23, section 2)
- Atkinson (1989, equation (5.1.9))
- Atkinson (1989, p. 285)
- (Rahman & Schmeisser 1990)
- Kress, Rainer (1998). Numerical Analysis, volume 181 of Graduate Texts in Mathematics. Springer-Verlag.
- "Euler-Maclaurin Summation Formula for Multiple Sums". math.stackexchange.com.
- Thompson, Nick. "Numerical Integration over Brillouin Zones". bandgap.io. Retrieved 19 December 2017.
- (Cruz-Uribe & Neugebauer 2002)
- (Weideman 2002)
- Atkinson, Kendall E. (1989), An Introduction to Numerical Analysis (2nd ed.), New York: John Wiley & Sons, ISBN 978-0-471-50023-0
- Rahman, Qazi I.; Schmeisser, Gerhard (December 1990), "Characterization of the speed of convergence of the trapezoidal rule", Numerische Mathematik, 57 (1): 123–138, doi:10.1007/BF01386402, ISSN 0945-3245
- Burden, Richard L.; Faires, J. Douglas (2000), Numerical Analysis (7th ed.), Brooks/Cole, ISBN 0-534-38216-9
- Weideman, J. A. C. (January 2002), "Numerical Integration of Periodic Functions: A Few Examples", The American Mathematical Monthly, 109 (1): 21–36, doi:10.2307/2695765, JSTOR 2695765
- Cruz-Uribe, D.; Neugebauer, C. J. (2002), "Sharp Error Bounds for the Trapezoidal Rule and Simpson's Rule" (PDF), Journal of Inequalities in Pure and Applied Mathematics, 3 (4)
|The Wikibook A-level Mathematics has a page on the topic of: Trapezium Rule|
- Trapezium formula. I.P. Mysovskikh, Encyclopedia of Mathematics, ed. M. Hazewinkel
- Notes on the convergence of trapezoidal-rule quadrature
- An implementation of trapezoidal quadrature provided by Boost.Math