Trapezium rule
|
|
Trapezoidal_rule_illustration.png
In mathematics, the trapezium rule (the British term) or trapezoidal rule (the American term) is a way to approximately calculate the definite integral
- <math> \int_{a}^{b} f(x)dx.<math>
The trapezium rule works by approximating the region under the graph of the function <math>f(x)<math> by a trapezium and calculating its area. It follows that
- <math> \int_{a}^{b} f(x) dx \approx (b-a)\frac{f(a) + f(b)}{2}.<math>
To calculate this integral more accurately, one first splits the interval of integration <math>[a, b]<math> into n smaller subintervals, and then applies the trapezium rule on each of them. One obtains the composite trapezium rule:
- <math>\int_a^b f(x)\,dx \approx \frac{b-a}{n} \left( {f(a) + f(b) \over 2} + \sum_{k=1}^{n-1} f \left( a+k \frac{b-a}{n} \right) \right).<math>
This can alternatively be written as:
- <math>\int_a^b f(x)\,dx \approx \frac{b-a}{2n} \left(f(x_0) + f(x_n) + 2f(x_1) + 2f(x_2)+\dots+2f(x_{n-1}) \right).<math>
The trapezium rule is one of a family of formulas for numerical integration called Newton-Cotes formulas. Simpson's rule is another, often more accurate, member of the same family. Simpson's rule and other like methods can be expected to improve on the trapezium rule for functions which are twice continuously differentiable; however for rougher functions the trapezium rule is likely to prove preferable. Moreover, the trapezium rule tends to become extremely accurate when periodic functions are integrated over their periods, a fact best understood in connection with the Euler-Maclaurin summation formula.
References
- Richard L. Burden, J. Douglas Faires (2000), "Numerical Analysis, (7th Ed)", Brooks/Cole. ISBN 0534382169fr:Méthode des trapčzes