Learning Objectives
- Recognize when a function of three variables is integrable over a rectangular box.
- Evaluate a triple integral by expressing it as an iterated integral.
- Recognize when a function of three variables is integrable over a closed and bounded region.
- Simplify a calculation by changing the order of integration of a triple integral.
Integrable Functions of Three Variables
We can define a rectangular box B in R3 as B={(x,y,z)∣a ≤ x ≤ b,c ≤ y ≤ d,e ≤ z ≤ f}. We follow a similar procedure to what we did in Double Integrals over Rectangular Regions. We divide the interval [a,b] into l subintervals [xi−1,xi] of equal length Δx=xi−xi−1l, divide the interval [c,d] into m subintervals [yi−1,yi] of equal length Δy=yj−yj−1m, and divide the interval [e,f] into n subintervals [zi−1,zi] of equal length Δz=zk−zk−1n. Then the rectangular box B is subdivided into lmn subboxes Bijk=[xi−1,xi] × [yi−1,yi] × [zi−1,zi], as shown in Figure 1.
Figure 1. A rectangular box in R3 divided into subboxes by planes parallel to the coordinate planes.
For each i, j, and k, consider a sample point (x∗ijk,y∗ijk,z∗ijk) in each sub-box Bijk. We see that its volume is ΔV=ΔxΔyΔz. Form the triple Riemann sum
l∑i=1m∑j=1n∑k=1f(x∗ijk,y∗ijk,z∗ijk)ΔxΔyΔz.
We define the triple integral in terms of the limit of a triple Riemann sum, as we did for the double integral in terms of a double Riemann sum.
definition
The triple integral of a function f(x,y,z) over a rectangular box B is defined as
liml,m,n→∞ l∑i=1m∑j=1n∑k=1 f(x∗ijk,y∗ijk,z∗ijk)ΔxΔyΔz=∭Bf(x,y,z)dV
if this limit exists.
When the triple integral exists on B, the function f(x,y,z) is said to be integrable on B. Also, the triple integral exists if f(x,y,z) is continuous on B. Therefore, we will use continuous functions for our examples. However, continuity is sufficient but not necessary; in other words, f is bounded on B and continuous except possibly on the boundary of B The sample point (x∗ijk,y∗ijk,z∗ijk) can be any point in the rectangular sub-box Bijk and all the properties of a double integral apply to a triple integral. Just as the double integral has many practical applications, the triple integral also has many applications, which we discuss in later sections.
Now that we have developed the concept of the triple integral, we need to know how to compute it. Just as in the case of the double integral, we can have an iterated triple integral, and consequently, a version of Fubini’s thereom for triple integrals exists.
theorem: fubini’s theorem for triple integrals
If f(x,y,z) is continuous on a rectangular box B=[a,b] × [c,d] × [e,f], then
∭Bf(x,y,z)dV=∫fe∫dc∫baf(x,y,z)dx dy dz.
This integral is also equal to any of the other five possible orderings for the iterated triple integral.
For a, b, c, d, e, and f real numbers, the iterated triple integral can be expressed in six different orderings:
∫fe∫dc∫baf(x,y,z)dx dy dz=∫fe(∫dc(∫baf(x,y,z)dx)dy)dz=∫dc(∫fe(∫baf(x,y,z)dx)dz)dy=∫ba(∫fe(∫dcf(x,y,z)dy)dz)dx=∫fe(∫ba(∫dcf(x,y,z)dy)dx)dz=∫ec(∫ba(∫fef(x,y,z)dz)dx)dy=∫ba(∫ec(∫fef(x,y,z)dz)dy)dx.
For a rectangular box, the order of integration does not make any significant difference in the level of difficulty in computation. We compute triple integrals using Fubini’s Theorem rather than using the Riemann sum definition. We follow the order of integration in the same way as we did for double integrals (that is, from inside to outside).
Example: evaluating a triple integral 1
Evaluate the triple integral ∫z=1z=0∫y=4y=2∫x=5x=−1(x+yz2)dx dy dz.
Show Solution
The order of integration is specified in the problem, so integrate with respect to x first, then y, and then z.
∫z=1z=0∫y=4y=2∫x=5x=−1(x+yz2)dx dy dz=∫z=1z=0∫y=4y=2[x22+xyz2∣∣∣x=5x=−1]dy dzIntegrate with respect to x.=∫z=1z=0∫y=4y=2[12+6yz2] dy dzEvaluate.=∫z=1z=0[12y+6y22x2∣∣∣y=4y=2]dzIntegrate with respect to y.=∫z=1z=0[24+36z2]dzEvaluate.=[24z+36z33]z=1z=0=36.
Example: evaluating a triple integral 2
Evaluate the triple integral ∭Bx2yzdV where
B={(x,y,z)∣−2 ≤ x ≤ 1,0 ≤ y ≤ 3,1 ≤ z ≤ 5} as shown in the following figure.
Figure 2. Evaluating a triple integral over a given rectangular box.
Show Solution
The order is not specified, but we can use the iterated integral in any order without changing the level of difficulty. Choose, say, to integrate y first, then x, and then z.
∭Bx2yzdV=∫51∫1−2∫30[x2yz]dy dx dz=∫51∫1−2[x2y22z∣∣∣30]dx dz=∫51∫1−292x2z dx dz=∫51[92x33z∣∣∣1−2]dz=∫51272z22∣∣∣51=162.
Now try to integrate in a different order just to see that we get the same answer. Choose to integrate with respect to x first, then z, and then y.
∭Bx2yzdV=∫30∫51∫1−2[x2yz]dx dz dy=∫30∫51[x33yz∣∣∣1−2]dz dy=∫30∫513yz dz dy=∫30[3yz22∣∣∣51]dy=∫3036y dy=36y22∣∣∣30=18(9−0)=162.
try it
Evaluate the triple integral ∭Bz sin x cos y dV where
B={(x,y,z)∣0 ≤ x ≤ π,3π2 ≤ y ≤ 2π,1 ≤ z ≤ 3}.
Show Solution
∭Bzsinxcosy dV=8
Watch the following video to see the worked solution to the above Try It
Triple Integrals over a General Bounded Region
We now expand the definition of the triple integral to compute a triple integral over a more general bounded region E in R3. The general bounded regions we will consider are of three types. First, let D be the bounded region that is a projection of E onto the xy-plane. Suppose the region E in R3 has the form
E={(x,y,z)∣(x,y) ∈ D,u1(x,y) ≤ z ≤ u2(x,y)}.
For two functions z=u1(x,y) and z=u2(x,y), such that u1(x,y) ≤ u2(x,y) for all (x,y) in D as shown in the following figure.
Figure 3. We can describe region E as the space between u1(x,y) and u2(x,y) above the projection D of E onto the xy-plane.
theorem: triple integral over a general region
The triple integral of a continuous function f(x,y,z) over a general three-dimensional region
E={(x,y,z)∣(x,y) ∈ D,u1(x,y) ≤ z ≤ u2(x,y)}
in R3, where D is the projection of E onto the xy–plane, is
∭Ef(x,y,z)dV=∬D[∫u2(x,y)u1(x,y)f(x,y,z)dz]dA.
Similarly, we can consider a general bounded region D in the xy-plane and two functions y=u1(x,z) and y=u2(x,z) such that u1(x,z) ≤ u2(x,z) for all (x,z) in D. Then we can describe the solid region E in R3 as
E={(x,y,z)∣(x,z) ∈ D,u1(x,z) ≤ z ≤ u2(x,z)}
where D is the projection of E onto the xy-plane and the triple integral is
∭Ef(x,y,z)dV=∬D[∫u2(x,z)u1(x,z)f(x,y,z)dy]dA.
Finally, if D is a general bounded region in the yz-plane and we have two functions x=u1(y,z) and x=u2(y,z) such that u1(y,z) ≤ u2(y,z) for all (y,z) in D, then the solid region E in R3 can be described as
E={(x,y,z)∣(y,z) ∈ D,u1(y,z) ≤ z ≤ u2(y,z)}
where D is the projection of E onto the yz-plane and the triple integral is
∭Ef(x,y,z)dV=∬D[∫u2(y,z)u1(y,z)f(x,y,z)dx]dA.
Note that the region D in any of the planes may be of Type I or Type II as described in Double Integrals over General Regions. If D in the xy-plane is of Type I (Figure 4), then
E={(x,y,z)∣a ≤ x ≤ b,g1(x) ≤ y ≤ g2(x),u1(x,y) ≤ z ≤ u2(x,y)}.
Figure 4. A box E where the projection D in the xy-plane is of Type I.
Then the triple integral becomes
∭Ef(x,y,z)dV=∫ba ∫g2(x)g1(x) ∫u2(x,y)u1(x,y)f(x,y,z)dzdydx.
If D in the xy-plane is of Type II (Figure 5), then
E={(x,y,z)∣c ≤ x ≤ d,h1(x) ≤ y ≤ h2(x),u1(x,y) ≤ z ≤ u2(x,y)}.
Figure 5. A box E where the projection D in the xy-plane is of Type II.
Then the triple integral becomes
∭Ef(x,y,z)dV=∫y=dy=c∫x=h2(y)x=h1(y)∫z=u2(x,y)z=u1(x,y)f(x,y,z)dzdxdy.
Example: evaluating a triple integral over a general bounded region
Evaluate the triple integral of the function f(x,y,z)=5x−3y over the solid tetrahedron bounded by the planes x=0, y=0, z=0, and x+y+z=1.
Show Solution
Figure 6 shows the solid tetrahedron E and its projection D on the xy-plane.
Figure 6. The solid E has a projection D on the xy-plane of Type I.
We can describe the solid region tetrahedron as
E={(x,y,z)∣0 ≤ x ≤ 1,0 ≤ y ≤ 1−x,0 ≤ z ≤ 1−x−y}.
Hence, the triple integral is
∭Ef(x,y,z)dV=∫x=1x=0∫y=1−xy=0∫z=1−x−yz=0(5x−3y)dzdydx.
To simplify the calculation, first evaluate the integral ∫z=1−x−yz=0(5x−3y)dz. We have
∫z=1−x−yz=0(5x−3y)dz=(5x−3y)(1−x−y).
Now evaluate the integral ∫y=1−xy=0(5x−3y)(1−x−y)dy, obtaining
∫y=1−xy=0(5x−3y)(1−x−y)dy=12(x−1)2(6x−1).
Finally, evaluate
∫x=1x=012(x−1)2(6x−1)dx=112.
Putting it all together, we have
∭Ef(x,y,z)dV=∫x=1x=0∫y=1−xy=0∫z=1−x−yz=0(5x−3y)dzdydx=112.
Just as we used the double integral ∬D1dA to find the area of a general bounded region D, we can use ∭E1dV to find the volume of a general solid bounded region E. The next example illustrates the method.
Example: finding a volume by evaluating a triple integral
Find the volume of a right pyramid that has the square base in the xy-plane [−1,1] × [−1,1] and vertex at the point (0,0,1) as shown in the following figure.
Figure 7. Finding the volume of a pyramid with a square base.
Show Solution
In this pyramid the value of z changes from 0 to 1 and at each height z, the cross section of the pyramid for any value of z is the square [−1+z,1−z] × [−1+z,1−z]. Hence, the volume of the pyramid is ∭E1dV where
E={(x,y,z)∣0 ≤ z ≤ 1,−1+z ≤ y ≤ 1−z,−1+z ≤ x ≤ 1−z}.
Thus, we have
∭E1dV=∫z=1z=0∫y=1−zy=−1+z∫x=1−zx=1+z1dxdydz=∫z=1z=0∫y=1−zy=−1+z(2−2z)dydz=∫z=1z=0(2−2z)2dz=43.
Hence, the volume of the pyramid is 43 cubic units.
try it
Consider the solid sphere E=(x,y,z)∣{x2+y2+z2=9}. Write the triple integral ∭Ef(x,y,z)dV for an arbitrary function f as an iterated integral. Then evaluate this triple integral with f(x,y,z)=1. Notice that this gives the volume of a sphere using a triple integral.
Show Solution
∭E1 dV=8∫x=3x=−3∫y=√9−x2y=−√9−x2∫z=√9−x2−y2z=−√9−x2−y21 dz dy dx=36π.
Watch the following video to see the worked solution to the above Try It
Changing the Order of Integration
As we have already seen in double integrals over general bounded regions, changing the order of the integration is done quite often to simplify the computation. With a triple integral over a rectangular box, the order of integration does not change the level of difficulty of the calculation. However, with a triple integral over a general bounded region, choosing an appropriate order of integration can simplify the computation quite a bit. Sometimes making the change to polar coordinates can also be very helpful. We demonstrate two examples here.
Example: changing the order of integration
Consider the iterated integral
∫x=1x=0 ∫y=x2y=0 ∫z=y2z=0f(x,y,z)dzdydx.
The order of integration here is first with respect to z, then y, and then x. Express this integral by changing the order of integration to be first with respect to x, then z, and then y. Verify that the value of the integral is the same if we let f(x,y,z)=xyz.
Show Solution
The best way to do this is to sketch the region E and its projections onto each of the three coordinate planes. Thus, let
E={(x,y,z)∣0 ≤ x ≤ 1,0 ≤ y ≤ x2,0 ≤ z ≤ y2}.
and
∫x=1x=0 ∫y=x2y=0 ∫z=y2z=0f(x,y,z)dzdydx=∭Ef(x,y,z)dV.
We need to express this triple integral as
∫y=dy=c ∫z=v2(y)z=v1(y) ∫x=u2(y,z)x=u1(y,z)f(x,y,z)dxdzdy.
Knowing the region E we can draw the following projections (Figure 8):
on the xy-plane is D1={(x,y)∣0 ≤ x ≤ 1,0 ≤ y ≤ x2}={(x,y)∣0 ≤ y ≤ 1,√y ≤ x ≤ 1},
on the yz-plane is D2={(y,z)∣0 ≤ y ≤ 1,0 ≤ z ≤ y2}, and
on the xz-plane is D3={(x,z)∣0 ≤ x ≤ 1,0 ≤ z ≤ x2}.
Figure 8. The three cross sections of E on the three coordinate planes.
Now we can describe the same region E as {(x,y,z)∣0 ≤ y ≤ 1,0 ≤ z ≤ y2,√y ≤ x ≤ 1}, and consequently, the triple integral becomes
∫y=dy=c ∫z=v2(y)z=v1(y) ∫x=u2(y,z)x=u1(y,z)f(x,y,z)dxdzdy=∫y=1y=0 ∫z=x2z=0 ∫x=1x=√yf(x,y,z)dxdzdy.
Now assume that f(x,y,z)=xyz in each of the integrals. Then we have
[latex]\begin{align}
\displaystyle\int^{x=1}_{x=0}\displaystyle\int^{y=x^2}_{y=0}\displaystyle\int^{z=y^2}_{z=0}xyz{dz}{dy}{dx}&=\displaystyle\int^{x=1}_{x=0}\displaystyle\int^{y=x^2}_{y=0}\left[xy\frac{z^2}2\bigg|_{z=0}^{z=y^2}\right]{dy}{dx} \\ &=\displaystyle\int^{x=1}_{x=0}\displaystyle\int^{y=x^2}_{y=0}\left(x\frac{y^5}2\right)dy \ dx \\ &=\displaystyle\int^{x=1}_{x=0}\left[x\frac{y^6}{12}\bigg|^{y=x^2}_{y=0}\right]dx \\ & =\displaystyle\int^{x=1}_{x=0}\frac{x^{13}}{12}dx =\frac1{168}, \\ \displaystyle\int^{y=1}_{y=0}\displaystyle\int^{z=y^2}_{z=0}\displaystyle\int^{x=1}_{x=\sqrt{y}}xyz{dx}{dz}{dy} &=\displaystyle\int^{y=1}_{y=0}\displaystyle\int^{z=y^2}_{z=0}\left[yz\frac{x^2}2\bigg|^{1}_{\sqrt{y}}\right]dz \ dy\\ &=\displaystyle\int^{y=1}_{y=0}\displaystyle\int^{z=y^2}_{z=0}\left(\frac{yz}2-\frac{y^2z}2\right)dz \ dy \\ &=\displaystyle\int^{y=1}_{y=0}\left[\frac{yz^2}4-\frac{y^2z^2}4\bigg|^{z=y^2}_{z=0}\right]dy \\ &=\displaystyle\int^{y=1}_{y=0}\left(\frac{y^5}4-\frac{y^6}4\right)dy=\frac1{168}. \end{align}[/latex]
The answers match.
try it
Write five different iterated integrals equal to the given integral ∫z=4z=0 ∫y=4−zy=0 ∫x=√yx=0f(x,y,z)dxdydz.
Show Solution
(i) ∫z=4z=0∫x=√x−4x=0∫y=4−zy=x2f(x,y,z)dy dx dz,
(ii)∫y=4y=0∫z=4−yz=0∫x=√yx=0f(x,y,z)dx dz dy,
(iii)∫y=4y=0∫x=√yx=0∫z=4−yz=0f(x,y,z)dz dx dy ,
(iv) ∫x=2x=0∫y=4y=x2∫z=4−yz=0f(x,y,z)dz dy dx,
(v) ∫x=2x=0∫z=4−x2z=0∫y=4−yy=x2f(x,y,z)dy dz dx.
Example: changing integration order and coordinate systems
Evaluate the triple integral ∭E√x2+z2dV, where E is the region bounded by the paraboloid y=x2+z2 (Figure 9) and the plane y=4.
Figure 9. Integrating a triple integral over a paraboloid.
Show Solution
The projection of the solid region E onto the xy-plane is the region bounded above by y=4 and below by the parabola y=x2 as shown.
Figure 10. Cross section in the xy-plane of the paraboloid in Figure 9.
Thus, we have
E={(x,y,z)∣−2 ≤ x ≤ 2,x2 ≤ y ≤ 4,−√y−x2 ≤ z ≤ √y−x2}.
The triple integral becomes
∭E√x2+z2dV=∫x=2x=−2 ∫y=4y=x2 ∫z=√y−x2z=−√y−x2√x2+z2dzdydx.
This expression is difficult to compute, so consider the projection of E onto the xz-plane. This is a circular disc x2+z2 ≤ 4. So we obtain
∭E√x2+z2dV=∫x=2x=−2 ∫y=4y=x2 ∫z=√y−x2z=−√y−x2√x2+z2dzdydx=∫x=2x=−2 ∫z=√4−x2z=−√4−x2 ∫y=4y=x2+z2√x2+z2dydzdx.
Here the order of integration changes from being first with respect to z, then y and then x to being first with respect to y, then to z, and then to x. It will soon be clear how this change can be beneficial for computation. We have
∫x=2x=−2 ∫z=√4−x2z=−√4−x2 ∫y=4y=x2+z2√x2+z2dydzdx=∫x=2x=−2 ∫z=√4−x2z=−√4−x2(4−x2−z2)√x2+z2dzdx.
Now use the polar substitution x=rcosθ,z=rsinθ, and dzdx=rdrdθ in the xz-plane. This is essentially the same thing as when we used polar coordinates in the xy-plane, except we are replacing y by z. Consequently the limits of integration change and we have, by using r2=x2+z2,
∫x=2x=−2 ∫z=√4−x2z=−√4−x2(4−x2−z2)√x2+z2dzdx=∫θ=2πθ=0 ∫r=2r=0(4−r2)rrdrdθ
=∫2π0 [4r33−r55∣∣∣20] dθ=∫2π0 6415dθ=128π15.
Candela Citations
CC licensed content, Original
CC licensed content, Shared previously