CoCalc Public Files05 - Volume part 3 Assignment / Volume part 3 Notes.sagewsOpen with one click!
Author: Aaron Tresham
Views : 48
License: Other -- explicitly state in your code
Compute Environment: Ubuntu 20.04 (Default)
This material was developed by Aaron Tresham at the University of Hawaii at Hilo and is Creative Commons License
licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.


  • Intro to Sage
  • Graphing and Solving Equations
  • Symbolic Integration

Volume, part 3: Cylindrical Shells

For some solids of revolution, the disk and washer methods we discussed last week are difficult, if not impossible. For these examples, there is another method: cylindrical shells.

Example 1

Find the volume of the solid obtained by rotating around the y-axis the region between f(x)=2x2x3f(x)=2x^2-x^3 and the x-axis from x=0x=0 to x=2x=2.

Here is the region:

f(x)=2*x^2-x^3 plot(f,xmin=0,xmax=2,fill='axis')

Here is the surface we get by rotating the curve. It's kind of like a bundt cake.

If we take cross sections perpendicular to the y-axis (the axis of rotation), we get annuli (or washers). But we run into problems with this method:

  • The inner and outer radii are given by the same function.

  • To find the limits of integration, we need to find the maximum y-value on this interval.

  • We need to know xx given yy; that is, we must solve y=2x2x3y=2x^2-x^3 for xx.

For this particular example, we could solve these problems (but it would take some work). Let's try another approach, which turns out to be easier.

Instead of slicing with a horizontal knife (to get cross sections), let's cut this "cake" using cylindrical "cookie cutters". That is, take a cylinder of radius rr (centered on the axis of rotation) and press it straight through the solid. Then take a cylinder with a slightly bigger radius (say RR) and press it down. What you have now cut out is a cylindrical shell, a slab of cake shaped like a circular cylinder.

Here's a sample picture I found online (it's from a YouTube video by Andrea Marchese of the City College of New York):

What's the volume of this cylinder of cake?

Slice the slab down the side and unroll it. What you have is a flat slab of cake. Its height is the height of the cake at that point (f(R)f(r))(f(R)\approx f(r)), its width is the circumference of the cylinder (2πR2πr)(2\pi R\approx 2\pi r), and its thickness is RrR-r. I have the approximations because our cylindrical shell is not exactly flat on top (because of the curve of the function).

So the volume of this slab is approximately f(r)2πr(Rr)f(r)\cdot2\pi r\cdot (R-r).

Now we need to keep cutting more slabs until we've cut up the whole cake, approximate the volume of each slab, and add them all up. That gives the approximate volume of the cake.

Of course, we make the approximation better by cutting thinner slabs. In other words, we let the number of cylindrical shells approach infinity.

Does this process sound familiar? It should; this is the limit of Riemann sums again - an integral!

Here's the final result:

The volume of the solid is V=022πxf(x)dx=16π5\displaystyle V=\int_0^2 2\pi x f(x)\, dx=\frac{16\pi}{5}.


In general, the volume of the solid obtained by rotating about the y-axis the region between a curve y=f(x)y=f(x) (f(x)0)(f(x)\ge0) and the x-axis from x=ax=a to x=bx=b is given by V=ab2πxf(x)dx=ab2π(shellradius)(shellheight)dx.V=\int_a^b 2\pi x f(x)\, dx=\int_a^b 2\pi\binom{\mathrm{shell}}{\mathrm{radius}}\binom{\mathrm{shell}}{\mathrm{height}}\, dx.

Example 2

Find the volume of the solid obtained by rotating around the y-axis the region between f(x)=xf(x)=\sqrt{x} and the x-axis from x=0x=0 to x=4x=4.

f(x)=sqrt(x) plot(f,xmin=0,xmax=4,fill='axis')

The volume is V=042πxxdx=128π5\displaystyle V=\int_0^4 2\pi x \sqrt{x}\, dx=\frac{128\pi}{5}.

[Note: we could have done this with washers as well.]

3D rendering not yet implemented

Example 3 - Rotating around a vertical line other than the y-axis

Find the volume of the solid obtained by rotating around the line x=6x=6 the region between f(x)=(x2)(x4)f(x)=-(x-2)(x-4) and the x-axis from x=2x=2 to x=4x=4.

f(x)=-(x-2)*(x-4) plot(f,xmin=2,xmax=4,fill='axis')+line([(6,0),(6,1)],color='black',linestyle='--')

Since we are rotating around the line x=6x=6 and not the y-axis, the shell radius is not given by xx; instead, the shell radius is 6x6-x. The shell height is given by the function, f(x)f(x).

So the volume is V=242π(6x)f(x)dx=8π\displaystyle V=\int_2^4 2\pi(6-x)f(x)\, dx=8\pi.

3D rendering not yet implemented

Example 4 - Rotating around a horizontal line

Find the volume of the solid obtained by rotating around the x-axis the region between the curves y=2x+1y=2\sqrt{x}+1 and y=x+1y=x+1 from x=0x=0 to x=4x=4.

f(x)=2*sqrt(x)+1 g(x)=x+1 plot(f,xmin=0,xmax=4,fill=g(x))+plot(g,xmin=0,xmax=4,ymin=0)

This could be done with washers, but let's try to use cylindrical shells. Here we are rotating around a horizontal line (the x-axis). That means our shell radius will be yy, and our shell height is the horizontal distance between the curves at that given yy.

We need to solve our equations for xx. For the first, we get x=(y1)24x=\frac{(y-1)^2}{4}. For the second, x=y1x=y-1. The shell height is then y1(y1)24y-1-\frac{(y-1)^2}{4}.

Thus, the volume is V=152πy(y1(y1)24)dy=16π\displaystyle V=\int_1^5 2\pi y\left(y-1-\frac{(y-1)^2}{4}\right)\, dy=16\pi.

[Note the limits of integration are y-values. When x=0x=0, then y=1y=1; when x=4x=4, then y=5y=5.]

%var y integral(2*pi*y*(y-1-(y-1)^2/4),y,1,5)
3D rendering not yet implemented

Example 5 - Two Functions with a Gap

Find the volume of the solid obtained by rotating around the x-axis the region between the curves y=x2y=x^2 and y=2x+4y=2x+4 from x=0x=0 to x=2x=2.

Here is a picture of the region. Notice the gap between the two functions.

f(x)=x^2 g(x)=2*x+4 plot(f,xmin=0,xmax=2)+plot(g,xmin=0,xmax=2,fill=f)

When using cylindrical shells to find the volume, the shell radius will be given by yy, but the shell height is not as simple as the last problem. From y=0y=0 to y=4y=4, the height is the distance from the y-axis to the curve y=x2y=x^2 (or x=yx=\sqrt{y}), and from y=4y=4 to y=8y=8, the height is the distance from the curve y=2x+4y=2x+4 (or x=12y2x=\frac{1}{2}y-2) to the vertical line x=2x=2.

This means we need two different integrals, to cover the two cases.

V=042πyydy+482πy(2(12y2))dy=1024π15V=\int_0^4 2\pi y \sqrt{y}\, dy+\int_4^8 2\pi y \left(2-\left(\frac{1}{2}y-2\right)\right)\,dy=\frac{1024\pi}{15}

%var y integral(2*pi*y*sqrt(y),y,0,4)+integral(2*pi*y*(2-(1/2*y-2)),y,4,8)
3D rendering not yet implemented