Banner A. The Calculus Lifesaver: All the Tools You Need to Excel at Calculus
Подождите немного. Документ загружается.
706 • Estimating Integrals
You still need to choose the numbers c
j
, but things are a lot simpler. For
example, let’s estimate our integral
PSfrag
replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shado
w
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror
(y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y =
(x − 1)
2
−
1
x
Same
height
−
x
Same
length,
opp
osite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y =
2
x
y =
10
x
y =
2
−x
y =
log
2
(x)
4
3
units
mirror
(x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
=
0 radians
θ
hyp
otenuse
opp
osite
adjacen
t
0
(≡ 2π)
π
2
π
3
π
2
I
I
I
I
II
IV
θ
(
x, y)
x
y
r
7
π
6
reference
angle
reference
angle =
π
6
sin
+
sin −
cos
+
cos −
tan
+
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this
angle is
5π
6
clo
ckwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
cos(x)
−
π
2
π
2
y =
tan(x), −
π
2
<
x <
π
2
0
−
π
2
π
2
y =
tan(x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y =
sec(x)
y =
csc(x)
y =
cot(x)
y = f(
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y =
tan
−1
(x)
π
2
π
y =
sin(
x)
x
,
x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 <
x < 0.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x +
2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f(
x)
a
tangen
t at x = a
b
tangen
t at x = b
c
tangen
t at x = c
y = x
2
tangen
t
at x = −
1
u
v
uv
u +
∆u
v +
∆v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f(
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f (
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y =
sin(x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin(
x)
x
π
2
π
1
−
1
x =
0
a =
0
x
> 0
a
> 0
x
< 0
a
< 0
rest
position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f(
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y =
ln(x)
y =
cosh(x)
y =
sinh(x)
y =
tanh(x)
y =
sech(x)
y =
csch(x)
y =
coth(x)
1
−
1
y = f(
x)
original
function
in
verse function
slop
e = 0 at (x, y)
slop
e is infinite at (y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
sin(x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y =
sin
−1
(x)
y =
cos(x)
π
π
2
y =
cos
−1
(x)
−
π
2
1
x
α
β
y =
tan(x)
y =
tan(x)
1
y =
tan
−1
(x)
y =
sec(x)
y =
sec
−1
(x)
y =
csc
−1
(x)
y =
cot
−1
(x)
1
y =
cosh
−1
(x)
y =
sinh
−1
(x)
y =
tanh
−1
(x)
y =
sech
−1
(x)
y =
csch
−1
(x)
y =
coth
−1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y =
1
2
(2
, 3)
y = f(
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Lo
cal maximum
Lo
cal minimum
Horizon
tal point of inflection
1
e
y = f
0
(
x)
y = f (
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2
+
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln
(x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2
+
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing
fence
new
fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallo
w
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f (
x)
11
y = L
(x)
y = f(
x)
F
P
a
a +
∆x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
d
f
∆
x
a
b
y = f(
x)
true
zero
starting
approximation
b
etter approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f(
x)
1
2
y = x
5
0
−
2
y =
1
a
b
y =
sin(x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
=
2
width
of each interval =
2
n
−
2
1
3
0
I
I
I
I
II
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
I
I
I
Ia
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f(
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
a
v
y = f
a
v
c
A
M
0
1
2
a
b
x
t
y = f (
t)
F (
x )
y = f (
t)
F (
x + h)
x + h
F (
x + h) − F (x)
f(
x)
1
2
y =
sin(x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y =
ln|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f(
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
ev
en smaller
y = g(
x)
infinite
area
finite
area
1
y =
1
x
y =
1
x
p
, p
< 1 (typical)
y =
1
x
p
, p
> 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f(
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t =
0
t = π
/6
t = π
/4
t = π
/3
t = π
/2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t =
0
t =
1/2
t =
3/2
t =
2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong
point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r =
3 sin(θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r =
1 + 2 cos(θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r =
1 + cos(θ)
r =
1 +
3
4
cos(
θ)
−
1
4
r =
sin(2θ)
r =
sin(3θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1
+ sin(θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−
4
−
5
4
5
f(
θ)
f(
θ + dθ)
θ
dθ
θ + dθ
appro
ximating region
exact
region
0 ≤ θ ≤ 2
π
r = |
1 + 2 cos(θ)|
2
i
2 − 3
i
−
1
θ =
0
θ =
π
4
θ =
π
2
θ =
2
π
3
θ = π
θ =
13
π
12
θ =
3
π
2
θ =
7
π
4
1
= e
0
e
i
π
4
i = e
i
π
2
e
i
2
π
3
−
1 = e
iπ
e
i
13
π
12
−
i = e
i
3π
2
e
i
7
π
4
i
−
i
1
θ
1 − i
2
i
−
2i
2
−
2
6
i
−
6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f(x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
Z
2
0
e
−x
2
dx
using 10 strips of equal width. The width of each strip is h = (2 − 0)/10, or
1/5, and n = 10, so we have
Z
2
0
e
−x
2
dx
∼
=
h ×
n
X
j=1
f(c
j
) =
1
5
10
X
j=1
e
−c
2
j
.
The intervals all have width
1
5
, so starting from 0, we see that they are par-
titioned as follows:
0 <
1
5
<
2
5
<
3
5
<
4
5
< 1 <
6
5
<
7
5
<
8
5
<
9
5
< 2.
If we let c
j
be at the right-hand endpoint in each case, then we’ll have c
1
=
1
5
,
c
2
=
2
5
, and so on up to c
10
= 2. Plugging these numbers into the above
formula, we have
Z
2
0
e
−x
2
dx
∼
=
1
5
e
−(1/5)
2
+ e
−(2/5)
2
+ ··· + e
−(9/5)
2
+ e
−2
2
.
There are ten terms in the sum. Since our function f is decreasing between
0 and 2, and we’ve used the right-hand endpoint for each strip, the above
estimate is an underestimate. (Can you see why?) In any case, you can use a
calculator or computer to find that the above sum is approximately 0.783670
(to six decimal places).
Now, what if you wanted to use the midpoint of each interval, rather than
the left-hand or right-hand boundary? Well, the midpoint of [0,
1
5
] is
1
10
, the
midpoint of [
1
5
,
2
5
] is
3
10
, and so on. So another possible approximation is given
by
Z
2
0
e
−x
2
dx
∼
=
1
5
e
−(1/10)
2
+ e
−(3/10)
2
+ ··· + e
−(17/10)
2
+ e
−(19/10)
2
.
This is approximately 0.882202.
B.2 The Trapezoidal Rule
There’s quite a bit of a burden involved in picking the numbers c
j
. Most of the
time, people choose either the left-hand endpoint or the right-hand endpoint,
but the midpoint is also a common (and reasonable) choice. Here’s another
method for estimating integrals that removes the element of choice (once you
decide to use the method, of course!) while giving even better estimates. It’s
called the trapezoidal rule.
The idea is very simple: allow the tops of the strips to be nonparallel
to the base. The top of each strip will be the line segment joining the two
corresponding points on the curve y = f(x). Here’s a picture illustrating the
difference in the two approaches:
Section B.2: The Trapezoidal Rule • 707
PSfrag replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shadow
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror (
y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y = (
x − 1)
2
−
1
x
Same height
−
x
Same length,
opposite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y = 2
x
y = 10
x
y = 2
−
x
y = log
2
(
x)
4
3 units
mirror (
x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
= 0 radians
θ
h
yp
otenuse
opposite
adjacent
0 (
≡ 2π)
π
2
π
3
π
2
I
II
III
IV
θ
(
x, y)
x
y
r
7
π
6
reference angle
reference angle =
π
6
sin +
sin −
cos +
cos −
tan +
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this angle is
5
π
6
clockwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y = sin(
x)
y = cos(
x)
−
π
2
π
2
y = tan(
x), −
π
2
< x <
π
2
0
−
π
2
π
2
y = tan(
x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y = sec(
x)
y = csc(
x)
y = cot(
x)
y = f(
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y = tan
−
1
(x)
π
2
π
y =
sin(
x)
x
, x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 < x < 0
.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f(
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x + 2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f(
x)
a
tangent at x = a
b
tangent at x = b
c
tangent at x = c
y = x
2
tangent
at x = −
1
u
v
uv
u + ∆
u
v + ∆
v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f(
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f(
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y = sin(
x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin(
x)
x
π
2
π
1
−
1
x = 0
a = 0
x > 0
a > 0
x < 0
a < 0
rest position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f(
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y = ln(
x)
y = cosh(
x)
y = sinh(
x)
y = tanh(
x)
y = sech(
x)
y = csch(
x)
y = coth(
x)
1
−
1
y = f(
x)
original function
inverse function
slope = 0 at (
x, y)
slope is infinite at (
y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y = sin(
x)
y = sin(
x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y = sin
−
1
(x)
y = cos(
x)
π
π
2
y = cos
−
1
(x)
−
π
2
1
x
α
β
y = tan(
x)
y = tan(
x)
1
y = tan
−
1
(x)
y = sec(
x)
y = sec
−
1
(x)
y = csc
−
1
(x)
y = cot
−
1
(x)
1
y = cosh
−
1
(x)
y = sinh
−
1
(x)
y = tanh
−
1
(x)
y = sech
−
1
(x)
y = csch
−
1
(x)
y = coth
−
1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y = 1
2
(2
, 3)
y = f(
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Local maximum
Local minimum
Horizontal point of inflection
1
e
y = f
0
(
x)
y = f(
x) = x ln(x)
−
1
e
?
y = f (
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2 +
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln(
x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2 +
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing fence
new fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallow
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f(
x)
11
y = L
(x)
y = f(
x)
F
P
a
a + ∆
x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
df
∆
x
a
b
y = f(
x)
true zero
starting approximation
better approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f(
x)
1
2
y = x
5
0
−
2
y = 1
a
b
y = sin(
x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
= 2
width of each interval =
2
n
−
2
1
3
0
I
II
III
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
II
IIa
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f(
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
av
y = f
av
c
A
M
0
1
2
a
b
x
t
y = f (
t)
F (
x )
y = f (
t)
F (
x + h)
x + h
F (
x + h) − F (x)
f(
x)
1
2
y = sin(
x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y = ln
|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f(
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
even smaller
y = g(
x)
infinite area
finite area
1
y =
1
x
y =
1
x
p
, p < 1 (typical)
y =
1
x
p
, p > 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f(
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t = 0
t = π/
6
t = π/
4
t = π/
3
t = π/
2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t = 0
t = 1
/2
t = 3
/2
t = 2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r = 3 sin(
θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r = 1 + 2 cos(
θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r = 1 + cos(
θ)
r = 1 +
3
4
cos(
θ)
−
1
4
r = sin(2
θ)
r = sin(3
θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1 + sin(
θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−
4
−
5
4
5
f(
θ)
f(
θ + dθ)
θ
dθ
θ + dθ
approximating region
exact region
0 ≤ θ ≤ 2
π
r = |
1 + 2 cos(θ)|
2
i
2 − 3
i
−
1
θ = 0
θ =
π
4
θ =
π
2
θ =
2
π
3
θ = π
θ =
13
π
12
θ =
3
π
2
θ =
7
π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2
π
3
−
1 = e
iπ
e
i
13
π
12
−
i = e
i
3π
2
e
i
7
π
4
i
−
i
1
θ
1 − i
2
i
−
2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f(x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0 0
1 1
11
2
2
−1
−1
−2
−2
Let’s take a closer look at one of the new strips:
PSfrag replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shadow
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror (
y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y = (
x − 1)
2
−
1
x
Same height
−
x
Same length,
opposite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y = 2
x
y = 10
x
y = 2
−
x
y = log
2
(
x)
4
3 units
mirror (
x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
= 0 radians
θ
hyp
otenuse
opposite
adjacent
0 (
≡ 2π)
π
2
π
3
π
2
I
II
III
IV
θ
(
x, y)
x
y
r
7
π
6
reference angle
reference angle =
π
6
sin +
sin −
cos +
cos −
tan +
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this angle is
5
π
6
clockwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y = sin(
x)
y = cos(
x)
−
π
2
π
2
y = tan(
x), −
π
2
< x <
π
2
0
−
π
2
π
2
y = tan(
x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y = sec(
x)
y = csc(
x)
y = cot(
x)
y = f (
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y = tan
−
1
(x)
π
2
π
y =
sin(
x)
x
, x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 < x < 0
.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x + 2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f (
x)
a
tangent at x = a
b
tangent at x = b
c
tangent at x = c
y = x
2
tangent
at x = −
1
u
v
uv
u + ∆
u
v + ∆
v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f (
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f(
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y = sin(
x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin(
x)
x
π
2
π
1
−
1
x = 0
a = 0
x > 0
a > 0
x < 0
a < 0
rest position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f (
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y = ln(
x)
y = cosh(
x)
y = sinh(
x)
y = tanh(
x)
y = sech(
x)
y = csch(
x)
y = coth(
x)
1
−
1
y = f (
x)
original function
inverse function
slope = 0 at (
x, y)
slope is infinite at (
y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y = sin(
x)
y = sin(
x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y = sin
−
1
(x)
y = cos(
x)
π
π
2
y = cos
−
1
(x)
−
π
2
1
x
α
β
y = tan(
x)
y = tan(
x)
1
y = tan
−
1
(x)
y = sec(
x)
y = sec
−
1
(x)
y = csc
−
1
(x)
y = cot
−
1
(x)
1
y = cosh
−
1
(x)
y = sinh
−
1
(x)
y = tanh
−
1
(x)
y = sech
−
1
(x)
y = csch
−
1
(x)
y = coth
−
1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y = 1
2
(2
, 3)
y = f (
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Local maximum
Local minimum
Horizontal point of inflection
1
e
y = f
0
(
x)
y = f(
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2 +
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln(
x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2 +
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing fence
new fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallow
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f(
x)
11
y = L
(x)
y = f (
x)
F
P
a
a + ∆
x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
df
∆
x
a
b
y = f (
x)
true zero
starting approximation
better approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f (
x)
1
2
y = x
5
0
−
2
y = 1
a
b
y = sin(
x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
= 2
width of each interval =
2
n
−
2
1
3
0
I
II
III
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
II
IIa
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f (
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
av
y = f
av
c
A
M
0
1
2
a
b
x
t
y = f(
t)
F (
x )
y = f(
t)
F (
x + h )
x + h
F (
x + h) − F (x)
f(
x)
1
2
y = sin(
x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y = ln
|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f (
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
even smaller
y = g(
x)
infinite area
finite area
1
y =
1
x
y =
1
x
p
, p < 1 (typical)
y =
1
x
p
, p > 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f (
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t = 0
t = π/
6
t = π/
4
t = π/
3
t = π/
2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t = 0
t = 1
/2
t = 3
/2
t = 2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r = 3 sin(
θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r = 1 + 2 cos(
θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r = 1 + cos(
θ)
r = 1 +
3
4
cos(
θ)
−
1
4
r = sin(2
θ)
r = sin(3
θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1 + sin(
θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−
4
−
5
4
5
f(
θ)
f(
θ + dθ)
θ
dθ
θ + dθ
approximating region
exact region
0 ≤ θ ≤ 2
π
r = |
1 + 2 cos(θ)|
2
i
2 − 3
i
−
1
θ = 0
θ =
π
4
θ =
π
2
θ =
2
π
3
θ = π
θ =
13
π
12
θ =
3
π
2
θ =
7
π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2π
3
−1 = e
iπ
e
i
13π
12
−i = e
i
3π
2
e
i
7π
4
i
−i
1
θ
1 − i
2i
−2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f (x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
Since there are two parallel sides, the strip is a trapezoid. The base length is
(x
j
− x
j−1
) units, whereas the heights of the parallel sides are f (x
j−1
) and
f(x
j
) units. By the formula for the area of a trapezoid, we see that the area
of this trapezoidal strip is
1
2
(f(x
j−1
) + f(x
j
))(x
j
− x
j−1
) square units. If we
also make sure that we always take our partition to be evenly spaced, then as
in the previous section, we see that x
j
−x
j−1
is just (b−a)/n. This is exactly
the strip width (in units), which we called h, so the area of one strip becomes
h
2
(f(x
j−1
) + f(x
j
))
square units. All that’s left is to add up the areas of all the trapezoidal strips.
We could just whack a sigma sign outside the above quantity, pulling out the
constant factor h/2, like this:
Z
b
a
f(x) dx
∼
=
h
2
n
X
j=1
(f(x
j−1
) + f(x
j
)).
Actually, we can simplify this expression quite a bit. You see, except for the
leftmost and rightmost strips, every other pair of adjacent strips shares an
edge, like this:
PSfrag replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shadow
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror (
y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y = (
x − 1)
2
−
1
x
Same height
−
x
Same length,
opposite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y = 2
x
y = 10
x
y = 2
−
x
y = log
2
(
x)
4
3 units
mirror (
x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
= 0 radians
θ
hyp
otenuse
opposite
adjacent
0 (
≡ 2π)
π
2
π
3
π
2
I
II
III
IV
θ
(
x, y)
x
y
r
7
π
6
reference angle
reference angle =
π
6
sin +
sin −
cos +
cos −
tan +
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this angle is
5
π
6
clockwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y = sin(
x)
y = cos(
x)
−
π
2
π
2
y = tan(
x), −
π
2
< x <
π
2
0
−
π
2
π
2
y = tan(
x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y = sec(
x)
y = csc(
x)
y = cot(
x)
y = f(
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y = tan
−
1
(x)
π
2
π
y =
sin(
x)
x
, x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 < x < 0
.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x + 2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f(
x)
a
tangent at x = a
b
tangent at x = b
c
tangent at x = c
y = x
2
tangent
at x = −
1
u
v
uv
u + ∆
u
v + ∆
v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f(
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f (
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y = sin(
x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin(
x)
x
π
2
π
1
−
1
x = 0
a = 0
x > 0
a > 0
x < 0
a < 0
rest position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f(
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y = ln(
x)
y = cosh(
x)
y = sinh(
x)
y = tanh(
x)
y = sech(
x)
y = csch(
x)
y = coth(
x)
1
−
1
y = f(
x)
original function
inverse function
slope = 0 at (
x, y)
slope is infinite at (
y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y = sin(
x)
y = sin(
x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y = sin
−
1
(x)
y = cos(
x)
π
π
2
y = cos
−
1
(x)
−
π
2
1
x
α
β
y = tan(
x)
y = tan(
x)
1
y = tan
−
1
(x)
y = sec(
x)
y = sec
−
1
(x)
y = csc
−
1
(x)
y = cot
−
1
(x)
1
y = cosh
−
1
(x)
y = sinh
−
1
(x)
y = tanh
−
1
(x)
y = sech
−
1
(x)
y = csch
−
1
(x)
y = coth
−
1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y = 1
2
(2
, 3)
y = f(
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Local maximum
Local minimum
Horizontal point of inflection
1
e
y = f
0
(
x)
y = f (
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2 +
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln(
x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2 +
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing fence
new fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallow
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f (
x)
11
y = L
(x)
y = f(
x)
F
P
a
a + ∆
x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
df
∆
x
a
b
y = f(
x)
true zero
starting approximation
better approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f(
x)
1
2
y = x
5
0
−
2
y = 1
a
b
y = sin(
x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
= 2
width of each interval =
2
n
−
2
1
3
0
I
II
III
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
II
IIa
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f(
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
av
y = f
av
c
A
M
0
1
2
a
b
x
t
y = f(
t)
F (
x )
y = f(
t)
F (
x + h)
x + h
F (
x + h) − F (x)
f(
x)
1
2
y = sin(
x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y = ln
|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f(
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
even smaller
y = g(
x)
infinite area
finite area
1
y =
1
x
y =
1
x
p
, p < 1 (typical)
y =
1
x
p
, p > 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f(
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t = 0
t = π/
6
t = π/
4
t = π/
3
t = π/
2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t = 0
t = 1
/2
t = 3
/2
t = 2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r = 3 sin(
θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r = 1 + 2 cos(
θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r = 1 + cos(
θ)
r = 1 +
3
4
cos(
θ)
−
1
4
r = sin(2
θ)
r = sin(3
θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1 + sin(
θ)
−
π
4
≤ θ ≤
5π
4
0 ≤ θ ≤ 2π
0 ≤ θ ≤ π
−4
−5
4
5
f(θ)
f(θ + dθ)
θ
dθ
θ + dθ
approximating region
exact region
0 ≤ θ ≤ 2π
r = |1 + 2 cos(θ)|
2i
2 − 3i
−1
θ = 0
θ =
π
4
θ =
π
2
θ =
2π
3
θ = π
θ =
13π
12
θ =
3π
2
θ =
7π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2π
3
−1 = e
iπ
e
i
13π
12
−i = e
i
3π
2
e
i
7π
4
i
−i
1
θ
1 − i
2i
−2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f(x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
f(x
j+1
)
x
j+1
708 • Estimating Integrals
This means that we can collect a lot of the terms together. In particular,
except for x
0
and x
n
, every term of the form f (x
j
) needs to be counted twice.
For example, if n = 4, we have
Z
b
a
f(x) dx
∼
=
h
2
((f(x
0
) + f(x
1
)) + (f(x
1
) + f (x
2
))
+ f (x
2
) + f(x
3
)) + (f(x
3
) + f(x
4
))) .
So we can group all but the first and last terms in the sum into pairs to get
Z
b
a
f(x) dx
∼
=
h
2
(f(x
0
) + 2f (x
1
) + 2f (x
2
) + 2f(x
3
) + f(x
4
)) .
The same trick works in general; we end up with:
Trapezoidal rule: if x
0
< x
1
< ··· < x
n
is an evenly spaced
partition of [a, b], and h = (b − a)/n is the strip width, then
Z
b
a
f(x) dx
∼
=
h
2
f(x
0
) + 2f(x
1
) + 2f(x
2
)
+ ··· + 2f(x
n−2
) + 2f(x
n−1
) + f(x
n
)
.
Let’s use this to find an approximate value of our old integral,
PSfrag
replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shado
w
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror
(y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y =
(x − 1)
2
−
1
x
Same
height
−
x
Same
length,
opp
osite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y =
2
x
y =
10
x
y =
2
−x
y =
log
2
(x)
4
3
units
mirror
(x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
=
0 radians
θ
hyp
otenuse
opp
osite
adjacen
t
0
(≡ 2π)
π
2
π
3
π
2
I
I
I
I
II
IV
θ
(
x, y)
x
y
r
7
π
6
reference
angle
reference
angle =
π
6
sin
+
sin −
cos
+
cos −
tan
+
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this
angle is
5π
6
clo
ckwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
cos(x)
−
π
2
π
2
y =
tan(x), −
π
2
<
x <
π
2
0
−
π
2
π
2
y =
tan(x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y =
sec(x)
y =
csc(x)
y =
cot(x)
y = f(
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y =
tan
−1
(x)
π
2
π
y =
sin(
x)
x
,
x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 <
x < 0.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x +
2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f(
x)
a
tangen
t at x = a
b
tangen
t at x = b
c
tangen
t at x = c
y = x
2
tangen
t
at x = −
1
u
v
uv
u +
∆u
v +
∆v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f(
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f (
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y =
sin(x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin(
x)
x
π
2
π
1
−
1
x =
0
a =
0
x
> 0
a
> 0
x
< 0
a
< 0
rest
position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f(
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y =
ln(x)
y =
cosh(x)
y =
sinh(x)
y =
tanh(x)
y =
sech(x)
y =
csch(x)
y =
coth(x)
1
−
1
y = f(
x)
original
function
in
verse function
slop
e = 0 at (x, y)
slop
e is infinite at (y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
sin(x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y =
sin
−1
(x)
y =
cos(x)
π
π
2
y =
cos
−1
(x)
−
π
2
1
x
α
β
y =
tan(x)
y =
tan(x)
1
y =
tan
−1
(x)
y =
sec(x)
y =
sec
−1
(x)
y =
csc
−1
(x)
y =
cot
−1
(x)
1
y =
cosh
−1
(x)
y =
sinh
−1
(x)
y =
tanh
−1
(x)
y =
sech
−1
(x)
y =
csch
−1
(x)
y =
coth
−1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y =
1
2
(2
, 3)
y = f(
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Lo
cal maximum
Lo
cal minimum
Horizon
tal point of inflection
1
e
y = f
0
(
x)
y = f (
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2
+
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln
(x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2
+
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing
fence
new
fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallo
w
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f (
x)
11
y = L
(x)
y = f(
x)
F
P
a
a +
∆x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
d
f
∆
x
a
b
y = f(
x)
true
zero
starting
approximation
b
etter approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f(
x)
1
2
y = x
5
0
−
2
y =
1
a
b
y =
sin(x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
=
2
width
of each interval =
2
n
−
2
1
3
0
I
I
I
I
II
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
I
I
I
Ia
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f(
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
a
v
y = f
a
v
c
A
M
0
1
2
a
b
x
t
y = f (
t)
F (
x )
y = f (
t)
F (
x + h)
x + h
F (
x + h) − F (x)
f(
x)
1
2
y =
sin(x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y =
ln|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f(
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
ev
en smaller
y = g(
x)
infinite
area
finite
area
1
y =
1
x
y =
1
x
p
, p
< 1 (typical)
y =
1
x
p
, p
> 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f(
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t =
0
t = π
/6
t = π
/4
t = π
/3
t = π
/2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t =
0
t =
1/2
t =
3/2
t =
2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong
point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r =
3 sin(θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r =
1 + 2 cos(θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r =
1 + cos(θ)
r =
1 +
3
4
cos(
θ)
−
1
4
r =
sin(2θ)
r =
sin(3θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1
+ sin(θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−
4
−
5
4
5
f(
θ)
f(
θ + dθ)
θ
dθ
θ + dθ
appro
ximating region
exact
region
0 ≤ θ ≤ 2
π
r = |
1 + 2 cos(θ)|
2
i
2 − 3
i
−
1
θ =
0
θ =
π
4
θ =
π
2
θ =
2
π
3
θ = π
θ =
13π
12
θ =
3π
2
θ =
7π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2π
3
−1 = e
iπ
e
i
13π
12
−i = e
i
3π
2
e
i
7π
4
i
−i
1
θ
1 − i
2i
−2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f(x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
f(x
j+1
)
x
j+1
Z
2
0
e
−x
2
dx.
We’ll take n = 5. Since [0, 2] has length 2 units, the width of each strip is
therefore h =
2
5
units, and the partition is
0 <
2
5
<
4
5
<
6
5
<
8
5
< 2.
According to the trapezoidal rule, we have
Z
2
0
e
−x
2
dx
∼
=
2/5
2
e
−0
2
+ 2e
−(2/5)
2
+ 2e
−(4/5)
2
+ 2e
−(6/5)
2
+ 2e
−(8/5)
2
+ e
−2
2
.
If you want, you can simplify the right-hand side to
1
5
1 + 2e
−4/25
+ 2e
−16/25
+ 2e
−36/25
+ 2e
−64/25
+ e
−4
.
You could also use a calculator or computer to see that the above number
is 0.881131 to six decimal places. This is somewhat less than the estimate
1.0865 which we found at the end of Section B.1 above, but quite close to the
estimate 0.882202 from the end of Section B.1.1.
Section B.3: Simpson’s Rule • 709
B.3 Simpson’s Rule
Why stop at trapezoids? They still have a clunky linear top. We can do
better by using a curve at the top of the strip, rather than a line segment.
Here’s how it’s done. Start by looking at a pair of adjacent strips, but instead
of connecting the tops with two line segments, use a quadratic curve, like this:
PSfrag replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shadow
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror (
y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y = (
x − 1)
2
−
1
x
Same height
−
x
Same length,
opposite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y = 2
x
y = 10
x
y = 2
−
x
y = log
2
(
x)
4
3 units
mirror (
x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
= 0 radians
θ
hyp
otenuse
opposite
adjacent
0 (
≡ 2π)
π
2
π
3
π
2
I
II
III
IV
θ
(
x, y)
x
y
r
7
π
6
reference angle
reference angle =
π
6
sin +
sin −
cos +
cos −
tan +
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this angle is
5
π
6
clockwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y = sin(
x)
y = cos(
x)
−
π
2
π
2
y = tan(
x), −
π
2
< x <
π
2
0
−
π
2
π
2
y = tan(
x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y = sec(
x)
y = csc(
x)
y = cot(
x)
y = f(
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y = tan
−
1
(x)
π
2
π
y =
sin(
x)
x
, x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 < x < 0
.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f(
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x + 2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f(
x)
a
tangent at x = a
b
tangent at x = b
c
tangent at x = c
y = x
2
tangent
at x = −
1
u
v
uv
u + ∆
u
v + ∆
v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f(
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f(
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y = sin(
x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin(
x)
x
π
2
π
1
−
1
x = 0
a = 0
x > 0
a > 0
x < 0
a < 0
rest position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f(
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y = ln(
x)
y = cosh(
x)
y = sinh(
x)
y = tanh(
x)
y = sech(
x)
y = csch(
x)
y = coth(
x)
1
−
1
y = f(
x)
original function
inverse function
slope = 0 at (
x, y)
slope is infinite at (
y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y = sin(
x)
y = sin(
x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y = sin
−
1
(x)
y = cos(
x)
π
π
2
y = cos
−
1
(x)
−
π
2
1
x
α
β
y = tan(
x)
y = tan(
x)
1
y = tan
−
1
(x)
y = sec(
x)
y = sec
−
1
(x)
y = csc
−
1
(x)
y = cot
−
1
(x)
1
y = cosh
−
1
(x)
y = sinh
−
1
(x)
y = tanh
−
1
(x)
y = sech
−
1
(x)
y = csch
−
1
(x)
y = coth
−
1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y = 1
2
(2
, 3)
y = f(
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y = sin(
x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Local maximum
Local minimum
Horizontal point of inflection
1
e
y = f
0
(
x)
y = f(
x) = x ln(x)
−
1
e
?
y = f (
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2 +
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln(
x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2 +
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing fence
new fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallow
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f(
x)
11
y = L
(x)
y = f(
x)
F
P
a
a + ∆
x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
df
∆
x
a
b
y = f(
x)
true zero
starting approximation
better approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f(
x)
1
2
y = x
5
0
−
2
y = 1
a
b
y = sin(
x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
= 2
width of each interval =
2
n
−
2
1
3
0
I
II
III
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
II
IIa
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f(
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
av
y = f
av
c
A
M
0
1
2
a
b
x
t
y = f (
t)
F (
x )
y = f (
t)
F (
x + h )
x + h
F (
x + h) − F (x)
f(
x)
1
2
y = sin(
x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y = ln
|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f(
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
even smaller
y = g(
x)
infinite area
finite area
1
y =
1
x
y =
1
x
p
, p < 1 (typical)
y =
1
x
p
, p > 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f(
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t = 0
t = π/
6
t = π/
4
t = π/
3
t = π/
2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t = 0
t = 1
/2
t = 3
/2
t = 2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r = 3 sin(
θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r = 1 + 2 cos(
θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r = 1 + cos(
θ)
r = 1 +
3
4
cos(
θ)
−
1
4
r = sin(2
θ)
r = sin(3
θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1 + sin(
θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−
4
−
5
4
5
f(
θ)
f(
θ + dθ)
θ
dθ
θ + dθ
approximating region
exact region
0 ≤ θ ≤ 2
π
r = |
1 + 2 cos(θ)|
2
i
2 − 3
i
−
1
θ = 0
θ =
π
4
θ =
π
2
θ =
2
π
3
θ = π
θ =
13
π
12
θ =
3
π
2
θ =
7
π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2π
3
−1 = e
iπ
e
i
13π
12
−i = e
i
3π
2
e
i
7π
4
i
−i
1
θ
1 − i
2i
−2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f(x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
f(x
j+1
)
x
j+1
f(x
0
)
x
0
f(x
1
)
x
1
f(x
2
)
x
2
h
h
As we’ll see in Section B.3.1 below, the shaded area is given by
h
3
(f(x
0
) + 4f (x
1
) + f(x
2
))
square units, where we have once again set h = (b−a)/n. Now if we repeat this
for each pair of strips and add up all the areas, we’ll get our approximation.
As in the case of the trapezoidal rule, adjacent pairs of strips share an edge,
so there’s some doubling up. For example, if there are four strips, then the
total would be
h
3
((f(x
0
) + 4f (x
1
) + f(x
2
)) + (f(x
2
) + 4f (x
3
) + f(x
4
))) ;
the two terms of the form f(x
2
) combine to give 2f(x
2
), so the total is
h
3
(f(x
0
) + 4f(x
1
) + 2f(x
2
) + 4f(x
3
) + f (x
4
)) .
The same pattern persists with more strips, so that the coefficient of f (x
j
) is
equal to 2 if j is even and 4 if j is odd—except for f(x
0
) and f (x
n
), which
both have coefficient 1. All in all, we have:
Simpson’s rule: if n is even, x
0
< x
1
< ··· < x
n
is an evenly
spaced partition of [a, b], and h = (b − a)/n, then
Z
b
a
f(x) dx
∼
=
h
3
f(x
0
) + 4f (x
1
) + 2f (x
2
) + 4f(x
3
)
+ ··· + 2f(x
n−2
) + 4f(x
n−1
) + f(x
n
)
.
Compare this to the trapezoidal rule from the previous section. Instead of
coefficients which look like 1, 2, 2, . . . , 2, 2, 1, this time the coefficients follow
710 • Estimating Integrals
the pattern 1, 4, 2, 4, 2, . . . , 2, 4, 2, 4, 1. Also note that the denominator of the
constant out front is 3 instead of 2.
Simpson’s rule is pretty easy to apply. Let’s go back to our old example
PSfrag
replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shado
w
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror
(y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y =
(x − 1)
2
−
1
x
Same
height
−
x
Same
length,
opp
osite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y =
2
x
y =
10
x
y =
2
−x
y =
log
2
(x)
4
3
units
mirror
(x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
=
0 radians
θ
hyp
otenuse
opp
osite
adjacen
t
0
(≡ 2π)
π
2
π
3
π
2
I
I
I
I
II
IV
θ
(
x, y)
x
y
r
7
π
6
reference
angle
reference
angle =
π
6
sin
+
sin −
cos
+
cos −
tan
+
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this
angle is
5π
6
clo
ckwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
cos(x)
−
π
2
π
2
y =
tan(x), −
π
2
<
x <
π
2
0
−
π
2
π
2
y =
tan(x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y =
sec(x)
y =
csc(x)
y =
cot(x)
y = f(
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y =
tan
−1
(x)
π
2
π
y =
sin(
x)
x
,
x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 <
x < 0.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x +
2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f(
x)
a
tangen
t at x = a
b
tangen
t at x = b
c
tangen
t at x = c
y = x
2
tangen
t
at x = −
1
u
v
uv
u +
∆u
v +
∆v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f(
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f (
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y =
sin(x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin(
x)
x
π
2
π
1
−
1
x =
0
a =
0
x
> 0
a
> 0
x
< 0
a
< 0
rest
position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f(
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y =
ln(x)
y =
cosh(x)
y =
sinh(x)
y =
tanh(x)
y =
sech(x)
y =
csch(x)
y =
coth(x)
1
−
1
y = f(
x)
original
function
in
verse function
slop
e = 0 at (x, y)
slop
e is infinite at (y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
sin(x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y =
sin
−1
(x)
y =
cos(x)
π
π
2
y =
cos
−1
(x)
−
π
2
1
x
α
β
y =
tan(x)
y =
tan(x)
1
y =
tan
−1
(x)
y =
sec(x)
y =
sec
−1
(x)
y =
csc
−1
(x)
y =
cot
−1
(x)
1
y =
cosh
−1
(x)
y =
sinh
−1
(x)
y =
tanh
−1
(x)
y =
sech
−1
(x)
y =
csch
−1
(x)
y =
coth
−1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y =
1
2
(2
, 3)
y = f(
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Lo
cal maximum
Lo
cal minimum
Horizon
tal point of inflection
1
e
y = f
0
(
x)
y = f (
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2
+
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln
(x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2
+
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing
fence
new
fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallo
w
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f (
x)
11
y = L
(x)
y = f(
x)
F
P
a
a +
∆x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
d
f
∆
x
a
b
y = f(
x)
true
zero
starting
approximation
b
etter approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f(
x)
1
2
y = x
5
0
−
2
y =
1
a
b
y =
sin(x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
=
2
width
of each interval =
2
n
−
2
1
3
0
I
I
I
I
II
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
I
I
I
Ia
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f(
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
a
v
y = f
a
v
c
A
M
0
1
2
a
b
x
t
y = f (
t)
F (
x )
y = f (
t)
F (
x + h)
x + h
F (
x + h) − F (x)
f(
x)
1
2
y =
sin(x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y =
ln|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f(
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
ev
en smaller
y = g(
x)
infinite
area
finite
area
1
y =
1
x
y =
1
x
p
, p
< 1 (typical)
y =
1
x
p
, p
> 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f(
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t =
0
t = π
/6
t = π
/4
t = π
/3
t = π
/2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t =
0
t =
1/2
t =
3/2
t =
2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong
point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r =
3 sin(θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r =
1 + 2 cos(θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r =
1 + cos(θ)
r =
1 +
3
4
cos(
θ)
−
1
4
r =
sin(2θ)
r =
sin(3θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1
+ sin(θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−
4
−
5
4
5
f(
θ)
f(
θ + dθ)
θ
dθ
θ + dθ
appro
ximating region
exact
region
0 ≤ θ ≤ 2
π
r = |
1 + 2 cos(θ)|
2
i
2 − 3
i
−
1
θ =
0
θ =
π
4
θ =
π
2
θ =
2
π
3
θ = π
θ =
13
π
12
θ =
3
π
2
θ =
7
π
4
1
= e
0
e
i
π
4
i = e
i
π
2
e
i
2
π
3
−
1 = e
iπ
e
i
13
π
12
−
i = e
i
3π
2
e
i
7
π
4
i
−
i
1
θ
1 − i
2
i
−
2i
2
−
2
6
i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f(x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
f(x
j+1
)
x
j+1
f(x
0
)
x
0
f(x
1
)
x
1
f(x
2
)
x
2
h
of
Z
2
0
e
−x
2
dx,
and use Simpson’s rule with n = 8. (We can’t do n = 5, since n has to be
even in order to use Simpson’s rule.) The length of each strip is h = (2 −0)/8
units, which is
1
4
, so the partition is
0 <
1
4
<
1
2
<
3
4
< 1 <
5
4
<
3
2
<
7
4
< 2.
By the above formula, we have
Z
2
0
e
−x
2
dx
∼
=
1/4
3
e
−0
2
+ 4e
−(1/4)
2
+ 2e
−(1/2)
2
+ 4e
−(3/4)
2
+ 2e
−1
2
+ 4e
−(5/4)
2
+ 2e
−(3/2)
2
+ 4e
−(7/4)
2
+ e
−2
2
.
This is approximately 0.882066, according to my calculator. This is quite
close to our estimate from the previous section; specifically, when we used the
trapezoidal rule with n = 5, we got the estimate 0.881131. For the record, my
computer program says that the correct value of the integral is 0.882081 to six
decimal places, so Simpson’s rule (with n = 8) is better than the trapezoidal
rule (with n = 5). Of course, a fairer comparison would be to use n = 8 in
both cases; I leave it to you to repeat the trapezoidal rule computation in this
PSfrag
replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shado
w
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror
(y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y =
(x − 1)
2
−
1
x
Same
height
−
x
Same
length,
opp
osite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y =
2
x
y =
10
x
y =
2
−x
y =
log
2
(x)
4
3
units
mirror
(x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
=
0 radians
θ
hyp
otenuse
opp
osite
adjacen
t
0
(≡ 2π)
π
2
π
3
π
2
I
I
I
I
II
IV
θ
(
x, y)
x
y
r
7
π
6
reference
angle
reference
angle =
π
6
sin
+
sin −
cos
+
cos −
tan
+
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this
angle is
5π
6
clo
ckwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
cos(x)
−
π
2
π
2
y =
tan(x), −
π
2
<
x <
π
2
0
−
π
2
π
2
y =
tan(x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y =
sec(x)
y =
csc(x)
y =
cot(x)
y = f(
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y =
tan
−1
(x)
π
2
π
y =
sin(
x)
x
,
x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 <
x < 0.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x +
2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f(
x)
a
tangen
t at x = a
b
tangen
t at x = b
c
tangen
t at x = c
y = x
2
tangen
t
at x = −
1
u
v
uv
u +
∆u
v +
∆v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f(
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f (
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y =
sin(x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin(
x)
x
π
2
π
1
−
1
x =
0
a =
0
x
> 0
a
> 0
x
< 0
a
< 0
rest
position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f(
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y =
ln(x)
y =
cosh(x)
y =
sinh(x)
y =
tanh(x)
y =
sech(x)
y =
csch(x)
y =
coth(x)
1
−
1
y = f(
x)
original
function
in
verse function
slop
e = 0 at (x, y)
slop
e is infinite at (y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
sin(x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y =
sin
−1
(x)
y =
cos(x)
π
π
2
y =
cos
−1
(x)
−
π
2
1
x
α
β
y =
tan(x)
y =
tan(x)
1
y =
tan
−1
(x)
y =
sec(x)
y =
sec
−1
(x)
y =
csc
−1
(x)
y =
cot
−1
(x)
1
y =
cosh
−1
(x)
y =
sinh
−1
(x)
y =
tanh
−1
(x)
y =
sech
−1
(x)
y =
csch
−1
(x)
y =
coth
−1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y =
1
2
(2
, 3)
y = f(
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Lo
cal maximum
Lo
cal minimum
Horizon
tal point of inflection
1
e
y = f
0
(
x)
y = f (
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2
+
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln
(x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2
+
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing
fence
new
fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallo
w
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f (
x)
11
y = L
(x)
y = f(
x)
F
P
a
a +
∆x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
d
f
∆
x
a
b
y = f(
x)
true
zero
starting
approximation
b
etter approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f(
x)
1
2
y = x
5
0
−
2
y =
1
a
b
y =
sin(x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
=
2
width
of each interval =
2
n
−
2
1
3
0
I
I
I
I
II
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
I
I
I
Ia
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f(
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
a
v
y = f
a
v
c
A
M
0
1
2
a
b
x
t
y = f (
t)
F (
x )
y = f (
t)
F (
x + h)
x + h
F (
x + h) − F (x)
f(
x)
1
2
y =
sin(x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y =
ln|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f(
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
ev
en smaller
y = g(
x)
infinite
area
finite
area
1
y =
1
x
y =
1
x
p
, p
< 1 (typical)
y =
1
x
p
, p
> 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f(
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t =
0
t = π
/6
t = π
/4
t = π
/3
t = π
/2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t =
0
t =
1/2
t =
3/2
t =
2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong
point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r =
3 sin(θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r =
1 + 2 cos(θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r =
1 + cos(θ)
r =
1 +
3
4
cos(
θ)
−
1
4
r =
sin(2θ)
r =
sin(3θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1
+ sin(θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−
4
−
5
4
5
f(
θ)
f(
θ + dθ)
θ
dθ
θ + dθ
appro
ximating region
exact
region
0 ≤ θ ≤ 2
π
r = |
1 + 2 cos(θ)|
2
i
2 − 3
i
−
1
θ =
0
θ =
π
4
θ =
π
2
θ =
2π
3
θ = π
θ =
13π
12
θ =
3π
2
θ =
7π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2π
3
−1 = e
iπ
e
i
13π
12
−i = e
i
3π
2
e
i
7π
4
i
−i
1
θ
1 − i
2i
−2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f(x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
f(x
j+1
)
x
j+1
f(x
0
)
x
0
f(x
1
)
x
1
f(x
2
)
x
2
h
case and compare the result with the corresponding Simpson’s rule estimate
which we just found.
B.3.1 Proof of Simpson’s rule
Let’s translate the picture over so that the middle line lies along the y-axis,
like this:
PSfrag
replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shado
w
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror
(y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y =
(x − 1)
2
−
1
x
Same
height
−
x
Same
length,
opp
osite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y =
2
x
y =
10
x
y =
2
−x
y =
log
2
(x)
4
3
units
mirror
(x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
=
0 radians
θ
hyp
otenuse
opp
osite
adjacen
t
0
(≡ 2π)
π
2
π
3
π
2
I
I
I
I
II
IV
θ
(
x, y)
x
y
r
7
π
6
reference
angle
reference
angle =
π
6
sin
+
sin −
cos
+
cos −
tan
+
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this
angle is
5π
6
clo
ckwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
cos(x)
−
π
2
π
2
y =
tan(x), −
π
2
<
x <
π
2
0
−
π
2
π
2
y =
tan(x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y =
sec(x)
y =
csc(x)
y =
cot(x)
y = f (
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y =
tan
−1
(x)
π
2
π
y =
sin(
x)
x
,
x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 <
x < 0.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x +
2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f (
x)
a
tangen
t at x = a
b
tangen
t at x = b
c
tangen
t at x = c
y = x
2
tangen
t
at x = −
1
u
v
uv
u +
∆u
v +
∆v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f (
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f (
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y =
sin(x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin
(x)
x
π
2
π
1
−
1
x =
0
a =
0
x
> 0
a
> 0
x
< 0
a
< 0
rest
position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f(
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y =
ln(x)
y =
cosh(x)
y =
sinh(x)
y =
tanh(x)
y =
sech(x)
y =
csch(x)
y =
coth(x)
1
−
1
y = f (
x)
original
function
in
verse function
slop
e = 0 at (x, y)
slop
e is infinite at (y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
sin(x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y =
sin
−1
(x)
y =
cos(x)
π
π
2
y =
cos
−1
(x)
−
π
2
1
x
α
β
y =
tan(x)
y =
tan(x)
1
y =
tan
−1
(x)
y =
sec(x)
y =
sec
−1
(x)
y =
csc
−1
(x)
y =
cot
−1
(x)
1
y =
cosh
−1
(x)
y =
sinh
−1
(x)
y =
tanh
−1
(x)
y =
sech
−1
(x)
y =
csch
−1
(x)
y =
coth
−1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y =
1
2
(2
, 3)
y = f (
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Lo
cal maximum
Lo
cal minimum
Horizon
tal point of inflection
1
e
y = f
0
(
x)
y = f(
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2
+
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln
(x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2
+
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing
fence
new
fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallo
w
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f (
x)
11
y = L
(x)
y = f (
x)
F
P
a
a +
∆x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
d
f
∆
x
a
b
y = f (
x)
true
zero
starting
approximation
b
etter approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f (
x)
1
2
y = x
5
0
−
2
y =
1
a
b
y =
sin(x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
=
2
width
of each interval =
2
n
−
2
1
3
0
I
I
I
I
II
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
I
I
I
Ia
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f (
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
a
v
y = f
a
v
c
A
M
0
1
2
a
b
x
t
y = f(
t)
F (
x )
y = f(
t)
F (
x + h )
x + h
F (
x + h) − F (x)
f(
x)
1
2
y =
sin(x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y =
ln|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f (
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
ev
en smaller
y = g(
x)
infinite
area
finite
area
1
y =
1
x
y =
1
x
p
, p
< 1 (typical)
y =
1
x
p
, p
> 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f (
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t =
0
t = π
/6
t = π
/4
t = π
/3
t = π
/2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t =
0
t =
1/2
t =
3/2
t =
2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong
point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r =
3 sin(θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r =
1 + 2 cos(θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r =
1 + cos(θ)
r =
1 +
3
4
cos(
θ)
−
1
4
r =
sin(2θ)
r =
sin(3θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1
+ sin(θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−4
−5
4
5
f(θ)
f(θ + dθ)
θ
dθ
θ + dθ
approximating region
exact region
0 ≤ θ ≤ 2π
r = |1 + 2 cos(θ)|
2i
2 − 3i
−1
θ = 0
θ =
π
4
θ =
π
2
θ =
2π
3
θ = π
θ =
13π
12
θ =
3π
2
θ =
7π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2π
3
−1 = e
iπ
e
i
13π
12
−i = e
i
3π
2
e
i
7π
4
i
−i
1
θ
1 − i
2i
−2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f (x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
f(x
j+1
)
x
j+1
f(x
0
)
x
0
f(x
1
)
x
1
f(x
2
)
x
2
h
P
Q
R
h−h
0
As you can see, this shifts the x-coordinates of the partition endpoints over
to −h, 0, and h. Instead of writing f (x
0
), f(x
1
), and f (x
2
), let’s just write
Section B.4: The Error in Our Approximations • 711
P , Q, and R, respectively. The top points are connected by some quadratic,
but we don’t have a clue what it is. Well, let’s call it g and suppose that
g(x) = Ax
2
+ Bx + C. We know that P = g(−h), Q = g(0) and R = g(h);
this means that
P = A(−h)
2
+ B(−h) + C,
Q = A(0)
2
+ B(0) + C,
R = Ah
2
+ Bh + C.
The middle equation says that C = Q; then you can rearrange the other two
equations to see that A = (P + R −2Q)/(2h
2
). (We don’t need to know what
B is!) Now, the shaded area we want is given by
Z
h
−h
(Ax
2
+ Bx + C) dx =
A
3
x
3
+
B
2
x
2
+ Cx
h
−h
=
2Ah
3
3
+ 2Ch
square units, after simplifying. Substituting the values of A and C from above,
the above expression reduces to
2h
3
3
×
P + R − 2Q
2h
2
+ 2Qh =
h
3
(P + 4Q + R).
Now all we have to do is translate to a more general position (which doesn’t
affect the area) and replace P , Q, and R by the function values f(x
0
), f(x
1
),
and f(x
2
), respectively, to get the prototype formula from the beginning of
the previous section.
B.4 The Error in Our Approximations
The whole point of making an approximation (or estimate, if you prefer that
word), is that you end up with something close to the exact quantity you’re
looking for. If you could actually pin your finger on the exact answer, you’d do
that, but sometimes it’s just too difficult. So an approximation at least gives
you a number close to the exact one. As we’ve seen a number of times, most
notably when we looked at linearization and also Taylor series (see Section 13.2
of Chapter 13 and Section 25.3 of Chapter 25), there’s one other important
consideration: how good is the approximation? Is your approximate answer
at least close to the actual one, or does it suck?
To quantify this, we look once again at the error in the approximation,
which is the difference between the actual quantity and the approximation
itself. So suppose that we use one of the above techniques—evenly spaced
strips, the trapezoidal rule or Simpson’s rule—to approximate the integral
R
b
a
f(x) dx. We’d get something like
Z
b
a
f(x) dx = A,
where A is our approximate value. The absolute value of the error is then
given by
|error| =
Z
b
a
f(x) dx − A
.
712 • Estimating Integrals
It turns out that we can get some idea of how big the error could possibly be
by using the derivatives of f, if they exist. In that case, we can let M
1
be the
maximum value of |f
0
(x)| on [a, b]. Similarly, let M
2
be the maximum value
of |f
00
(x)| on [a, b], and finally let M
4
be the maximum value of |f
(4)
(x)| on
[a, b]. Then one can show the following bounds on the error term, depending
on the method used:
for evenly spaced strips, |error| ≤
1
2
M
1
(b − a)h,
for the trapezoidal rule, |error| ≤
1
12
M
2
(b − a)h
2
,
for Simpson’s rule, |error| ≤
1
180
M
4
(b − a)h
4
.
Here h is the strip width (b − a)/n, as usual. Although the above formulas
are pretty similar, there are some differences. First, the coefficients out front
are different. Second, different derivatives are involved: for strips, the first
derivative comes up (in the form of M
1
); for the trapezoidal rule, it’s the
second derivative; while for Simpson’s rule, it’s the fourth derivative. The
most significant difference, however, is the power of h which appears. This
shows how much the error decreases as the strip width becomes smaller, which
of course happens when you take more strips. As h becomes small, h
4
gets
smaller much more quickly than h
2
or h, so Simpson’s rule should kick some
serious butt over the other methods when you use lots of strips.
B.4.1 Examples of estimating the error
Let’s see how these errors turn out for the example
PSfrag
replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shado
w
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror
(y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y =
(x − 1)
2
−
1
x
Same
height
−
x
Same
length,
opp
osite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y =
2
x
y =
10
x
y =
2
−x
y =
log
2
(x)
4
3
units
mirror
(x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
=
0 radians
θ
hyp
otenuse
opp
osite
adjacen
t
0
(≡ 2π)
π
2
π
3
π
2
I
I
I
I
II
IV
θ
(
x, y)
x
y
r
7
π
6
reference
angle
reference
angle =
π
6
sin
+
sin −
cos
+
cos −
tan
+
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this
angle is
5π
6
clo
ckwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
cos(x)
−
π
2
π
2
y =
tan(x), −
π
2
<
x <
π
2
0
−
π
2
π
2
y =
tan(x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y =
sec(x)
y =
csc(x)
y =
cot(x)
y = f(
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y =
tan
−1
(x)
π
2
π
y =
sin(
x)
x
,
x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 <
x < 0.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x +
2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f(
x)
a
tangen
t at x = a
b
tangen
t at x = b
c
tangen
t at x = c
y = x
2
tangen
t
at x = −
1
u
v
uv
u +
∆u
v +
∆v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f(
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f (
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y =
sin(x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin(
x)
x
π
2
π
1
−
1
x =
0
a =
0
x
> 0
a
> 0
x
< 0
a
< 0
rest
position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f(
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y =
ln(x)
y =
cosh(x)
y =
sinh(x)
y =
tanh(x)
y =
sech(x)
y =
csch(x)
y =
coth(x)
1
−
1
y = f(
x)
original
function
in
verse function
slop
e = 0 at (x, y)
slop
e is infinite at (y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
sin(x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y =
sin
−1
(x)
y =
cos(x)
π
π
2
y =
cos
−1
(x)
−
π
2
1
x
α
β
y =
tan(x)
y =
tan(x)
1
y =
tan
−1
(x)
y =
sec(x)
y =
sec
−1
(x)
y =
csc
−1
(x)
y =
cot
−1
(x)
1
y =
cosh
−1
(x)
y =
sinh
−1
(x)
y =
tanh
−1
(x)
y =
sech
−1
(x)
y =
csch
−1
(x)
y =
coth
−1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y =
1
2
(2
, 3)
y = f(
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Lo
cal maximum
Lo
cal minimum
Horizon
tal point of inflection
1
e
y = f
0
(
x)
y = f (
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2
+
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln
(x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2
+
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing
fence
new
fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallo
w
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f (
x)
11
y = L
(x)
y = f(
x)
F
P
a
a +
∆x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
d
f
∆
x
a
b
y = f(
x)
true
zero
starting
approximation
b
etter approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f(
x)
1
2
y = x
5
0
−
2
y =
1
a
b
y =
sin(x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
=
2
width
of each interval =
2
n
−
2
1
3
0
I
I
I
I
II
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
I
I
I
Ia
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f(
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
a
v
y = f
a
v
c
A
M
0
1
2
a
b
x
t
y = f (
t)
F (
x )
y = f (
t)
F (
x + h)
x + h
F (
x + h) − F (x)
f(
x)
1
2
y =
sin(x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y =
ln|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f(
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
ev
en smaller
y = g(
x)
infinite
area
finite
area
1
y =
1
x
y =
1
x
p
, p
< 1 (typical)
y =
1
x
p
, p
> 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f(
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t =
0
t = π
/6
t = π
/4
t = π
/3
t = π
/2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t =
0
t =
1/2
t =
3/2
t =
2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong
point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r =
3 sin(θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r =
1 + 2 cos(θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r =
1 + cos(θ)
r =
1 +
3
4
cos(
θ)
−
1
4
r =
sin(2θ)
r =
sin(3θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1
+ sin(θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−
4
−
5
4
5
f(
θ)
f(
θ + dθ)
θ
dθ
θ + dθ
appro
ximating region
exact
region
0 ≤ θ ≤ 2
π
r = |
1 + 2 cos(θ)|
2
i
2 − 3
i
−
1
θ =
0
θ =
π
4
θ =
π
2
θ =
2π
3
θ = π
θ =
13π
12
θ =
3π
2
θ =
7π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2π
3
−1 = e
iπ
e
i
13π
12
−i = e
i
3π
2
e
i
7π
4
i
−i
1
θ
1 − i
2i
−2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f(x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
f(x
j+1
)
x
j+1
f(x
0
)
x
0
f(x
1
)
x
1
f(x
2
)
x
2
h
P
Q
R
h
−h
0
Z
2
0
e
−x
2
dx,
which we’ve looked at earlier in this appendix. First, we’ll set f (x) = e
−x
2
,
and then calculate that
f
0
(x) = −2xe
−x
2
, f
00
(x) = (4x
2
− 2)e
−x
2
, f
(3)
(x) = −4x(2x
2
− 3)e
−x
2
,
and f
(4)
(x) = 4(4x
4
− 12x
2
+ 3)e
−x
2
.
Let’s find M
1
first. This means that we have to find the maximum value of
|f
0
(x)|, which is actually −f
0
(x), on [0, 2]. Since the second derivative f
00
(x)
is 0 at x = 1/
√
2, and changes sign from negative to positive there, we have
a local minimum for f
0
(x) at 1/
√
2. This means that the minimum value of
f
0
(x) on [0, 2] is −
√
2e
−1/2
, so the maximum value of |f
0
(x)| is
√
2e
−1/2
. That
is, M
1
=
√
2e
−1/2
.
Now we can go back to our estimates for the integral from Section B.1.1
above. There we used 10 evenly spaced strips to estimate our integral. Since
a = 0, b = 2, and h = (2 − 0)/10 =
1
5
, we have
|error using 10 strips| ≤
1
2
M
1
(b − a)h =
1
2
×
√
2e
−1/2
(2 − 0)
1
5
=
√
2
5
e
−1/2
.
Section B.4.1: Examples of estimating the error • 713
This is approximately 0.171553. Note that it doesn’t matter if you use the
left-hand endpoint, the right-hand endpoint, or somewhere in between as your
choice of c
n
. (In Section B.1.1, we used the right-hand endpoints and then
the midpoints to get two different estimates, but they are both accurate to
about ±0.171553.)
Let’s move on to the trapezoidal rule estimate. In Section B.2 above, we
used 5 trapezoids of width h = 2/5 to estimate our integral (so n = 5). To see
how big the error could be, we’ll need to find M
2
by maximizing the value of
|f
00
(x)| on [0, 2]. To do this, look back at the above formulas for f
(2)
(x) and
f
(3)
(x). The zeroes of f
(3)
(x) which lie in [0, 2] are at x = 0 and x =
p
3/2,
so these are the critical points of f
(2)
(x). (Remember, the third derivative is
the derivative of the second derivative!) So we can test the values of f
00
(0)
and f
00
(
p
3/2), as well as the value f(2) at the other endpoint 2. We find
f
00
(0) = −2, f
00
(
p
3/2) = 4e
−3/2
and f
00
(2) = 14e
−4
. The largest of these, in
absolute value, is f
00
(0). This means that M
2
= 2. Now we can estimate the
error (remembering that h = 2/5):
|error using 5 trapezoids| ≤
1
12
M
2
(b − a)h
2
=
1
12
× 2(2 − 0)
2
5
2
=
4
75
,
which is 0.053333 . . .. This is a lot less than the error using 10 strips, even
though we only used 5 trapezoids! Since our previous estimate was approxi-
mately 0.881131, we have shown that the approximation
Z
2
0
e
−x
2
dx
∼
=
0.881131
is accurate to about ±0.053333. (This is certainly consistent with my ob-
servation at the end of Section B.3 above that the correct value is actually
0.882081, to six decimal places.)
Finally, we’ll estimate the error using Simpson’s rule. In Section B.3 above,
we used Simpson’s rule with n = 8 to show that
Z
2
0
e
−x
2
dx
∼
=
0.882066.
We’ll need M
4
, which is the maximum value of |f
(4)
(x)| on [0, 2]. This could
be very messy, since f
(4)
(x) = 4(4x
4
− 12x
2
+ 3)e
−x
2
. Let’s cheat by finding
the maximum value of each of the three factors. No problem with 4, and e
−x
2
is positive and is maximized at x = 0 (with a value of 1); so we only have to
find the maximum of |4x
4
− 12x
2
+ 3| on [0, 2]. We have
d
dx
(4x
4
− 12x
2
+ 3) = 16x
3
− 24x = 8x(2x
2
− 3),
so the maximum we’re looking for could only occur at one of the critical points
x = 0 and x =
p
3/2, or at the other endpoint x = 2. Plugging these numbers
in, we find that the greatest value 19 occurs at x = 2, which means that
|4x
4
− 12x
2
+ 3| ≤ 19
714 • Estimating Integrals
on [0, 2]. Putting everything together, we can say that
M
4
≤ 4 × 19 × 1 = 76.
(Actually, M
4
= 12, but you need to look at the fifth derivative of f to
see this, and enough is enough!) Now we can finally use our formula (with
h = (2 − 0)/8 = 1/4):
|error using Simpson’s rule with n = 8| ≤
1
180
M
4
(b − a)h
4
≤
1
180
× 76(2 − 0)
1
4
4
=
19
5760
.
This is about 0.003299, which is much lower than the previous two errors we
calculated.
B.4.2 Proof of an error term inequality
The proofs of the last two of the three error inequalities in Section B.4 above
are a little beyond the scope of this book, but it’s not too hard to show that
the first one is true:
|error using n evenly spaced strips of width h| ≤
1
2
M
1
(b − a)h,
where M
1
is the maximum value of f
0
(x) in [a, b]. Suppose that we use the
left-hand endpoints for our estimate. Let’s look at just one of the strips. If its
base is the interval [q, q + h] (for some q), then it looks something like this:
PSfrag
replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shado
w
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror
(y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y =
(x − 1)
2
−
1
x
Same
height
−
x
Same
length,
opp
osite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y =
2
x
y =
10
x
y =
2
−x
y =
log
2
(x)
4
3
units
mirror
(x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
=
0 radians
θ
hypotenuse
opp
osite
adjacen
t
0
(≡ 2π)
π
2
π
3
π
2
I
I
I
I
II
IV
θ
(
x, y)
x
y
r
7
π
6
reference
angle
reference
angle =
π
6
sin
+
sin −
cos
+
cos −
tan
+
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this
angle is
5π
6
clo
ckwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
cos(x)
−
π
2
π
2
y =
tan(x), −
π
2
<
x <
π
2
0
−
π
2
π
2
y =
tan(x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y =
sec(x)
y =
csc(x)
y =
cot(x)
y = f (
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y =
tan
−1
(x)
π
2
π
y =
sin(
x)
x
,
x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 <
x < 0.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x +
2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f (z))
z
y = f (
x)
a
tangen
t at x = a
b
tangen
t at x = b
c
tangen
t at x = c
y = x
2
tangen
t
at x = −
1
u
v
uv
u +
∆u
v +
∆v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f (
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f (
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y =
sin(x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin
(x)
x
π
2
π
1
−
1
x =
0
a =
0
x
> 0
a
> 0
x
< 0
a
< 0
rest
position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f (
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y =
ln(x)
y =
cosh(x)
y =
sinh(x)
y =
tanh(x)
y =
sech(x)
y =
csch(x)
y =
coth(x)
1
−
1
y = f (
x)
original
function
in
verse function
slop
e = 0 at (x, y)
slop
e is infinite at (y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
sin(x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y =
sin
−1
(x)
y =
cos(x)
π
π
2
y =
cos
−1
(x)
−
π
2
1
x
α
β
y =
tan(x)
y =
tan(x)
1
y =
tan
−1
(x)
y =
sec(x)
y =
sec
−1
(x)
y =
csc
−1
(x)
y =
cot
−1
(x)
1
y =
cosh
−1
(x)
y =
sinh
−1
(x)
y =
tanh
−1
(x)
y =
sech
−1
(x)
y =
csch
−1
(x)
y =
coth
−1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y =
1
2
(2
, 3)
y = f (
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Lo
cal maximum
Lo
cal minimum
Horizon
tal point of inflection
1
e
y = f
0
(
x)
y = f (
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2
+
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln
(x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2
+
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing
fence
new
fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallo
w
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f (
x)
11
y = L
(x)
y = f (
x)
F
P
a
a +
∆x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
d
f
∆
x
a
b
y = f (
x)
true
zero
starting
approximation
b
etter approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f (
x)
1
2
y = x
5
0
−
2
y =
1
a
b
y =
sin(x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
=
2
width
of each interval =
2
n
−
2
1
3
0
I
I
I
I
II
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
I
I
I
Ia
5
3
0
1
2
a
b
y = f(
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f (
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
a
v
y = f
a
v
c
A
M
0
1
2
a
b
x
t
y = f(
t)
F (
x )
y = f(
t)
F (
x + h )
x + h
F (
x + h) − F (x)
f(
x)
1
2
y =
sin(x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y =
ln|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f (
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
ev
en smaller
y = g(
x)
infinite
area
finite
area
1
y =
1
x
y =
1
x
p
, p
< 1 (typical)
y =
1
x
p
, p
> 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f (
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t =
0
t = π
/6
t = π
/4
t = π
/3
t = π
/2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t =
0
t =
1/2
t =
3/2
t =
2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong
point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r =
3 sin(θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r =
1 + 2 cos(θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r =
1 + cos(θ)
r = 1 +
3
4
cos(θ)
−
1
4
r = sin(2θ)
r = sin(3θ)
r =
1
π
θ
0 ≤ θ ≤ 4π
r =
2
1 + sin(θ)
−
π
4
≤ θ ≤
5π
4
0 ≤ θ ≤ 2π
0 ≤ θ ≤ π
−4
−5
4
5
f(θ)
f(θ + dθ)
θ
dθ
θ + dθ
approximating region
exact region
0 ≤ θ ≤ 2π
r = |1 + 2 cos(θ)|
2i
2 − 3i
−1
θ = 0
θ =
π
4
θ =
π
2
θ =
2π
3
θ = π
θ =
13π
12
θ =
3π
2
θ =
7π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2π
3
−1 = e
iπ
e
i
13π
12
−i = e
i
3π
2
e
i
7π
4
i
−i
1
θ
1 − i
2i
−2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f (x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
f(x
j+1
)
x
j+1
f(x
0
)
x
0
f(x
1
)
x
1
f(x
2
)
x
2
h
P
Q
R
h
−h
0
q q + h
f(q)
y = f (x)
The approximating rectangle has height f(q) and width h units, so the ap-
proximate area is hf(q) square units. How bad could this be, in general?
It all depends on how much the graph of f deviates from the constant line
y = f (q). Here are the two worst-case scenarios:
Section B.4.2: Proof of an error term inequality • 715
PSfrag
replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shado
w
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror
(y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y =
(x − 1)
2
−
1
x
Same
height
−
x
Same
length,
opp
osite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y =
2
x
y =
10
x
y =
2
−x
y =
log
2
(x)
4
3
units
mirror
(x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
=
0 radians
θ
hypotenuse
opp
osite
adjacen
t
0
(≡ 2π)
π
2
π
3
π
2
I
I
I
I
II
IV
θ
(
x, y)
x
y
r
7
π
6
reference
angle
reference
angle =
π
6
sin
+
sin −
cos
+
cos −
tan
+
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this
angle is
5π
6
clo
ckwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
cos(x)
−
π
2
π
2
y =
tan(x), −
π
2
<
x <
π
2
0
−
π
2
π
2
y =
tan(x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y =
sec(x)
y =
csc(x)
y =
cot(x)
y = f (
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y =
tan
−1
(x)
π
2
π
y =
sin(
x)
x
,
x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 <
x < 0.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x +
2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f (z))
z
y = f (
x)
a
tangen
t at x = a
b
tangen
t at x = b
c
tangen
t at x = c
y = x
2
tangen
t
at x = −
1
u
v
uv
u +
∆u
v +
∆v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f (
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f (
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y =
sin(x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin
(x)
x
π
2
π
1
−
1
x =
0
a =
0
x
> 0
a
> 0
x
< 0
a
< 0
rest
position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f (
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y =
ln(x)
y =
cosh(x)
y =
sinh(x)
y =
tanh(x)
y =
sech(x)
y =
csch(x)
y =
coth(x)
1
−
1
y = f (
x)
original
function
in
verse function
slop
e = 0 at (x, y)
slop
e is infinite at (y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
sin(x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y =
sin
−1
(x)
y =
cos(x)
π
π
2
y =
cos
−1
(x)
−
π
2
1
x
α
β
y =
tan(x)
y =
tan(x)
1
y =
tan
−1
(x)
y =
sec(x)
y =
sec
−1
(x)
y =
csc
−1
(x)
y =
cot
−1
(x)
1
y =
cosh
−1
(x)
y =
sinh
−1
(x)
y =
tanh
−1
(x)
y =
sech
−1
(x)
y =
csch
−1
(x)
y =
coth
−1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y =
1
2
(2
, 3)
y = f (
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Lo
cal maximum
Lo
cal minimum
Horizon
tal point of inflection
1
e
y = f
0
(
x)
y = f (
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2
+
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln
(x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2
+
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing
fence
new
fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallo
w
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f (
x)
11
y = L
(x)
y = f (
x)
F
P
a
a +
∆x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
d
f
∆
x
a
b
y = f (
x)
true
zero
starting
approximation
b
etter approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f (
x)
1
2
y = x
5
0
−
2
y =
1
a
b
y =
sin(x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
=
2
width
of each interval =
2
n
−
2
1
3
0
I
I
I
I
II
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
I
I
I
Ia
5
3
0
1
2
a
b
y = f(
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f (
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
a
v
y = f
a
v
c
A
M
0
1
2
a
b
x
t
y = f(
t)
F (
x )
y = f(
t)
F (
x + h )
x + h
F (
x + h) − F (x)
f(
x)
1
2
y =
sin(x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y =
ln|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f (
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
ev
en smaller
y = g(
x)
infinite
area
finite
area
1
y =
1
x
y =
1
x
p
, p
< 1 (typical)
y =
1
x
p
, p
> 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f (
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t =
0
t = π
/6
t = π
/4
t = π
/3
t = π
/2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t =
0
t =
1/2
t =
3/2
t =
2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong
point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r =
3 sin(θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r =
1 + 2 cos(θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r =
1 + cos(θ)
r =
1 +
3
4
cos(
θ)
−
1
4
r =
sin(2θ)
r =
sin(3θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1
+ sin(θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−
4
−
5
4
5
f(
θ)
f(
θ + dθ)
θ
dθ
θ + dθ
appro
ximating region
exact
region
0 ≤ θ ≤ 2
π
r = |
1 + 2 cos(θ)|
2
i
2 − 3
i
−
1
θ =
0
θ =
π
4
θ =
π
2
θ =
2π
3
θ = π
θ =
13π
12
θ =
3π
2
θ =
7π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2π
3
−1 = e
iπ
e
i
13π
12
−i = e
i
3π
2
e
i
7π
4
i
−i
1
θ
1 − i
2i
−2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f (x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
f(x
j+1
)
x
j+1
f(x
0
)
x
0
f(x
1
)
x
1
f(x
2
)
x
2
h
P
Q
R
h
h
−h
0
qq q + hq + h
f(q)f(q)
y = f (x)
slope M
1
slope −M
1
M
1
h
M
1
h
The first shows a line segment of slope M
1
beginning at (q, f(q)), while the
second shows a line segment of slope −M
1
beginning at the same point. In
fact, the function must be trapped between these two extremes. Indeed, the
first line has equation y = f(q) + M
1
(x − q). If f(x) ever gets above this line
(for x in the interval [q, q + h]), then we have f(x) > f(q) + M
1
(x − q), or
f(x) − f (q)
x − q
> M
1
.
By the Mean Value Theorem (see Section 11.3 in Chapter 11), the left-hand
side is equal to f
0
(c) for some c in [q, x], so f
0
(c) > M
1
. This is impossible,
since M
1
is the maximum value of |f
0
(x)| on [a, b]. A similar argument shows
that y = f(x) always lies above the down-sloping line.
Now we can look at the error. In the first worst-case scenario, the actual
region includes the strip plus a triangle of side lengths h by M
1
h units; in
the second worst-case scenario, the same triangle is actually removed from
the strip. In either case, our area could be off by the area of the triangle,
which is
1
2
M
1
h
2
square units. All that’s left is to multiply this error by the
number of strips, which is n, to see that our approximation can’t be worse
than
1
2
M
1
h
2
n. In fact, we can steal away one of the factors of h and use the
equation nh = (b − a) to rewrite the above expression as
1
2
M
1
(b − a)h. This
is exactly what we want! Now I leave it to you to repeat the above argument
PSfrag
replacements
(
a, b)
[
a, b]
(
a, b]
[
a, b)
(
a, ∞)
[
a, ∞)
(
−∞, b)
(
−∞, b]
(
−∞, ∞)
{
x : a < x < b}
{
x : a ≤ x ≤ b}
{
x : a < x ≤ b}
{
x : a ≤ x < b}
{
x : x ≥ a}
{
x : x > a}
{
x : x ≤ b}
{
x : x < b}
R
a
b
shado
w
0
1
4
−
2
3
−
3
g(
x) = x
2
f(
x) = x
3
g(
x) = x
2
f(
x) = x
3
mirror
(y = x)
f
−
1
(x) =
3
√
x
y = h
(x)
y = h
−
1
(x)
y =
(x − 1)
2
−
1
x
Same
height
−
x
Same
length,
opp
osite signs
y = −
2x
−
2
1
y =
1
2
x − 1
2
−
1
y =
2
x
y =
10
x
y =
2
−x
y =
log
2
(x)
4
3
units
mirror
(x-axis)
y = |
x|
y = |
log
2
(x)|
θ radians
θ units
30
◦
=
π
6
45
◦
=
π
4
60
◦
=
π
3
120
◦
=
2
π
3
135
◦
=
3
π
4
150
◦
=
5
π
6
90
◦
=
π
2
180
◦
= π
210
◦
=
7
π
6
225
◦
=
5
π
4
240
◦
=
4
π
3
270
◦
=
3
π
2
300
◦
=
5
π
3
315
◦
=
7
π
4
330
◦
=
11
π
6
0
◦
=
0 radians
θ
hyp
otenuse
opp
osite
adjacen
t
0
(≡ 2π)
π
2
π
3
π
2
I
I
I
I
II
IV
θ
(
x, y)
x
y
r
7
π
6
reference
angle
reference
angle =
π
6
sin
+
sin −
cos
+
cos −
tan
+
tan −
A
S
T
C
7
π
4
9
π
13
5
π
6
(this
angle is
5π
6
clo
ckwise)
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
cos(x)
−
π
2
π
2
y =
tan(x), −
π
2
<
x <
π
2
0
−
π
2
π
2
y =
tan(x)
−
2π
−
3π
−
5
π
2
−
3
π
2
−
π
−
π
2
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
y =
sec(x)
y =
csc(x)
y =
cot(x)
y = f(
x)
−
1
1
2
y = g(
x)
3
y = h
(x)
4
5
−
2
f(
x) =
1
x
g(
x) =
1
x
2
etc.
0
1
π
1
2
π
1
3
π
1
4
π
1
5
π
1
6
π
1
7
π
g(
x) = sin
1
x
1
0
−
1
L
10
100
200
y =
π
2
y = −
π
2
y =
tan
−1
(x)
π
2
π
y =
sin(
x)
x
,
x > 3
0
1
−
1
a
L
f(
x) = x sin (1/x)
(0 <
x < 0.3)
h
(x) = x
g(
x) = −x
a
L
lim
x
→a
+
f(x) = L
lim
x
→a
+
f(x) = ∞
lim
x
→a
+
f(x) = −∞
lim
x
→a
+
f(x) DNE
lim
x
→a
−
f(x) = L
lim
x
→a
−
f(x) = ∞
lim
x
→a
−
f(x) = −∞
lim
x
→a
−
f(x) DNE
M
}
lim
x
→a
−
f(x) = M
lim
x
→a
f(x) = L
lim
x
→a
f(x) DNE
lim
x
→∞
f(x) = L
lim
x
→∞
f(x) = ∞
lim
x
→∞
f(x) = −∞
lim
x
→∞
f(x) DNE
lim
x
→−∞
f(x) = L
lim
x
→−∞
f(x) = ∞
lim
x
→−∞
f(x) = −∞
lim
x
→−∞
f(x) DNE
lim
x →a
+
f(
x) = ∞
lim
x →a
+
f(
x) = −∞
lim
x →a
−
f(
x) = ∞
lim
x →a
−
f(
x) = −∞
lim
x →a
f(
x) = ∞
lim
x →a
f(
x) = −∞
lim
x →a
f(
x) DNE
y = f (
x)
a
y =
|
x|
x
1
−
1
y =
|
x + 2|
x +
2
1
−
1
−
2
1
2
3
4
a
a
b
y = x sin
1
x
y = x
y = −
x
a
b
c
d
C
a
b
c
d
−
1
0
1
2
3
time
y
t
u
(
t, f(t))
(
u, f(u))
time
y
t
u
y
x
(
x, f(x))
y = |
x|
(
z, f(z))
z
y = f(
x)
a
tangen
t at x = a
b
tangen
t at x = b
c
tangen
t at x = c
y = x
2
tangen
t
at x = −
1
u
v
uv
u +
∆u
v +
∆v
(
u + ∆u)(v + ∆v)
∆
u
∆
v
u
∆v
v∆
u
∆
u∆v
y = f(
x)
1
2
−
2
y = |
x
2
− 4|
y = x
2
− 4
y = −
2x + 5
y = g(
x)
1
2
3
4
5
6
7
8
9
0
−
1
−
2
−
3
−
4
−
5
−
6
y = f (
x)
3
−
3
3
−
3
0
−
1
2
easy
hard
flat
y = f
0
(
x)
3
−
3
0
−
1
2
1
−
1
y =
sin(x)
y = x
x
A
B
O
1
C
D
sin(
x)
tan(
x)
y =
sin(
x)
x
π
2
π
1
−
1
x =
0
a =
0
x
> 0
a
> 0
x
< 0
a
< 0
rest
position
+
−
y = x
2
sin
1
x
N
A
B
H
a
b
c
O
H
A
B
C
D
h
r
R
θ
1000
2000
α
β
p
h
y = g(
x) = log
b
(x)
y = f(
x) = b
x
y = e
x
5
10
1
2
3
4
0
−
1
−
2
−
3
−
4
y =
ln(x)
y =
cosh(x)
y =
sinh(x)
y =
tanh(x)
y =
sech(x)
y =
csch(x)
y =
coth(x)
1
−
1
y = f(
x)
original
function
in
verse function
slop
e = 0 at (x, y)
slop
e is infinite at (y, x)
−
108
2
5
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
y =
sin(x)
y =
sin(x), −
π
2
≤ x ≤
π
2
−
2
−
1
0
2
π
2
−
π
2
y =
sin
−1
(x)
y =
cos(x)
π
π
2
y =
cos
−1
(x)
−
π
2
1
x
α
β
y =
tan(x)
y =
tan(x)
1
y =
tan
−1
(x)
y =
sec(x)
y =
sec
−1
(x)
y =
csc
−1
(x)
y =
cot
−1
(x)
1
y =
cosh
−1
(x)
y =
sinh
−1
(x)
y =
tanh
−1
(x)
y =
sech
−1
(x)
y =
csch
−1
(x)
y =
coth
−1
(x)
(0
, 3)
(2
, −1)
(5
, 2)
(7
, 0)
(
−1, 44)
(0
, 1)
(1
, −12)
(2
, 305)
y =
1
2
(2
, 3)
y = f(
x)
y = g(
x)
a
b
c
a
b
c
s
c
0
c
1
(
a, f(a))
(
b, f(b))
1
2
1
2
3
4
5
6
0
−
1
−
2
−
3
−
4
−
5
−
6
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
3
π
5
π
2
2
π
3
π
2
π
π
2
y =
sin(x)
1
0
−
1
−
3π
−
5
π
2
−
2π
−
3
π
2
−
π
−
π
2
3
π
5
π
2
2
π
2
π
3
π
2
π
π
2
c
OR
Lo
cal maximum
Lo
cal minimum
Horizon
tal point of inflection
1
e
y = f
0
(
x)
y = f (
x) = x ln(x)
−
1
e
?
y = f(
x) = x
3
y = g(
x) = x
4
x
f(
x)
−
3
−
2
−
1
0
1
2
1
2
3
4
+
−
?
1
5
6
3
f
0
(
x)
2 −
1
2
√
6
2
+
1
2
√
6
f
00
(
x)
7
8
g
00
(
x)
f
00
(
x)
0
y =
(
x − 3)(x − 1)
2
x
3
(
x + 2)
y = x ln
(x)
1
e
−
1
e
5
−
108
2
α
β
2 −
1
2
√
6
2
+
1
2
√
6
y = x
2
(
x − 5)
3
−
e
−
1/2
√
3
e
−
1/2
√
3
−
e
−3/2
e
−
3/2
−
1
√
3
1
√
3
−
1
1
y = xe
−
3x
2
/2
y =
x
3
− 6
x
2
+ 13x − 8
x
28
2
600
500
400
300
200
100
0
−
100
−
200
−
300
−
400
−
500
−
600
0
10
−
10
5
−
5
20
−
20
15
−
15
0
4
5
6
x
P
0
(
x)
+
−
−
existing
fence
new
fence
enclosure
A
h
b
H
99
100
101
h
dA/dh
r
h
1
2
7
shallo
w
deep
LAND
SEA
N
y
z
s
t
3
11
9
L
(11)
√
11
y = L
(x)
y = f (
x)
11
y = L
(x)
y = f(
x)
F
P
a
a +
∆x
f(
a + ∆x)
L
(a + ∆x)
f(
a)
error
d
f
∆
x
a
b
y = f(
x)
true
zero
starting
approximation
b
etter approximation
v
t
3
5
50
40
60
4
20
30
25
t
1
t
2
t
3
t
4
t
n
−2
t
n
−1
t
0
= a
t
n
= b
v
1
v
2
v
3
v
4
v
n
−1
v
n
−
30
6
30
|
v|
a
b
p
q
c
v(
c)
v(
c
1
)
v(
c
2
)
v(
c
3
)
v(
c
4
)
v(
c
5
)
v(
c
6
)
t
1
t
2
t
3
t
4
t
5
c
1
c
2
c
3
c
4
c
5
c
6
t
0
=
a
t
6
=
b
t
16
=
b
t
10
=
b
a
b
x
y
y = f(
x)
1
2
y = x
5
0
−
2
y =
1
a
b
y =
sin(x)
π
−
π
0
−
1
−
2
0
2
4
y = x
2
0
1
2
3
4
2
n
4
n
6
n
2(
n−2)
n
2(
n−1)
n
2
n
n
=
2
width
of each interval =
2
n
−
2
1
3
0
I
I
I
I
II
IV
4
y
dx
y = −
x
2
− 2x + 3
3
−
5
y = |−
x
2
− 2x + 3|
I
I
I
I
Ia
5
3
0
1
2
a
b
y = f (
x)
y = g(
x)
y = x
2
a
b
5
3
0
1
2
y =
√
x
2
√
2
2
2
dy
x
2
a
b
y = f(
x)
y = g(
x)
M
m
1
2
−
1
−
2
0
y = e
−
x
2
1
2
e
−
1/4
f
a
v
y = f
a
v
c
A
M
0
1
2
a
b
x
t
y = f (
t)
F (
x )
y = f (
t)
F (
x + h)
x + h
F (
x + h) − F (x)
f(
x)
1
2
y =
sin(x)
π
−
π
−
1
−
2
y =
1
x
y = x
2
1
2
1
−
1
y =
ln|x|
θ
a
x
a
x
p
a
2
− x
2
3
x
p
9 − x
2
p
x
2
+ a
2
x
a
p
x
2
+ 15
x
√
15
x
p
x
2
− a
2
a
x
p
x
2
− 4
2
x
−
p
x
2
− a
2
a
x
−
p
x
2
− 4
2
y = f(
x)
a
b
a + ε
ε
Z
b
a
+ε
f(x) dx
small
ev
en smaller
y = g(
x)
infinite
area
finite
area
1
y =
1
x
y =
1
x
p
, p
< 1 (typical)
y =
1
x
p
, p
> 1 (typical)
a
1
a
2
a
3
a
4
a
5
a
6
a
7
a
8
1
2
3
4
5
6
7
8
n
a
n
x
y
y = f(
x)
(
a, f(a))
a
−
1
0
1
a
6
1
2
7
1
2
7
?
−
2
−
1
−
2
t =
0
t = π
/6
t = π
/4
t = π
/3
t = π
/2
3
0
t = −
2
t = −
3/2
t = ±
1
t = −
1/2
t =
0
t =
1/2
t =
3/2
t =
2
12
−
12
θ
r
P
θ
r
P
11
π
6
2
(
−1, −1)
wrong
point
π
4
5
π
4
√
2
(0
, 1)
(0
, −3)
(
−2, 0)
π
2
3
π
2
π
r =
3 sin(θ)
3
π
2
θ
2
π
1
0
−
1
−
2
−
3
0
3
2
−
3
2
0
r =
1 + 2 cos(θ)
2
π
3
4
π
3
0
π
0
pi
−
3
2
3
π
2
1
2
3
0
−
1
−
2
−
3
0 ≤ θ ≤
2
π
3
0 ≤ θ ≤ π
0 ≤ θ ≤ 2
π
r =
1 + cos(θ)
r =
1 +
3
4
cos(
θ)
−
1
4
r =
sin(2θ)
r =
sin(3θ)
r =
1
π
θ
0 ≤ θ ≤ 4
π
r =
2
1
+ sin(θ)
−
π
4
≤ θ ≤
5
π
4
0 ≤ θ ≤ 2
π
0 ≤ θ ≤ π
−
4
−
5
4
5
f(
θ)
f(
θ + dθ)
θ
dθ
θ + dθ
approximating region
exact region
0 ≤ θ ≤ 2π
r = |1 + 2 cos(θ)|
2i
2 − 3i
−1
θ = 0
θ =
π
4
θ =
π
2
θ =
2π
3
θ = π
θ =
13π
12
θ =
3π
2
θ =
7π
4
1 = e
0
e
i
π
4
i = e
i
π
2
e
i
2π
3
−1 = e
iπ
e
i
13π
12
−i = e
i
3π
2
e
i
7π
4
i
−i
1
θ
1 − i
2i
−2i
2
−2
6i
−6i
6
−6
−
√
3
R
ϕ
2
1/5
θ =
π
6
θ =
17π
30
θ =
29π
30
θ =
41π
30
θ =
53π
30
z
0
z
1
z
2
z
3
z
4
−
√
3
2
√
3
2
1
2
i
−i
19π
6
−i
7π
6
i
5π
6
i
17π
6
i
29π
6
ln(2)
−
7π
4
−
3π
4
π
4
5π
4
9π
4
3
2
i
0
1
2
3
4
dx
y
x
y =
p
1 − (x − 3)
2
2πx
a
b
y = f(x)
A
B
y =
√
x
1
y = 2x
3
y = x
4
(2, 16)
−5
5
6
y = h
y−h
h−y
x = h
y
x−h
radius of shell = x−h
h−x
radius of shell = h−x
8
P
h
P
(slice)
(axis)
l
L
1
2
Base
Cross-section
Area = A
Area = A(x)
y = e
x
A
B
dx
dy
x + dx
a
b
p
(dx)
2
+ (dy)
2
P
t
600000
500000
400000
300000
200000
100000
−100000
−200000
−300000
−400000
−500000
−600000
0
1
2
3
4
5
6
7
L
ε
L + ε
L − ε
a
M
a + δ
a − δ
your move
my move
N
M
a + δ
a − δ
your move
my move
FIRST MOVE
SECOND MOVE
L + ε
L − ε
1
−1
L −
1
2
L +
1
2
L = 0
L = −
1
4
L =
7
8
L = 0
b
each segment is either the right half
or the left half of the one below it
infinitely many marked points
lie below each segment
0
1
2
−1
−2
y = e
−x
2
0
1
2
−1
−2
f(x
j
)
f(x
j−1
)
x
j
x
j−1
h
f(x
j+1
)
x
j+1
f(x
0
)
x
0
f(x
1
)
x
1
f(x
2
)
x
2
h
P
Q
R
h
−h
0
q
q + h
f(q)
y = f(x)
slope M
1
slope −M
1
M
1
h
in the case where we don’t necessarily choose the left-hand endpoint. (In
fact, if you use the midpoint, you can show that the error is actually only
1
4
M
1
(b − a)h.)