Chapter 6: Techniques of Integration
Section 6.7: Numeric Methods
Essentials
Table 6.7.1 lists five common fixed-step methods for the numeric evaluation of the definite integral ∫abfx dx. If h=b−a/n, the values of the integrand fx at the evenly-spaced nodes xk=a+k h,k=0,1,…,n, are denoted by fk rather than by fxk. Of course, x0=a and xn=b. The third column provides an upper bound on the absolute value of the error made when evaluating an integral by one of the listed methods. The maximum of the absolute value of the appropriate derivative is taken over x∈a,b.
Method
Rule
Error Bound
Left Endpoint
h∑k=0n−1fk
h2maxx(|f′|) b−a
Right Endpoint
h∑k=1nfk
Midpoint
h∑k=1nfxk−1+xk2
h224 maxx(|f″|) b−a
Trapezoid
h2f0+fn+2∑k=1n−1fk
h212 maxx(|f″|) b−a
Simpson-13
h3f0+fn+4∑k=1n/2f2 k−1+2∑k=1n/2−1f2 k
h4180 maxxf4b−a
Table 6.7.1 Common fixed-step methods for numeric evaluation of definite integrals
The Left Endpoint, Right Endpoint, and Midpoint rules are examples of "Rectangular" rules because rectangles are used to approximate the area under the graph of the function being integrated. The Trapezoid rule uses trapezoids, and Simpson's rule approximates the function with a sequence of quadratics and finds the areas under these arcs.
The Rectangular rules evaluate the function at n points, but Simpson's rule (as stated in Table 6.7.1) and the Trapezoid rule evaluate the function at n+1 points.
The traditional statement of Simpson's rule given in Table 6.7.1 requires n to be an even integer. This form of the rule is more easily remembered if described as "h/3 times the sum of the first and last function values plus four times the sum of the odd-indexed values plus two times the sum of the even-indexed values." Maple has both this form and an alternate form in which the function is evaluated at the midpoints of the n subintervals generated by the n nodes. In this form, the function is evaluated at 2 n+1 points, so when making comparisons of the accuracy of the various methods, it is important to do so with methods that expend comparable computational energy.
Examples
Example 6.7.1
Use the Trapezoid rule with n=10 to approximate the definite integral ∫141+sinxlnx+1 ⅆx.
Example 6.7.2
If λ is the "actual value" of the definite integral in Example 6.7.1, determine empirically (trial-and-error) the smallest partition for which λ and the Trapezoid rule's value agree when rounded to four places.
Example 6.7.3
Use the error bound in Table 6.7.1 to estimate the value of the partition n for which the Trapezoid rule makes an absolute error of no more than 10−4 when estimating λ, the value of the definite integral in Example 6.7.1.
What is the actual value of n for which the Trapezoid rule achieves this accuracy?
Example 6.7.4
Use Simpson's rule with n=10 to approximate the definite integral ∫141+sinxlnx+1 ⅆx.
Example 6.7.5
If λ is the "actual value" of the definite integral in Example 6.7.1, determine empirically (trial-and-error) the smallest partition for which λ and the value given by Simpson's rule agree when rounded to four places.
Example 6.7.6
Use the error bound in Table 6.7.1 to estimate the value of the partition n for which Simpson's rule makes an absolute error of no more than 10−4 when estimating λ, the value of the definite integral in Example 6.7.1.
What is the actual value of n for which Simpson's rule achieves this accuracy?
Example 6.7.7
Use the Trapezoid rule to approximate the area under the curve determined by the following data points.
k
0
1
2
3
4
5
6
7
8
9
10
xk
−1
−0.8
−0.6
−0.4
−0.2
0.2
0.4
0.6
0.8
fk
13
7.523
4.679
3.528
3.436
4.971
6.175
7.440
8.516
Example 6.7.8
Sketch a derivation of the Trapezoid rule.
Example 6.7.9
Sketch a derivation of Simpson's rule.
<< Previous Section Table of Contents Next Chapter >>
© Maplesoft, a division of Waterloo Maple Inc., 2024. All rights reserved. This product is protected by copyright and distributed under licenses restricting its use, copying, distribution, and decompilation.
For more information on Maplesoft products and services, visit www.maplesoft.com
Download Help Document