Active Calculus · 11: Multiple Integrals

11.3: Double Integrals over General Regions

Motivating Questions
  • How do we define a double integral over a non-rectangular region?
  • What general form does an iterated integral over a non-rectangular region have?

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

∬Rf(x,y) dA=lim⁡m,n→∞∑j=1n∑i=1mf(xij∗,yij∗)⋅ΔA, \iint_R f(x,y) \, dA = \lim_{m,n \to \infty} \sum_{j=1}^n \sum_{i=1}^m f(x_{ij}^*, y_{ij}^*) \cdot \Delta A,

where the notation is as described in Section 11.1. Furthermore, we have seen that we can evaluate a double integral ∬Rf(x,y) dA\iint_R f(x,y) \, dA over RR as an iterated integral of either of the forms

∫ab∫cdf(x,y) dy dx      or      ∫cd∫abf(x,y) dx dy. \int_a^b \int_c^d f(x,y) \, dy \, dx \ \ \ \ \ \text{ or } \ \ \ \ \ \int_c^d \int_a^b f(x,y) \, dx \, dy.

It is natural to wonder how we might define and evaluate a double integral over a non-rectangular region; we explore one such example in the following preview activity.

Preview Activity 11.3.111.3.1

A tetrahedron is a three-dimensional figure with four faces, each of which is a triangle. A picture of the tetrahedron TT with vertices (0,0,0),(0,0,0)\text{,} (1,0,0),(1,0,0)\text{,} (0,1,0),(0,1,0)\text{,} and (0,0,1)(0,0,1) is shown at left in Figure 11.3.111.3.1. If we place one vertex at the origin and let vectors a,\mathbf{a}\text{,} b,\mathbf{b}\text{,} and c\mathbf{c} be determined by the edges of the tetrahedron that have one end at the origin, then a formula that tells us the volume VV of the tetrahedron is

V=16∣a⋅(b×c)∣.(11.3.1) V = \frac{1}{6} \lvert \mathbf{a} \cdot (\mathbf{b} \times \mathbf{c}) \rvert.\tag{11.3.111.3.1}

fig_11_3_tetrahedron.svg fig_11_3_tetra_project.svg

Figure 11.3.111.3.1. Left: The tetrahedron T.T\text{.} Right: Projecting TT onto the xyxy-plane.
  1. Use the formula (11.3.111.3.1) to find the volume of the tetrahedron T.T\text{.}
  2. Instead of memorizing or looking up the formula for the volume of a tetrahedron, we can use a double integral to calculate the volume of the tetrahedron T.T\text{.} To see how, notice that the top face of the tetrahedron TT is the plane whose equation is
    z=1−(x+y). z = 1-(x+y).

    Provided that we can use an iterated integral on a non-rectangular region, the volume of the tetrahedron will be given by an iterated integral of the form

    ∫x=?x=?∫y=?y=?1−(x+y) dy dx. \int_{x=?}^{x=?} \int_{y=?}^{y=?} 1-(x+y) \, dy \, dx.

    The issue that is new here is how we find the limits on the integrals; note that the outer integral's limits are in x,x\text{,} while the inner ones are in y,y\text{,} since we have chosen dA=dy dx.dA = dy \, dx\text{.} To see the domain over which we need to integrate, think of standing way above the tetrahedron looking straight down on it, which means we are projecting the entire tetrahedron onto the xyxy-plane. The resulting domain is the triangular region shown at right in Figure 11.3.111.3.1. Explain why we can represent the triangular region with the inequalities

    0≤y≤1−x    and    0≤x≤1. 0 \leq y \leq 1-x \ \ \ \text{ and } \ \ \ 0 \leq x \leq 1.

    (Hint: Consider the cross sectional slice shown at right in Figure 11.3.111.3.1.)

  3. Explain why it makes sense to now write the volume integral in the form
    ∫x=?x=?∫y=?y=?1−(x+y) dy dx=∫x=0x=1∫y=0y=1−x1−(x+y) dy dx. \int_{x=?}^{x=?} \int_{y=?}^{y=?} 1-(x+y) \, dy \, dx = \int_{x=0}^{x=1} \int_{y=0}^{y=1-x} 1-(x+y) \, dy \, dx.
  4. Use the Fundamental Theorem of Calculus to evaluate the iterated integral
    ∫x=0x=1∫y=0y=1−x1−(x+y) dy dx \int_{x=0}^{x=1} \int_{y=0}^{y=1-x} 1-(x+y) \, dy \, dx

    and compare to your result from part (a). (As with iterated integrals over rectangular regions, start with the inner integral.)

