Skip to main content
\(\newcommand{\R}{\mathbb{R}} \newcommand{\va}{\mathbf{a}} \newcommand{\vb}{\mathbf{b}} \newcommand{\vc}{\mathbf{c}} \newcommand{\vC}{\mathbf{C}} \newcommand{\vd}{\mathbf{d}} \newcommand{\ve}{\mathbf{e}} \newcommand{\vi}{\mathbf{i}} \newcommand{\vj}{\mathbf{j}} \newcommand{\vk}{\mathbf{k}} \newcommand{\vn}{\mathbf{n}} \newcommand{\vm}{\mathbf{m}} \newcommand{\vr}{\mathbf{r}} \newcommand{\vs}{\mathbf{s}} \newcommand{\vu}{\mathbf{u}} \newcommand{\vv}{\mathbf{v}} \newcommand{\vw}{\mathbf{w}} \newcommand{\vx}{\mathbf{x}} \newcommand{\vy}{\mathbf{y}} \newcommand{\vz}{\mathbf{z}} \newcommand{\vzero}{\mathbf{0}} \newcommand{\vF}{\mathbf{F}} \newcommand{\vG}{\mathbf{G}} \newcommand{\vH}{\mathbf{H}} \newcommand{\vR}{\mathbf{R}} \newcommand{\vT}{\mathbf{T}} \newcommand{\vN}{\mathbf{N}} \newcommand{\vL}{\mathbf{L}} \newcommand{\vB}{\mathbf{B}} \newcommand{\vS}{\mathbf{S}} \newcommand{\proj}{\text{proj}} \newcommand{\comp}{\text{comp}} \newcommand{\nin}{} \newcommand{\vecmag}[1]{|#1|} \newcommand{\grad}{\nabla} \DeclareMathOperator{\curl}{curl} \DeclareMathOperator{\divg}{div} \newcommand{\lt}{<} \newcommand{\gt}{>} \newcommand{\amp}{&} \)

Section3.2Iterated Integrals

Objectives

  • How do we evaluate a double integral over a rectangle as an iterated integral, and why does this process work?

Recall that in Section3.1 we defined the double integral of a continuous function \(f = f(x,y)\) over a rectangle \(R = [a,b] \times [c,d]\) as a limit of double Riemann sums, but while this definition tells us exactly what a double integral is, it is not very helpful for determining the value of a double integral. Fortunately, there is a way to view a double integral as an iterated integral, which will make computations feasible in many cases.

The viewpoint of an iterated integral is closely connected to an important idea from single-variable calculus. When we studied solids of revolution, such as the one shown in Figure3.2.1, we saw that in some circumstances we could slice the solid perpendicular to an axis and have each slice be approximately a circular disk. From there, we were able to find the volume of each disk, and then use an integral to add the volumes of the slices. In what follows, we are able to use single integrals to generalize this approach to handle even more general geometric shapes.

Figure3.2.1A solid of revolution.
Exploration3.2.1

Let \(f(x,y) = 25-x^2-y^2\) on the rectangular domain \(R = [-3,3] \times [-4,4]\text{.}\)

As with partial derivatives, we may treat one of the variables in \(f\) as constant and think of the resulting function as a function of a single variable. Now we investigate what happens if we integrate instead of differentiate.

  1. Choose a fixed value of \(x\) in the interior of \([-3,3]\text{.}\) Let

    \begin{equation*} A(x) = \int_{-4}^4 f(x,y) \, dy. \end{equation*}

    What is the geometric meaning of the value of \(A(x)\) relative to the surface defined by \(f\text{.}\) (Hint: Think about the trace determined by the fixed value of \(x\text{,}\) and consider how \(A(x)\) is related to the image at left in Figure3.2.2.)

    Figure3.2.2Left: A cross section with fixed \(x\text{.}\) Right: A cross section with fixed \(x\) and \(\Delta x\text{.}\)

  2. For a fixed value of \(x\text{,}\) say \(x_i^*\text{,}\) what is the geometric meaning of \(A(x_i^*) \ \Delta x\text{?}\) (Hint: Consider how \(A(x_i^*) \Delta x\) is related to the image at right in Figure3.2.2.)

  3. Since \(f\) is continuous on \(R\text{,}\) we can define the function \(A = A(x)\) at every value of \(x\) in \([-3,3]\text{.}\) Now think about subdividing the \(x\)-interval \([-3,3]\) into \(m\) subintervals, and choosing a value \(x_i^*\) in each of those subintervals. What will be the meaning of the sum \(\sum_{i=1}^m A(x_i^*) \ \Delta x\text{?}\)

  4. Explain why \(\int_{-3}^3 A(x) \, dx\) will determine the exact value of the volume under the surface \(z = f(x,y)\) over the rectangle \(R\text{.}\)

Subsection3.2.1Iterated Integrals

The ideas that we explored in Preview Activity3.2.1 work more generally and lead to the idea of an iterated integral. Let \(f\) be a continuous function on a rectangular domain \(R = [a,b] \times [c,d]\text{,}\) and let

\begin{equation*} A(x) = \int_c^d f(x,y) \, dy. \end{equation*}

The function \(A = A(x)\) determines the value of the cross sectional area (by area we mean signed area) in the \(y\) direction for the fixed value of \(x\) of the solid bounded between the surface defined by \(f\) and the \(xy\)-plane.

Figure3.2.3Summing volumes of cross section slices.

The value of this cross sectional area is determined by the input \(x\) in \(A\text{.}\) Since \(A\) is a function of \(x\text{,}\) it follows that we can integrate \(A\) with respect to \(x\text{.}\) In doing so, we use a partition of \([a, b]\) and make an approximation to the integral given by

\begin{equation*} \int_a^b A(x) \, dx \approx \sum_{i=1}^m A(x_i^*) \Delta x, \end{equation*}

where \(x_i^*\) is any number in the subinterval \([x_{i-1},x_i]\text{.}\) Each term \(A(x_i^*) \Delta x\) in the sum represents an approximation of a fixed cross sectional slice of the surface in the \(y\) direction with a fixed width of \(\Delta x\) as illustrated in Figure3.2.3. We add the signed volumes of these slices as shown in the frames in Figure3.2.3 to obtain an approximation of the total signed volume.

As we let the number of subintervals in the \(x\) direction approach infinity, we can see that the Riemann sum \(\sum_{i=1}^m A(x_i^*) \Delta x\) approaches a limit and that limit is the sum of signed volumes bounded by the function \(f\) on \(R\text{.}\) Therefore, since \(A(x)\) is itself determined by an integral, we have

\begin{equation*} \iint_R f(x,y) \, dA = \lim_{m \to \infty} \sum_{i=1}^m A(x_i^*) \Delta x = \int_a^b A(x) \, dx = \int_a^b \left( \int_c^d f(x,y) \, dy \right) \, dx. \end{equation*}

Hence, we can compute the double integral of \(f\) over \(R\) by first integrating \(f\) with respect to \(y\) on \([c, d]\text{,}\) then integrating the resulting function of \(x\) with respect to \(x\) on \([a, b]\text{.}\) The nested integral

\begin{equation*} \int_a^b \left( \int_c^d f(x,y) \, dy \right) \, dx = \int_a^b \int_c^d f(x,y) \, dy \, dx \end{equation*}

is called an iterated integral, and we see that each double integral may be represented by two single integrals.

We made a choice to integrate first with respect to \(y\text{.}\) The same argument shows that we can also find the double integral as an iterated integral integrating with respect to \(x\) first, or

\begin{equation*} \iint_R f(x,y) \, dA = \int_c^d \left( \int_a^b f(x,y) \, dx \right) \, dy = \int_c^d \int_a^b f(x,y) \, dx \, dy. \end{equation*}

The fact that integrating in either order results in the same value is known as Fubini's Theorem.

Fubini's Theorem

If \(f = f(x,y)\) is a continuous function on a rectangle \(R = [a,b] \times [c,d]\text{,}\) then

\begin{equation*} \iint_R f(x,y) \, dA = \int_c^d \int_a^b f(x,y) \, dx \, dy = \int_a^b \int_c^d f(x,y) \, dy \, dx. \end{equation*}

Fubini's theorem enables us to evaluate iterated integrals without resorting to the limit definition. Instead, working with one integral at a time, we can use the Fundamental Theorem of Calculus from single-variable calculus to find the exact value of each integral, starting with the inner integral.

Activity3.2.2

Let \(f(x,y) = 25-x^2-y^2\) on the rectangular domain \(R = [-3,3] \times [-4,4]\text{.}\)

  1. Viewing \(x\) as a fixed constant, use the Fundamental Theorem of Calculus to evaluate the integral

    \begin{equation*} A(x) = \int_{-4}^4 f(x,y) \, dy. \end{equation*}

    Note that you will be integrating with respect to \(y\text{,}\) and holding \(x\) constant. Your result should be a function of \(x\) only.

  2. Next, use your result from (a) along with the Fundamental Theorem of Calculus to determine the value of \(\int_{-3}^3 A(x) \, dx\text{.}\)

  3. What is the value of \(\iint_R f(x,y) \, dA\text{?}\) What are two different ways we may interpret the meaning of this value?

Example3.2.4

Integrate the function \(xe^{xy}\) over the rectangle \([0,1]\times [0,1]\) by first writing it as an iterated integral. Does one order of integration offer any advantages over the other?

Solution

Notice that if we set this up as a \(dxdy\) integral, we'd have to integrate by parts:

\begin{equation*} \int_{y=0}^1 \int_{x=0}^1 xe^{xy} \, dxdy \end{equation*}

So let's try the other order instead.

\begin{align*} \iint_{[0,1]\times [0,1]} xe^{xy}\,dA \amp = \int_{x=0}^1 \left[ \int_{y=0}^1 xe^{xy}\,dy\right]\,dx \\ \amp = \int_{x=0}^1 \left[ \frac{x}{x} e^{xy} \bigg|_{y=0}^1 \right] \,dx \\ \amp = \int_{x=0}^1 e^{x} - 1 \,dx \\ \amp = e-2 \end{align*}
Example3.2.5

Let \(f(x,y) = x+y^2\) on the rectangle \(R = [0,2] \times [1,3]\text{.}\)

  1. Evaluate \(\iint_R f(x,y) \, dA\) using an iterated integral. Choose an order for integration by deciding whether you want to integrate first with respect to \(x\) or \(y\text{.}\)

  2. Evaluate \(\iint_R f(x,y) \, dA\) using the iterated integral whose order of integration is the opposite of the order you chose in (a).

Solution

  1. We'll choose to integrate with respect to \(y\) first:

    \begin{align*} \iint_{[0,2]\times [1,3]} x + y^2\,dA \amp = \int_{x=0}^2 \int_{y=1}^3 x + y^2\,dy dx\\ \amp = \int_{x=0}^2 \left[ xy + \frac{y^3}{3} \bigg|_{y=1}^3 \right]\,dx\\ \amp = \int_{x=0}^2 (3x + \frac{27}{3} - x - \frac{1}{3})\,dx\\ \amp = \int_{x=0}^2 2x - \frac{26}{3}\,dx\\ \amp = 2^2 - \frac{26}{3}\cdot 2 = -\frac{40}{3} \end{align*}

  2. By Fubini's Theorem, you should get the same answer by integrating in the opposite order.

Subsection3.2.2Summary

  • We can evaluate the double integral \(\iint_R f(x,y) \, dA\) over a rectangle \(R = [a,b] \times [c,d]\) as an iterated integral in one of two ways:

    • \(\int_a^b \left( \int_c^d f(x,y) \, dy \right) \, dx\text{,}\) or

    • \(\int_c^d \left( \int_a^b f(x,y) \, dx \right) \, dy\text{.}\)

    This process works because each inner integral represents a cross-sectional (signed) area and the outer integral then sums all of the cross-sectional (signed) areas. Fubini's Theorem guarantees that the resulting value is the same, regardless of the order in which we integrate.

Subsection3.2.3Exercises

Evaluate each of the following double or iterated integrals exactly.

  1. \(\int_1^3 \left( \int_2^5 xy \, dy \right) \, dx\)

  2. \(\int_0^{\pi/4} \left( \int_0^{\pi/3} \sin(x) \cos(y) \, dx \right) \, dy\)

  3. \(\int_0^1 \left( \int_0^1 e^{-2x - 3y} \, dy \right) \, dx\)

  4. \(\iint_R \sqrt{2x + 5y} \, dA\text{,}\) where \(R = [0,2]\times[0,3]\text{.}\)

The temperature at any point on a metal plate in the \(xy\) plane is given by \(T(x,y) = 100-4x^2 - y^2\text{,}\) where \(x\) and \(y\) are measured in inches and \(T\) in degrees Celsius. Consider the portion of the plate that lies on the rectangular region \(R = [1,5] \times [3,6]\text{.}\)

  1. Write an iterated integral whose value represents the volume under the surface \(T\) over the rectangle \(R\text{.}\)

  2. Evaluate the iterated integral you determined in (a).

  3. Find the area of the rectangle, \(R\text{.}\)

  4. Determine the exact average temperature, \(T_{\operatorname{AVG}(R)}\text{,}\) over the region \(R\text{.}\)

Consider the box with a sloped top that is given by the following description: the base is the rectangle \(R = [1,4] \times [2,5]\text{,}\) while the top is given by the plane \(z = p(x,y) = 30 - x - 2y\text{.}\)

  1. Write an iterated integral whose value represents the volume under \(p\) over the rectangle \(R\text{.}\)

  2. Evaluate the iterated integral you determined in (a).

  3. What is the exact average value of \(p\) over \(R\text{?}\)

  4. If you wanted to build a rectangular box (with an identical base) that has the same volume as the box with the sloped top described here, how tall would the rectangular box have to be?