TRAPEZOIDAL METHOD Let f x) have two continuous derivatives on
TRAPEZOIDAL METHOD ERROR FORMULA
Theorem Let f (x) have two continuous derivatives on
the interval a x b. Then
EnT
(f
)
Zb a
f
(x)
dx
-
Tn(f )
=
h2 -
(b - 12
a) f
00
(cn)
for some cn in the interval [a, b].
Later I will say something about the proof of this result, as it leads to some other useful formulas for the error.
The above formula says that the error decreases in a manner that is roughly proportional to h2. Thus doubling n (and halving h) should cause the error to decrease by a factor of approximately 4. This is what we observed with a past example from the preceding section.
Example. Consider evaluating Z 2 dx
I = 0 1 + x2
using the trapezoidal method Tn(f ). How large should n be chosen in order to ensure that
???EnT (f )??? 5 ? 10-6
We begin by calculating the derivatives:
f 0(x)
=
-2x
?
1
+
x2?2
,
f 00(x)
=
-2 + 6x2
?
1
+
x2?3
From a graph of f 00(x), max ???f 00(x)??? = 2
0x2
Recall that b - a = 2. Therefore,
EnT (f )
=
h2 -
(b - 12
a) f
00
(cn)
???EnT (f )???
h2 (2)
h2
?2=
12
3
EnT (f )
=
h2 -
(b - 12
a) f
00
(cn)
We
bound???E??fnT0(0 f(c)n??? )??sinhc1e222w?e2
= do
h2 3 not
know
cn,
and
therefore we must assume the worst possible case, that
which makes the error formula largest. That is what
has been done above.
When do we have
???EnT (f )??? 5 ? 10-6
(1)
To ensure this, we choose h so small that
h2 5 ? 10-6 3 This is equivalent to choosing h and n to satisfy
h .003873 2 n = 516.4 h
Thus n 517 will imply (1).
DERIVING THE ERROR FORMULA
There are two stages in deriving the error: (1) Obtain the error formula for the case of a single subinterval (n = 1); (2) Use this to obtain the general error formula given earlier.
For the trapezoidal method with only a single subin-
terval, we have
Z +h f (x) dx - h [f () + f ( + h)] = -h3f 00(c)
2
12
for some c in the interval [, + h].
A sketch of the derivation of this error formula is given in the problems.
Recall that the general trapezoidal rule Tn(f ) was obtained by applying the simple trapezoidal rule to a subdivision of the original interval of integration. Recall defining and writing
h
=
b
-
a ,
n
xj = a + j h, j = 0, 1, ..., n
xZn
I = f (x) dx
x0
xZ1
xZ2
= f (x) dx + f (x) dx + ? ? ?
x0
x1
xZn
+ f (x) dx
xn-1
I
h 2
[f (x0)
+
f (x1)]
+
h 2
[f (x1)
+
f (x2)]
+???
+
h 2
[f
(xn-2)
+
f (xn-1)]
+
h 2
[f (xn-1)
+
f
(xn)]
................
................
In order to avoid copyright disputes, this page is only a partial summary.
To fulfill the demand for quickly locating and searching documents.
It is intelligent file search solution for home and business.
Related download
- techniques of integration whitman college
- solution to math 2433 calculus iii term exam 3
- integration by substitution
- expected value the expected value of a random variable
- table of integrals
- the squeeze theorem
- homework assignment 4 elizabethtown college
- chapter 2 simple linear regression
- square roots via newton s method mit mathematics
- trapezoidal method let f x have two continuous derivatives on