8.4: Double and Triple Integrals
- Page ID
- 106851
\( \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}\)We can extend the idea of a definite integral to more dimensions. If \(f(x,y)\) is continuous over the rectangle \(R=[a,b]\times[c,d]\) then,
\[\label{c2v:eq:doubleint} \underset{R}{\iint}f(x,y)dA=\int_{a}^{b}\int_{c}^{d}f{(x,y)\, dy\, dx}=\int_{c}^{d}\int_{a}^{b}f{(x,y)\, dx\, dy} \]
If \(f(x,y)\geq0\), then the double integral represents the volume \(V\) of the solid that lies above the rectangle \(R\) and below the surface \(z=f(x,y)\) (Figure \(\PageIndex{1}\)).
![doubleintegral.jpg](https://chem.libretexts.org/@api/deki/files/141547/doubleintegral.jpg?revision=1&size=bestfit&width=478&height=417)
We can compute the double integral of Equation \ref{c2v:eq:doubleint} as:
\[\underset{R}{\iint}f(x,y)dA=\int_{a}^{b}\left[ \int_{c}^{d}f{(x,y)\, dy} \right]dx \nonumber \]
meaning that we will first compute
\[\int_{c}^{d}f{(x,y)\, dy} \nonumber \]
holding \(x\) constant and integrating with respect to \(y\). The result will be a function containing only \(x\), which we will integrate between \(a\) and \(b\) with respect to \(x\).
For example, let’s solve \(\int_{0}^{3}\int_{1}^{2}{x^2 y\; dy\; dx}\). We’ll start by solving \(\int_{1}^{2}{x^2 y\; dy}\) holding \(x\) constant:
\[\int_{1}^{2}{x^2 y\; dy}=x^2\int_{1}^{2}{y\; dy}=\dfrac{3}{2}x^2 \nonumber \]
Now we integrate this function from 0 to 3 with respect to \(x\):
\[\int_{0}^{3}\int_{1}^{2}{x^2 y\; dy\; dx}=\int_{0}^{3}{\dfrac{3}{2}x^2 \; dx}=\dfrac{27}{2} \nonumber \]
You can of course integrate from 0 to 3 first with respect to \(x\) holding \(y\) constant, and then integrate the result with respect to \(y\) from 1 to 2. Try it this way and verify you get the same result.
Triple integrals work in the same way. If \(f(x,y,z)\) is continuous on the rectangular box \(B=[a,b]\times[c,d]\times[r,s]\), then
\[\label{c2v:eq:tripleint} \underset{B}{\iiint}f(x,y,z)dV=\int_{r}^{s}\int_{c}^{d}\int_{a}^{b}f{(x,y)\, dx\, dy\, dz} \]
This iterated integral means that we integrate first with respect to \(x\) (keeping \(y\) and \(z\) fixed), then we integrate with respect to \(y\) (keeping \(z\) fixed), and finally we integrate with respect to \(z\). There are five other possible orders in which we can integrate, all of which give the same value.
Do you need a refresher on double and triple integrals? Check the videos below before moving on to the physical chemistry examples.
- Example 1: http://www.youtube.com/watch?v=RqD89-afGS0
- Example 2: http://www.youtube.com/watch?v=CPR0ZD0IYVE (check the example that starts around 3:45 min. and ends at 5:07 min)
Triple integrals are used very often in physical chemistry to normalize probability density functions. For example, in quantum mechanics, the absolute square of the wave function, \(\left | \psi (x,y,z) \right |^2\), is interpreted as a probability density, the probability that the particle is inside the volume \(dx.dy.dz\). Since the probability of finding the particle somewhere in space is 1, we require that:
\[\label{c2v:eq:calculus2v_normalization} \int_{-\infty }^{\infty }\int_{-\infty }^{\infty }\int_{-\infty }^{\infty }{\left | \psi (x,y,z) \right |}^2\; dx \;dy \;dz=1 \]
We already mentioned wave functions in Section 2.3, where we showed that
\[\left | \psi (x,y,z) \right |^2=\psi^*(x,y,z) \psi(x,y,z) \nonumber \]
The normalization condition, therefore, can also be written as
\[\label{c2v:eq:calculus2v_normalization2} \int_{-\infty }^{\infty }\int_{-\infty }^{\infty }\int_{-\infty }^{\infty }{ \psi^*\psi }\; dx \;dy \;dz=1 \]
In quantum mechanics, the lowest energy state of a particle confined in a three-dimensional box is represented by
\[\psi (x,y,z)=A\sin{\dfrac{\pi x}{a}}\sin{\dfrac{\pi y}{b}}\sin{\dfrac{\pi z}{c}}\; if\; \left\{ \begin{matrix} 0< x< a\\ 0< y< b\\ 0< z< c\\ \end{matrix} \right. \nonumber \]
and
\(\psi (x,y,z)=0\) otherwise (outside the box).
Here, \(A\) is a normalization constant, and \(a\),\(b\) and \(c\) are the lengths of the sides of the box. Since the probability of finding the particle somewhere in space is 1, we require that
\[\int_{-\infty }^{\infty }\int_{-\infty }^{\infty }\int_{-\infty }^{\infty }{\left | \psi (x,y,z) \right |}^2\; dx \;dy \;dz=1 \nonumber \]
Find the normalization constant \(A\) in terms of \(a,b,c\) and other constants.
Solution
Because \(\psi(x,y,z)=0\) outside the box,
\[\int_{-\infty }^{\infty }\int_{-\infty }^{\infty }\int_{-\infty }^{\infty }{\left | \psi (x,y,z) \right |}^2\; dx \;dy \;dz=\int_{0 }^{c }\int_{0 }^{b }\int_{0 }^{a }{\left | \psi (x,y,z) \right |}^2\; dx \;dy \;dz = \left | \psi (x,y,z) \right |^2=\psi^* (x,y,z) \psi(x,y,z) \nonumber \]
However, because in this case the function is real,
\[\left | \psi (x,y,z) \right |^2=\left ( \psi (x,y,z) \right )^2 \nonumber \]
\[\int_{0 }^{c }\int_{0 }^{b }\int_{0 }^{a }{\left | \psi (x,y,z) \right |}^2\; dx \;dy \;dz=\int_{0 }^{c }\int_{0 }^{b }\int_{0 }^{a }{A^2 \sin^2\left( \dfrac{\pi x}{a}\right)\sin^2\left( \dfrac{\pi y}{b}\right)\sin^2\left( \dfrac{\pi z}{c}\right)}\; dx \;dy \;dz=1 \nonumber \]
\[\int_{0 }^{c }\int_{0 }^{b }\int_{0 }^{a }{A^2 \sin^2\left( \dfrac{\pi x}{a}\right)\sin^2\left( \dfrac{\pi y}{b}\right)\sin^2\left( \dfrac{\pi z}{c}\right)}\; dx \;dy \;dz= A^2\int_{0 }^{a }{\sin^2\left( \dfrac{\pi x}{a}\right)}dx\int_{0 }^{b }{\sin^2\left( \dfrac{\pi y}{b}\right)}dy\int_{0 }^{c }{\sin^2\left( \dfrac{\pi z}{c}\right)}dz \nonumber \]
Using the formula sheet, we get
\[\int_{0 }^{a }{\sin^2\left( \dfrac{\pi x}{a}\right)}dx=a/2 \nonumber \]
And therefore,
\[A^2\int_{0 }^{a }{\sin^2\left( \dfrac{\pi x}{a}\right)}dx\int_{0 }^{b }{\sin^2\left( \dfrac{\pi y}{b}\right)}dy\int_{0 }^{c }{\sin^2\left( \dfrac{\pi z}{c}\right)}dz=A^2 \dfrac{a}{2} \dfrac{b}{2}\dfrac{c}{2}=\dfrac{A^2abc}{8}=1 \nonumber \]
Solving for \(A\):
\[\displaystyle{\color{Maroon}A=\left( \dfrac{8}{abc}\right)^{1/2}} \nonumber \]