15.5: Triple Integrals
- Page ID
- 149585
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)It will come as no surprise that we can also do triple integrals---integrals over a three-dimensional region. The simplest application allows us to compute volumes in an alternative way.
To approximate a volume in three dimensions, we can divide the three-dimensional region into small rectangular boxes, each \(\Delta x\times\Delta y\times\Delta z\) with volume \(\Delta x\Delta y\Delta z\). Then we add them all up and take the limit, to get an integral: \[\int_{x_0}^{x_1} \int_{y_0}^{y_1} \int_{z_0}^{z_1} \, dz\,dy\,dx. \nonumber\] If the limits are constant, we are simply computing the volume of a rectangular box.
We use an integral to compute the volume of the box with opposite corners at \((0,0,0)\) and \((1,2,3)\).
Solution
\[\int_0^1 \int_0^2 \int_0^3 \;dz\,dy\,dx= \int_0^1 \int_0^2 \left.z\right|_0^3 \,dy\,dx =\int_0^1 \int_0^2 3\,dy\,dx =\int_0^1 \left.3y\right|_0^2 \,dx =\int_0^1 6\,dx = 6. \nonumber\]
Of course, this is more interesting and useful when the limits are not constant.
Find the volume of the tetrahedron with corners at \((0,0,0)\), \((0,3,0)\), \((2,3,0)\), and \((2,3,5)\).
Solution
The whole problem comes down to correctly describing the region by inequalities: \(0\le x\le 2,\;\) \(3x/2\le y\le 3,\;\) \(0\le z\le 5x/2\). The lower \(y\) limit comes from the equation of the line \(y=3x/2\) that forms one edge of the tetrahedron in the \(x\)-\(y\) plane; the upper \(z\) limit comes from the equation of the plane \(z=5x/2\) that forms the "upper" side of the tetrahedron; see figure \(\PageIndex{1}\). Now the volume is
\[\eqalign{ \int_0^2 \int_{3x/2}^3 \int_0^{5x/2}\;dz\,dy\,dx&= \int_0^2 \int_{3x/2}^3 \left.z\right|_{z=0}^{5x/2} \,dy\,dx\cr&= \int_0^2 \int_{3x/2}^3 \dfrac{5x}{2}\,dy\,dx\cr&= \int_0^2 \left. \dfrac{5x}{2}y \right|_{y=3x/2}^3 \,dx\cr&= \int_0^2 \left(\dfrac{15x}{2}-\dfrac{15x^2}{4} \right)\,dx\cr&= \left| \dfrac{15x^2}{4}- \dfrac{15x^3}{12} \right|_0^2\cr&= 15-10=5.\cr \nonumber}\]
Figure \(\PageIndex{1}\): A tetrahedron.
Pretty much just the same way as we did for two dimensions, we can use triple integration to compute mass, center of mass, and various average quantities.
Suppose the temperature at a point is given by \(T=xyz\). Find the average temperature in the cube which has opposite corners at \((0,0,0)\) and \((2,2,2)\).
Solution
In two dimensions we add up the temperature at "each" point and divide by the area; here we add up the temperatures and divide by the volume of the cube, which is \(8\): \[\eqalign{ \dfrac{1}{8} \int_{0}^2 \int_{0}^2 \int_{0}^2 xyz\,dz\,dy\,dx&= \dfrac{1}{8} \int_{0}^2 \int_{0}^2 \left.\dfrac{xyz^2}{2}\right|_0^2\,dy\,dx= \dfrac{1}{16} \int_{0}^2 \int_{0}^2 xy\,dy\,dx\cr&= \dfrac{1}{4} \int_{0}^2 \left. \dfrac{xy^2}{2} \right|_0^2\,dx= \dfrac{1}{8} \int_{0}^2 4x\,dx= \dfrac{1}{2} \left.\dfrac{x^2}{2} \right|_0^2= 1.\cr \nonumber}\]
Suppose the density of an object is given by \(xz\), and the object occupies the tetrahedron with corners \((0,0,0)\), \((0,1,0)\), \((1,1,0)\), and \((0,1,1)\). Find the mass and center of mass of the object.
Solution
As usual, the mass is the integral of density over the region:
\[\eqalign{ M&= \int_0^1 \int_{x}^1 \int_0^{y-x} xz\;dz\,dy\,dx= \int_0^1 \int_x^1 \dfrac{x(y-x)^2}{2}\,dy\,dx= \dfrac{1}{2} \int_0^1 \dfrac{x(1-x)^3}{3}\,dx\cr&= \dfrac{1}{6} \int_0^1 (x-3x^2+3x^3-x^4)\,dx= \dfrac{1}{120}.\cr \nonumber}\]
We compute moments as before, except now there is a third moment: \[\eqalign{ M_{xy} &= \int_0^1 \int_{x}^1 \int_0^{y-x} xz^2\,dz\,dy\,dx= \dfrac{1}{360},\cr M_{xz} &= \int_0^1 \int_x^1 \int_0^{y-x} xyz\,dz\,dy\,dx= \dfrac{1}{144},\cr M_{yz} &= \int_0^1 \int_x^1 \int_0^{y-x} x^2z\,dz\,dy\,dx= \dfrac{1}{360}.\cr} \nonumber\]
Finally, the coordinates of the center of mass are \(\bar x=M_{yz}/M=1/3,\;\) \(\bar y=M_{xz}/M=5/6,\;\) and \(\;\bar z=M_{xy}/M=1/3\).
Exercises \(\PageIndex{}\)
Evaluate \(\displaystyle \int_0^1 \int_0^x \int_0^{x+y} (2x+y-1)\,dz\,dy\,dx.\)
- Answer
-
\(\dfrac{11}{24}\)
Evaluate \(\displaystyle \int_0^2 \int_{-1}^{x^2} \int_1^y xyz \,dz\,dy\,dx.\)
- Answer
-
\(\dfrac{623}{60}\)
Evaluate \(\displaystyle \int_0^1 \int_0^x \int_0^{\ln y} e^{x+y+z} \,dz\,dy\,dx.\)
- Answer
-
\(2e-\dfrac{3(e^2+1)}{4}\)
Evaluate \(\displaystyle \int_0^{\pi/2} \int_0^{\sin\theta} \int_0^{r\cos\theta} r^2 \,dz\,dr\,d\theta.\)
- Answer
-
\(\dfrac{1}{20}\)
Evaluate \(\displaystyle \int_0^{\pi} \int_0^{\sin\theta} \int_0^{r\sin\theta} r\cos^2\theta \,dz\,dr\,d\theta.\)
- Answer
-
\(\dfrac{\pi}{48}\)
Evaluate \(\displaystyle \int_0^1 \int_0^{y^2} \int_0^{x+y} x \,dz\,dx\,dy.\)
- Answer
-
\(\dfrac{11}{84}\)
Evaluate \(\displaystyle \int_1^2 \int_y^{y^2} \int_0^{\ln(y+z)} e^x \,dx\,dz\,dy.\)
- Answer
-
\(\dfrac{151}{60}\)
Evaluate \(\displaystyle \int_0^{\pi} \int_0^{\pi/2} \int_0^1 z(\sin x + \cos y) \,dz\,dy\,dx.\)
- Answer
-
\(\pi\)
For each of the integrals in the previous exercises, give a description of the volume (both algebraic and geometric) that is the domain of integration.
Compute \(\displaystyle \iiint (x+y+z)\,dV\) over the region \(x^2+y^2+z^2 \le 1\) in the first octant.
- Answer
-
\(\dfrac{3\pi}{16}\)
Find the mass of a cube with edge length \(2\) and density equal to the square of the distance from one corner.
- Answer
-
\(32\)
Find the mass of a cube with edge length \(2\) and density equal to the square of the distance from one edge.
- Answer
-
\(\dfrac{64}{3}\)
An object occupies the volume of the upper hemisphere of \(x^2+y^2+z^2=4\) and has density \(z\) at \((x,y,z)\). Find the center of mass.
- Answer
-
\((\bar x,\bar y,\bar z)=(0,0,16/15)\)
An object occupies the volume of the pyramid with corners at \((1,1,0), (1,−1,0), (−1,−1,0), (−1,1,0),\) and \((0,0,2),\) and has density \(x^2+y^2\) at \((x,y,z)\). Find the center of mass.
- Answer
-
\((\bar x,\bar y,\bar z)=(0,0,1/3)\)
Verify the moments \(M_{xy}, M_{xz},\) and \(M_{yz}\) of example \(\PageIndex{4}\) by evaluating the integrals.
Find the region \(E\) for which \(\displaystyle \iiint\limits_{E} (1−x^2−y^2−z^2)\,dV\) is a maximum.
Contributors
Integrated by Justin Marshall.