



Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
The limitations of using Newton-Cotes formulas and high degree polynomials for numerical integration on large intervals. The main idea is to divide the integration interval into subintervals and apply simple integration rules, such as Simpson's rule or the Trapezoidal rule, to each subinterval. examples and comparisons of the accuracy of these methods, as well as error terms for the Composite Trapezoidal rule, Composite Simpson's rule, and Composite Midpoint rule.
Typology: Summaries
1 / 5
This page cannot be seen from the preview
Don't miss anything!




4.4 Composite Numerical Integration
Motivation: 1) on large interval, use Newton-Cotes formulas are not accurate.
2) on large interval, interpolation using high degree polynomial is unsuitable because of oscillatory nature of high degree
polynomials.
Main idea: divide integration interval into subintervals and use simple integration rule for each subinterval.
Example a) Use Simpson’s rule to approximate ∫
. b) Divide into. Use
Simpson’s rule to approximate ∫
and ∫
. Then approximate ∫
by adding approximations
for ∫
and ∫
. Compare with accurate value.
Solution:
a) ∫
Error= | - | = 3.
b) ∫ ∫ ∫ ∫ ∫
Error=| - | = 0.
b) is much more accurate than a).
Composite Trapezoidal rule
Let , and for.
On each subinterval [ ], for for , apply Trapezoidal rule:
Figure 1 Composite Trapezoidal Rule
Error, which can be simplified
( )
( )
( )
( )
( )
Theorem 4.4 Let , and for each. There exists a ( )
for which Composite Simpson’s rule with its error term is
Error Term
Error, which can be simplified
𝑏
𝑎
𝑗
(
𝑛
)
𝑗
𝑗
(
𝑛
)
𝑗
( )
Composite Midpoint rule
Theorem 4.6 Let , and ( ) for each. There exists a
( ) for which Composite Midpoint rule with its error term is
Error Term
Figure 3 Composite Midpoint rule
𝑏
𝑎
𝑗
(
𝑛
)
𝑗