11.2: Iterated Integrals

Motivating Questions

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

Recall that we defined the double integral of a continous function \(f=f(x,y)\) over a rectanlge \(R = [a,b] \times [c,d]\) as

\[ \iint_R f(x,y) \, dA = \lim_{m,n \to \infty} \sum_{j=1}^n \sum_{i=1}^m f\left(x_{ij}^*, y_{ij}^*\right) \cdot \Delta A, \nonumber \]

where the different variables and notation are as described in Section 11.1. Thus \(\iint_R f(x,y) \, dA\) is 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 Figure \(\PageIndex{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.

Figure \(\PageIndex{1}\). A solid of revolution.

Preview Activity \(\PageIndex{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
   \[ A(x) = \int_{-4}^4 f(x,y) \, dy. \nonumber \]

   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 Figure \(\PageIndex{2}\).)

   

   

   

   Figure \(\PageIndex{2}\). Left: 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 Figure \(\PageIndex{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{.}\)

Iterated Integrals

The ideas that we explored in Preview Activity \(\PageIndex{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

\[ A(x) = \int_c^d f(x,y) \, dy. \nonumber \]

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.

  

Figure \(\PageIndex{3}\). Summing 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

\[ \int_a^b A(x) \, dx \approx \sum_{i=1}^m A(x_i^*) \Delta x, \nonumber \]

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 Figure \(\PageIndex{3}\). We add the signed volumes of these slices as shown in the frames in Figure \(\PageIndex{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

\[ \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. \nonumber \]

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

\[ \int_a^b \left( \int_c^d f(x,y) \, dy \right) \, dx = \int_a^b \int_c^d f(x,y) \, dy \, dx \nonumber \]

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

\[ \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. \nonumber \]

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

\[ \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. \nonumber \]

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.

Activity \(\PageIndex{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
   \[ A(x) = \int_{-4}^4 f(x,y) \, dy. \nonumber \]

   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?

Activity 4

Let \(f(x,y) = x+y^2\) on the rectangle \(R = [0,2] \times [0,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).

Summary

• 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.