So far, we have learned that a double integral over a rectangular region may be interpreted in one of two ways:

  • ∬Rf(x,y) dA\iint_R f(x,y) \, dA tells us the volume of the solids the graph of ff bounds above the xyxy-plane over the rectangle RR minus the volume of the solids the graph of ff bounds below the xyxy-plane under the rectangle R;R\text{;}
  • 1A(R)∬Rf(x,y) dA,\frac{1}{A(R)} \iint_R f(x,y) \, dA\text{,} where A(R)A(R) is the area of RR tells us the average value of the function ff on R.R\text{.} If f(x,y)≥0f(x, y) \geq 0 on R,R\text{,} we can interpret this average value of ff on RR as the height of the box with base RR that has the same volume as the volume of the surface defined by ff over R.R\text{.}

As we saw in Preview Activity 11.1.1, a function f=f(x,y)f = f(x,y) may be considered over regions other than rectangular ones, and thus we want to understand how to set up and evaluate double integrals over non-rectangular regions. Note that if we can, then the two interpretations of the double integral noted above will naturally extend to solid regions with non-rectangular bases.

So, suppose ff is a continuous function on a closed, bounded domain D.D\text{.} For example, consider DD as the circular domain shown at left in Figure 11.3.211.3.2.

fig_11_3_circular.svgfig_11_3_circular_rect.svg 

Figure 11.3.211.3.2. Left: A non-rectangular domain. Right: Enclosing this domain in a rectangle.

We can enclose DD in a rectangular domain RR as shown at right in Figure 11.3.211.3.2 and extend the function ff to be defined over RR in order to be able to use the definition of the double integral over a rectangle. We extend ff in such a way that its values at the points in RR that are not in DD contribute 0 to the value of the integral. In other words, define a function F=F(x,y)F = F(x, y) on RR as

