
Section 29.1.2: The shell method • 621
This weird solid is a pretty lumpy bagel half, but its volume is fairly close to
what we’re looking for. The thinner the maximum shell thickness, the better
the approximation. As before, integrating takes care of both the addition of
all the shell volumes and also taking the limit as the maximum shell thickness
goes to zero.
First we need to find the volume of one generic shell. The easiest way to
do this is to think of the shell as a really thin metal can without a top or
bottom. As you can see from the picture of the shell on the previous page,
the height of the can is y units, the radius is x units, and the thickness is dx
units. Imagine cutting the can down the side with some sharp scissors, then
unfolding it and flattening it out into a thin rectangle-like piece of metal. It’s
not actually a rectangle, of course. You see, a rectangle is a 2-dimensional
object, whereas the unrolled can is 3-dimensional—although the can is pretty
thin, it still has some thickness. (Even a piece of paper has some thickness,
or else a ream of paper would be really really thin.) Now it’s actually not
even a rectangular prism, since the inner radius of the can isn’t exactly the
same as the outer radius. But the point is, it’s almost a rectangular prism.
The thinner the can gets, the closer it is to a rectangular prism, and when we
take limits in the end (using the integral), everything will work out.
∗
So, the
idealized version of the unfolded can looks 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
2πx
The thickness is dx units, and the side we cut along is still the height of the
cylindrical shell, that is, y units. How about the long side? Well, that is equal
to the circumference of the shell (think about it!) which is 2πx units, since
the radius of the shell is basically x units. So, the volume of the shell is very
close to 2πxy dx cubic units. Now all we have to do is integrate from x = 2
to x = 4 to see that the volume of the bagel half (in cubic units) is
Z
4
2
2πxy dx = 2π
Z
4
2
x
p
1 − (x − 3)
2
dx.
Great—we’ve now reduced the problem to evaluating a definite integral, but
it’s a bit of a messy one. Start off by substituting t = x −3, so dt = dx; also,
when x = 2, we have t = −1, and when x = 4, we see that t = 1. So in t-land,
the integral becomes
2π
Z
1
−1
(t + 3)
p
1 − t
2
dt = 2π
Z
1
−1
t
p
1 − t
2
dt + 3
Z
1
−1
p
1 − t
2
dt
.
∗
More formally, we can view the volume of the shell as the difference in volumes of
the outer shell (of radius x + dx units) and the inner shell (of radius x units). Both shells
have height y units, so the volume of the shell is πy((x + dx)
2
− x
2
), which simplifies to
2πxy dx + πy(dx)
2
cubic units. When this is integrated, the second term vanishes due to
the negligible quantity (dx)
2
.