Transcript Chapter 8

Chapter 8
Numerical Integration
Lecture (I)1
Ref.: “Applied Numerical Methods with MATLAB for Engineers and
Scientists”, Steven Chapra, 2nd ed., Ch. 17, McGraw Hill, 2008.
1
Dr. Jie Zou PHY3320
1
Outline

Introduction




What is integration?
When do we need numerical integration?
Applications of integration in engineering and
science
Newton-cotes formulas

(1) The trapezoidal rule



Error of the Trapezoidal rule
The composite trapezoidal rule
Implementation in MATLAB
Dr. Jie Zou PHY3320
2
Introduction

What is integration?




Ref. Fig. 17.1 Graphical
representation of the
integral
Mathematically: A definite
integration is
b
represented by I   f x dx .
a
It means: The total value, or summation, of
f(x)dx over the range x = a to b.
Graphical representation: For functions lying
above the x axis, the integral corresponds to
the area under the curve of f(x) between x =
a and b.
When do we need numerical integration
(also referred to as quadrature)?


Functions that are difficult to or cannot be
integrated analytically.
Only a table of discrete data are available.
Dr. Jie Zou PHY3320
3
Applications of integration in
engineering and science


Examples related to “the integral as
the area under a curve”:
Examples related to the analogy
between integration and summation:

Ref. Fig. 17.3
An example: To determine the mean of a
continuous function
4
Newton-cotes formulas

Ref. Fig. 17.4 The approx. of an integral by the
area under (a) a straight line and (b) a parabola
Ref. Fig. 17.5 The
approx. of an integral
by the area under
three straight-line
segments
Basic strategy:
Replacing a
complicated function or
tabulated data with a
polynomial that is
easy to integrate.
b
b
I   f x dx   f n x dx
a
a


fn(x) = a0 + a1x + … +
an-1xn-1+anxn
n: The order of the
polynomial.
Dr. Jie Zou PHY3320
5
Newton-cotes formulas: (1)
The trapezoidal rule

Ref. Fig. 17.7 Single application

Basic idea: Replacing the complicated
function or tabulated data with a
polynomial or a series of polynomials of
the first order (linear).
Single and Composite applications

Single application formula:
f a   f b
I  b  a 

2
Width 
Average Height

Composite application formula:
n 1
f x0   2 f xi   f xn 
Ref. Fig. 17.9 Composite application
i 1
I  b  a 

2n

Width 
Average Height
6
Error of the trapezoidal rule

For single applications, an estimate
for the error:
1
1
3
3





Et  
f  b  a ; Ea  
f b  a 
12
12
If the function being integrated is linear,
Et = 0; otherwise, Et  0.


Ref. Fig. 17.8 Truncation
error for a single application
of the trapezoidal rule
For composite applications, an
estimate for the error:
n
Et   Et ,i
i 1


3

b  a

12n
3
n
 f   ; E
i 1
i
a
3

b  a

12n
2
If the number of segments is doubled, Et
is approximately quartered.
b
Here, f   1
f  x dx
b  a a
f 
7
Example: Composite application
of the trapezoidal rule

Example 17.2 (Ref.): Use the
two-segment trapezoidal rule
to estimate the integral of f(x)
= 0.2 + 25x – 200x2 + 675x3
– 900x4 + 400x5 from a = 0 to
b = 0.8. Also, find the true
error Et and the approximate
error, Ea.

x0 = a
x1
Two segments n = 2;
x2 = b

(1) By hand.
(2) Implement on a computerwrite an M-file.
Dr. Jie Zou PHY3320
8
Results
n
2
3
4
5
6
7
8
9
10
h
0.4
0.2667
0.2
0.16
0.1333
0.1143
0.1
0.0889
0.08
Dr. Jie Zou PHY3320
I
1.0688
1.3695
1.4848
1.5399
1.5703
1.5887
1.6008
1.6091
1.6150
t (%)
34.9
16.5
9.5
6.1
4.3
3.2
2.4
1.9
1.6
9
Implementation of composite
trapezoidal rule on a computer


Write an M-file called
My_Trapezoidal_Rule.m to do
Example 17.2.
A copy of the code will be handed out
later.
Dr. Jie Zou PHY3320
10