F(x,y)={f(x,y), if (x,y)∈D,0, if (x,y)∉D. F(x,y) = \begin{cases}f(x,y), & \text{ if } (x,y) \in D, \\[4pt] 0, & \text{ if } (x,y) \not\in D \end{cases} .

We then say that the double integral of ff over DD is the same as the double integral of FF over R,R\text{,} and thus

∬Df(x,y) dA=∬RF(x,y) dA. \iint_D f(x,y) \, dA = \iint_R F(x,y) \, dA.

In practice, we just ignore everything that is in RR but not in D,D\text{,} since these regions contribute 0 to the value of the integral.

Just as with double integrals over rectangles, a double integral over a domain DD can be evaluated as an iterated integral. If the region DD can be described by the inequalities g1(x)≤y≤g2(x)g_1(x) \leq y \leq g_2(x) and a≤x≤b,a \leq x \leq b\text{,} where g1=g1(x)g_1=g_1(x) and g2=g2(x)g_2=g_2(x) are functions of only x,x\text{,} then

∬Df(x,y) dA=∫x=ax=b∫y=g1(x)y=g2(x)f(x,y) dy dx. \iint_D f(x,y) \, dA = \int_{x=a}^{x=b} \int_{y=g_1(x)}^{y=g_2(x)} f(x,y) \, dy \, dx.

Alternatively, if the region DD is described by the inequalities h1(y)≤x≤h2(y)h_1(y) \leq x \leq h_2(y) and c≤y≤d,c \leq y \leq d\text{,} where h1=h1(y)h_1=h_1(y) and h2=h2(y)h_2=h_2(y) are functions of only y,y\text{,} we have

∬Df(x,y) dA=∫y=cy=d∫x=h1(y)x=h2(y)f(x,y) dx dy. \iint_D f(x,y) \, dA = \int_{y=c}^{y=d} \int_{x=h_1(y)}^{x=h_2(y)} f(x,y) \, dx \, dy.

The structure of an iterated integral is of particular note:

In an iterated double integral:

  • the limits on the outer integral must be constants;
  • the limits on the inner integral must be constants or in terms of only the remaining variable — that is, if the inner integral is with respect to y,y\text{,} then its limits may only involve xx and constants, and vice versa.

We next consider a detailed example.

Example 11.3.311.3.3

Let f(x,y)=x2yf(x,y) = x^2y be defined on the triangle DD with vertices (0,0),(0,0)\text{,} (2,0),(2,0)\text{,} and (2,3)(2,3) as shown at left in Figure 11.3.311.3.3.

fig_11_3_triangular.svg fig_11_3_triangular_v.svgfig_11_3_triangular_h.svg 

Figure 11.3.311.3.3. A triangular domain and slices in the yy and xx directions.

To evaluate ∬Df(x,y) dA,\iint_D f(x,y) \, dA\text{,} we must first describe the region DD in terms of the variables xx and y.y\text{.} We take two approaches.

Approach 1: Integrate first with respect to y.y\text{.}

In this case we choose to evaluate the double integral as an iterated integral in the form

∬Dx2y dA=∫x=ax=b∫y=g1(x)y=g2(x)x2y dy dx, \iint_D x^2y \, dA = \int_{x=a}^{x=b} \int_{y=g_1(x)}^{y=g_2(x)} x^2y \, dy \, dx,

and therefore we need to describe DD in terms of inequalities

g1(x)≤y≤g2(x)      and      a≤x≤b. g_1(x) \leq y \leq g_2(x) \ \ \ \ \ \text{ and } \ \ \ \ \ a \leq x \leq b.

Since we are integrating with respect to yy first, the iterated integral has the form

∬Dx2y dA=∫x=ax=bA(x) dx, \iint_D x^2y \, dA =\int_{x=a}^{x=b} A(x) \, dx,

where A(x)A(x) is a cross sectional area in the yy direction. So we are slicing the domain perpendicular to the xx-axis and want to understand what a cross sectional area of the overall solid will look like. Several slices of the domain are shown in the middle image in Figure 11.3.411.3.4. On a slice with fixed xx value, the yy values are bounded below by 0 and above by the yy coordinate on the hypotenuse of the right triangle. Thus, g1(x)=0;g_1(x) = 0\text{;} to find y=g2(x),y = g_2(x)\text{,} we need to write the hypotenuse as a function of x.x\text{.} The hypotenuse connects the points (0,0) and (2,3) and hence has equation y=32x.y = \frac{3}{2}x\text{.} This gives the upper bound on yy as g2(x)=32x.g_2(x) = \frac{3}{2}x\text{.} The leftmost vertical cross section is at x=0x=0 and the rightmost one is at x=2,x=2\text{,} so we have a=0a=0 and b=2.b=2\text{.} Therefore,

∬Dx2y dA=∫x=0x=2∫y=0y=32xx2y dy dx. \iint_D x^2y \, dA = \int_{x=0}^{x=2} \int_{y=0}^{y = \frac32 x} x^2y \, dy \, dx.

We evaluate the iterated integral by applying the Fundamental Theorem of Calculus first to the inner integral, and then to the outer one, and find that

∫x=0x=2∫y=0y=32xx2y dy dx=∫x=0x=2[x2⋅y22]∣y=0y=32x dx=∫x=0x=298x4 dx=98x55∣x=0x=2=(98)(325)=365.\begin{align*} \int_{x=0}^{x=2} \int_{y=0}^{y=\frac32 x} x^2y \, dy \, dx & = \int_{x=0}^{x=2} \left[x^2 \cdot \frac{y^2}{2}\right]\biggm|_{y=0}^{y=\frac32 x} \, dx\\[4pt] & = \int_{x=0}^{x=2} \frac{9}{8}x^4\, dx\\[4pt] & = \frac{9}{8}\frac{x^5}{5}\biggm|_{x=0}^{x=2}\\[4pt] & = \left(\frac{9}{8}\right) \left(\frac{32}{5}\right)\\[4pt] & = \frac{36}{5}. \end{align*}
Approach 2: Integrate first with respect to x.x\text{.}

In this case, we choose to evaluate the double integral as an iterated integral in the form

∬Dx2y dA=∫y=cy=d∫x=h1(y)x=h2(y)x2y dx dy \iint_D x^2y \, dA = \int_{y=c}^{y=d} \int_{x=h_1(y)}^{x=h_2(y)} x^2y \, dx \, dy

and thus need to describe DD in terms of inequalities

h1(y)≤x≤h2(y)      and      c≤y≤d. h_1(y) \leq x \leq h_2(y) \ \ \ \ \ \text{ and } \ \ \ \ \ c \leq y \leq d.

Since we are integrating with respect to xx first, the iterated integral has the form

∬Dx2y dA=∫cdA(y) dy, \iint_D x^2y \, dA = \int_c^d A(y) \, dy,

where A(y)A(y) is a cross sectional area of the solid in the xx direction. Several slices of the domain — perpendicular to the yy-axis — are shown at right in Figure 11.3.311.3.3. On a slice with fixed yy value, the xx values are bounded below by the xx coordinate on the hypotenuse of the right triangle and above by 2. So h2(y)=2;h_2(y) = 2\text{;} to find h1(y),h_1(y)\text{,} we need to write the hypotenuse as a function of y.y\text{.} Solving the earlier equation we have for the hypotenuse (y=32xy = \frac32 x) for xx gives us x=23y.x = \frac{2}{3}y\text{.} This makes h1(y)=23y.h_1(y) = \frac{2}{3}y\text{.} The lowest horizontal cross section is at y=0y=0 and the uppermost one is at y=3,y=3\text{,} so we have c=0c=0 and d=3.d=3\text{.} Therefore,

∬Dx2y dA=∫y=0y=3∫x=(2/3)yx=2x2y dx dy. \iint_D x^2y \, dA = \int_{y=0}^{y=3} \int_{x=(2/3)y}^{x=2} x^2y \, dx \, dy.

We evaluate the resulting iterated integral as before by twice applying the Fundamental Theorem of Calculus, and find that

∫y=0y=3∫x=23y2x2y dx dy=∫y=0y=3[x33]∣x=23yx=2y dx=∫y=0y=3[83y−881y4] dy=[83y22−881y55]∣y=0y=3=(83)(92)−(881)(2435)=12−245=365.\begin{align*} \int_{y=0}^{y=3} \int_{x=\frac{2}{3}y}^{2} x^2y \, dx \, dy & = \int_{y=0}^{y=3} \left[\frac{x^3}{3}\right]\biggm|_{x=\frac{2}{3}y}^{x=2}y \, dx\\[4pt] & = \int_{y=0}^{y=3} \left[\frac{8}{3}y - \frac{8}{81}y^4 \right] \, dy\\[4pt] & = \left[\frac{8}{3}\frac{y^2}{2} - \frac{8}{81}\frac{y^5}{5}\right]\biggm|_{y=0}^{y=3}\\[4pt] & = \left(\frac{8}{3}\right) \left(\frac{9}{2}\right) - \left(\frac{8}{81}\right) \left(\frac{243}{5}\right)\\[4pt] & = 12 - \frac{24}{5}\\[4pt] & = \frac{36}{5}. \end{align*}

We see, of course, that in the situation where DD can be described in two different ways, the order in which we choose to set up and evaluate the double integral doesn't matter, and the same value results in either case.

The meaning of a double integral over a non-rectangular region, D,D\text{,} parallels the meaning over a rectangular region. In particular,

  • ∬Df(x,y) dA\iint_D f(x,y) \, dA tells us the volume of the solids the graph of ff bounds above the xyxy-plane over the closed, bounded region DD minus the volume of the solids the graph of ff bounds below the xyxy-plane under the region D;D\text{;}
  • 1A(D)∬Rf(x,y) dA,\frac{1}{A(D)} \iint_R f(x,y) \, dA\text{,} where A(D)A(D) is the area of DD tells us the average value of the function ff on D.D\text{.} If f(x,y)≥0f(x, y) \geq 0 on D,D\text{,} we can interpret this average value of ff on DD as the height of the solid with base DD and constant cross-sectional area DD that has the same volume as the volume of the surface defined by ff over D.D\text{.}
Activity 11.3.211.3.2

Consider the double integral ∬D(4−x−2y) dA,\iint_D (4-x-2y) \, dA\text{,} where DD is the triangular region with vertices (0,0), (4,0), and (0,2).

  1. Write the given integral as an iterated integral of the form ∬D(4−x−2y) dy dx.\iint_D (4-x-2y) \, dy \, dx\text{.} Draw a labeled picture of DD with relevant cross sections.
  2. Write the given integral as an iterated integral of the form ∬D(4−x−2y) dx dy.\iint_D (4-x-2y) \, dx \, dy\text{.} Draw a labeled picture of DD with relevant cross sections.
  3. Evaluate the two iterated integrals from (a) and (b), and verify that they produce the same value. Give at least one interpretation of the meaning of your result.
Activity 11.3.311.3.3

Consider the iterated integral ∫x=0x=1∫y=xy=x(4x+10y) dy dx.\int_{x=0}^{x=1} \int_{y=x}^{y=\sqrt{x}} (4x+10y) \, dy \, dx\text{.}

  1. Sketch the region of integration, D,D\text{,} for which
    ∬D(4x+10y) dA=∫x=0x=1∫y=xy=x(4x+10y) dy dx. \iint_D (4x + 10y) \, dA = \int_{x=0}^{x=1} \int_{y=x}^{y=\sqrt{x}} (4x+10y) \, dy \, dx.
  2. Determine the equivalent iterated integral that results from integrating in the opposite order (dx dy,dx \, dy\text{,} instead of dy dxdy \, dx). That is, determine the limits of integration for which
    ∬D(4x+10y) dA=∫y=?y=?∫x=?x=?(4x+10y) dx dy. \iint_D (4x + 10y) \, dA = \int_{y=?}^{y=?} \int_{x=?}^{x=?} (4x+10y) \, dx \, dy.
  3. Evaluate one of the two iterated integrals above. Explain what the value you obtained tells you.
  4. Set up and evaluate a single definite integral to determine the exact area of D,D\text{,} A(D).A(D)\text{.}
  5. Determine the exact average value of f(x,y)=4x+10yf(x,y) = 4x + 10y over D.D\text{.}
Activity 11.3.411.3.4

Consider the iterated integral ∫x=0x=4∫y=x/2y=2ey2 dy dx.\int_{x=0}^{x=4} \int_{y=x/2}^{y=2} e^{y^2} \, dy \, dx\text{.}

  1. Explain why we cannot find a simple antiderivative for ey2e^{y^2} with respect to y,y\text{,} and thus are unable to evaluate ∫x=0x=4∫y=x/2y=2ey2 dy dx\int_{x=0}^{x=4} \int_{y=x/2}^{y=2} e^{y^2} \, dy \, dx in the indicated order using the Fundamental Theorem of Calculus.
  2. Given that ∬Dey2 dA=∫x=0x=4∫y=x/2y=2ey2 dy dx,\iint_D e^{y^2} \, dA = \int_{x=0}^{x=4} \int_{y=x/2}^{y=2} e^{y^2} \, dy \, dx\text{,} sketch the region of integration, D.D\text{.}
  3. Rewrite the given iterated integral in the opposite order, using dA=dx dy.dA = dx \, dy\text{.} (Hint: You may need more than one integral.)
  4. Use the Fundamental Theorem of Calculus to evaluate the iterated integral you developed in (c). Write one sentence to explain the meaning of the value you found.
  5. What is the important lesson this activity offers regarding the order in which we set up an iterated integral?

Summary

  • For a double integral ∬Df(x,y) dA\iint_D f(x,y) \, dA over a non-rectangular region D,D\text{,} we enclose DD in a rectangle RR and then extend integrand ff to a function FF so that F(x,y)=0F(x,y) = 0 at all points in RR outside of DD and F(x,y)=f(x,y)F(x,y) = f(x,y) for all points in D.D\text{.} We then define ∬Df(x,y) dA\iint_D f(x,y) \, dA to be equal to ∬RF(x,y) dA.\iint_R F(x,y) \, dA\text{.}
  • In an iterated double integral, the limits on the outer integral must be constants while the limits on the inner integral must be constants or in terms of only the remaining variable. In other words, an iterated double integral has one of the following forms (which result in the same value):
    ∫x=ax=b∫y=g1(x)y=g2(x)f(x,y) dy dx, \int_{x=a}^{x=b} \int_{y=g_1(x)}^{y=g_2(x)} f(x,y) \, dy \, dx,

    where g1=g1(x)g_1=g_1(x) and g2=g2(x)g_2=g_2(x) are functions of xx only and the region DD is described by the inequalities g1(x)≤y≤g2(x)g_1(x) \leq y \leq g_2(x) and a≤x≤ba \leq x \leq b or

    ∫y=cy=d∫x=h1(y)x=h2(y)f(x,y) dx dy, \int_{y=c}^{y=d} \int_{x=h_1(y)}^{x=h_2(y)} f(x,y) \, dx \, dy,

    where h1=h1(y)h_1=h_1(y) and h2=h2(y)h_2=h_2(y) are functions of yy only and the region DD is described by the inequalities h1(y)≤x≤h2(y)h_1(y) \leq x \leq h_2(y) and c≤y≤d.c \leq y \leq d\text{.}

Credit and licence

This page titled 11.3: Double Integrals over General Regions was last modified on Fri, 29 Sep 2023 14:34:47 GMT and is shared under a CC BY-SA 4.0 license and was authored, remixed, and/or curated by Matthew Boelkins, David Austin, and Steven Schlicker ( ScholarWorks @Grand Valley State University ) via source content that was edited to the style and standards of the LibreTexts platform.

Source: https://math.libretexts.org/Bookshelves/Calculus/Active_Calculus_(Boelkins_et_al.)/11%3A_Multiple_Integrals/11.03%3A_Double_Integrals_over_General_Regions

Licence: CC BY-SA 4.0 (creativecommons.org/licenses/by-sa/4.0/). You may share and adapt it under the same licence, with this credit. WorldClass adapted the formatting only: the text is the authors', with interactive parts and links to other sites left out.

Build status · Flash Infinity AI™ is a trademark of AI-Flash Inc.. WorldClass is not affiliated with, approved or accredited by any exam board, curriculum provider or university unless a page says so.