Lecture 26 Numerical Integration

1
Lecture 26 - Numerical Integration

• CVEN 302
• October 26, 2001

2
Lectures Goals

• Trapezoidal Rule
• Simpsons Rule
• 1/3 Rule
• 3/8 Rule
• Midpoint
• Gaussian Quadrature

Basic Numerical Integration
3
Basic Numerical Integration
We want to find integration of functions of
various forms of the equation known as the Newton
Cotes integration formulas.
4
Basic Numerical Integration

• Weighted sum of function values

f(x)
x
x0
x1
xn
xn-1
5
Numerical Integration
Idea is to do integral in small parts, like the
way you first learned integration - a
summation Numerical methods just try to make it
faster and more accurate
6
Numerical Integration

• Newton-Cotes Closed Formulae -- Use both end
  points
• Trapezoidal Rule Linear
• Simpsons 1/3-Rule Quadratic
• Simpsons 3/8-Rule Cubic
• Booles Rule Fourth-order
• Newton-Cotes Open Formulae -- Use only interior
  points
• midpoint rule

7
Trapezoid Rule

• Straight-line approximation

f(x)
L(x)
x
x0
x1
8
ExampleTrapezoid Rule

• Evaluate the integral
• Exact solution
• Trapezoidal Rule

9
Simpsons 1/3-Rule

• Approximate the function by a parabola

L(x)
f(x)
x
x0
x1
x2
h
h
10
Simpsons 3/8-Rule

• Approximate by a cubic polynomial

f(x)
L(x)
x
x0
x1
x2
h
h
x3
h
11
Example Simpsons Rules

• Evaluate the integral
• Simpsons 1/3-Rule
• Simpsons 3/8-Rule

12
Midpoint Rule

• Newton-Cotes Open Formula

f(x)
x
a
b
xm
13
Two-point Newton-Cotes Open Formula

• Approximate by a straight line

f(x)
x
x0
x1
x2
h
h
x3
h
14
Three-point Newton-Cotes Open Formula

• Approximate by a parabola

f(x)
x
x0
x1
x2
h
h
x3
h
h
x4
15
Better Numerical Integration

• Composite integration
• Composite Trapezoidal Rule
• Composite Simpsons Rule
• Richardson Extrapolation
• Romberg integration

16
Apply trapezoid rule to multiple segments over
integration limits
17
Composite Trapezoid Rule
f(x)
x
x0
x1
x2
h
h
x3
h
h
x4
18
Composite Trapezoid Rule

• Evaluate the integral

19
Composite Trapezoid Example
20
Composite Trapezoid Rule with Unequal Segments

• Evaluate the integral
• h1 2, h2 1, h3 0.5, h4 0.5

21
Composite Simpsons Rule
Piecewise Quadratic approximations
f(x)
...
x
x0
x2
x4
h
h
xn-2
h
xn
h
x3
x1
xn-1
22
Composite Simpsons Rule

• Multiple applications of Simpsons rule

23
Composite Simpsons Rule

• Evaluate the integral
• n 2, h 2
• n 4, h 1

24
Composite Simpsons Example
25
Composite Simpsons Rule with Unequal Segments

• Evaluate the integral
• h1 1.5, h2 0.5

26
Summary

• Integration Techniques
• Trapezoidal Rule Linear
• Simpsons 1/3-Rule Quadratic
• Simpsons 3/8-Rule Cubic

27
Homework

• Check the Homework webpage