Section 15.5 : Triple Integrals
Now that we know how to integrate over a two-dimensional region we need to move on to integrating over a three-dimensional region. We used a double integral to integrate over a two-dimensional region and so it shouldn’t be too surprising that we’ll use a triple integral to integrate over a three dimensional region. The notation for the general triple integrals is,
∭Ef(x,y,z)dV∭Ef(x,y,z)dVLet’s start simple by integrating over the box,
B=[a,b]×[c,d]×[r,s]Note that when using this notation we list the x’s first, the y’s second and the z’s third.
The triple integral in this case is,
∭Bf(x,y,z)dV=∫sr∫dc∫baf(x,y,z)dxdydzNote that we integrated with respect to x first, then y, and finally z here, but in fact there is no reason to the integrals in this order. There are 6 different possible orders to do the integral in and which order you do the integral in will depend upon the function and the order that you feel will be the easiest. We will get the same answer regardless of the order however.
Let’s do a quick example of this type of triple integral.
Just to make the point that order doesn’t matter let’s use a different order from that listed above. We’ll do the integral in the following order.
∭B8xyzdV=∫21∫32∫108xyzdzdxdy=∫21∫324xyz2|10dxdy=∫21∫324xydxdy=∫212x2y|32dy=∫2110ydy=15Before moving on to more general regions let’s get a nice geometric interpretation about the triple integral out of the way so we can use it in some of the examples to follow.
Fact
The volume of the three-dimensional region E is given by the integral,
V=∭EdVLet’s now move on the more general three-dimensional regions. We have three different possibilities for a general region. Here is a sketch of the first possibility.

In this case we define the region E as follows,
E={(x,y,z)|(x,y)∈D,u1(x,y)≤z≤u2(x,y)}where (x,y)∈D is the notation that means that the point (x,y) lies in the region D from the xy-plane. In this case we will evaluate the triple integral as follows,
∭Ef(x,y,z)dV=∬D[∫u2(x,y)u1(x,y)f(x,y,z)dz]dAwhere the double integral can be evaluated in any of the methods that we saw in the previous couple of sections. In other words, we can integrate first with respect to x, we can integrate first with respect to y, or we can use polar coordinates as needed.
We should first define octant. Just as the two-dimensional coordinates system can be divided into four quadrants the three-dimensional coordinate system can be divided into eight octants. The first octant is the octant in which all three of the coordinates are positive.
Here is a sketch of the plane in the first octant.


We now need to determine the region D in the xy-plane. We can get a visualization of the region by pretending to look straight down on the object from above. What we see will be the region D in the xy-plane. So D will be the triangle with vertices at (0,0), (3,0), and (0,2). Here is a sketch of D.

Now we need the limits of integration. Since we are under the plane and in the first octant (so we’re above the plane z=0) we have the following limits for z.
0≤z≤6−2x−3yWe can integrate the double integral over D using either of the following two sets of inequalities.
0≤x≤30≤y≤−23x+20≤x≤−32y+30≤y≤2Since neither really holds an advantage over the other we’ll use the first one. The integral is then,
∭E2xdV=∬D[∫6−2x−3y02xdz]dA=∬D2xz|6−2x−3y0dA=∫30∫−23x+202x(6−2x−3y)dydx=∫30(12xy−4x2y−3xy2)|−23x+20dx=∫3043x3−8x2+12xdx=(13x4−83x3+6x2)|30=9Let’s now move onto the second possible three-dimensional region we may run into for triple integrals. Here is a sketch of this region.

For this possibility we define the region E as follows,
E={(x,y,z)|(y,z)∈D,u1(y,z)≤x≤u2(y,z)}So, the region D will be a region in the yz-plane. Here is how we will evaluate these integrals.
∭Ef(x,y,z)dV=∬D[∫u2(y,z)u1(y,z)f(x,y,z)dx]dAAs with the first possibility we will have two options for doing the double integral in the yz-plane as well as the option of using polar coordinates if needed.
In this case we’ve been given D and so we won’t have to really work to find that. Here is a sketch of the region D as well as a quick sketch of the plane and the curves defining D projected out past the plane so we can get an idea of what the region we’re dealing with looks like.


Now, the graph of the region above is all okay, but it doesn’t really show us what the region is. So, here is a sketch of the region itself.


Here are the limits for each of the variables.
0≤y≤434y≤z≤32√y0≤x≤8−y−zThe volume is then,
V=∭EdV=∬D[∫8−y−z0dx]dA=∫40∫3√y/23y/48−y−zdzdy=∫40(8z−yz−12z2)|3√y23y4dy=∫4012y12−578y−32y32+3332y2dy=(8y32−5716y2−35y52+1132y3)|40=495We now need to look at the third (and final) possible three-dimensional region we may run into for triple integrals. Here is a sketch of this region.

In this final case E is defined as,
E={(x,y,z)|(x,z)∈D,u1(x,z)≤y≤u2(x,z)}and here the region D will be a region in the xz-plane. Here is how we will evaluate these integrals.
∭Ef(x,y,z)dV=∬D[∫u2(x,z)u1(x,z)f(x,y,z)dy]dAwhere we will can use either of the two possible orders for integrating D in the xz-plane or we can use polar coordinates if needed.
Here is a sketch of the solid E.


The region D in the xz-plane can be found by “standing” in front of this solid and we can see that D will be a disk in the xz-plane. This disk will come from the front of the solid and we can determine the equation of the disk by setting the elliptic paraboloid and the plane equal.
2x2+2z2=8⇒x2+z2=4This region, as well as the integrand, both seems to suggest that we should use something like polar coordinates. However, we are in the xz-plane and we’ve only seen polar coordinates in the xy-plane. This is not a problem. We can always “translate” them over to the xz-plane with the following definition.
x=rcosθz=rsinθSince the region doesn’t have y’s we will let z take the place of y in all the formulas. Note that these definitions also lead to the formula,
x2+z2=r2With this in hand we can arrive at the limits of the variables that we’ll need for this integral.
2x2+2z2≤y≤80≤r≤20≤θ≤2πThe integral is then,
∭E√3x2+3z2dV=∬D[∫82x2+2z2√3x2+3z2dy]dA=∬D(y√3x2+3z2)|82x2+2z2dA=∬D√3(x2+z2)(8−(2x2+2z2))dANow, since we are going to do the double integral in polar coordinates let’s get everything converted over to polar coordinates. The integrand is,
√3(x2+z2)(8−(2x2+2z2))=√3r2(8−2r2)=√3r(8−2r2)=√3(8r−2r3)The integral is then,
∭E√3x2+3z2dV=∬D√3(8r−2r3)dA=√3∫2π0∫20(8r−2r3)rdrdθ=√3∫2π0(83r3−25r5)|20dθ=√3∫2π012815dθ=256√3π15