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.4Applications of Double Integrals

Objectives

  • If we have a mass density function for a lamina (thin plate), how does a double integral determine the mass of the lamina?

  • How may a double integral be used to find the area between two curves?

  • Given a mass density function on a lamina, how can we find the lamina's center of mass?

  • What is a joint probability density function? How do we determine the probability of an event if we know a probability density function?

So far, we have interpreted the double integral of a function \(f\) over a domain \(D\) in two different ways. First, \(\iint_D f(x,y) \, dA\) tells us a difference of volumes the volume the surface defined by \(f\) bounds above the \(xy\)-plane on \(D\) minus the volume the surface bounds below the \(xy\)-plane on \(D\text{.}\) In addition, \(\frac{1}{A(D)} \iint_D f(x,y) \, dA\) determines the average value of \(f\) on \(D\text{.}\) In this section, we investigate several other applications of double integrals, using the integration process as seen in Preview Activity3.4.1: we partition into small regions, approximate the desired quantity on each small region, then use the integral to sum these values exactly in the limit.

The following preview activity explores how a double integral can be used to determine the density of a thin plate with a mass density distribution. Recall that in single-variable calculus, we considered a similar problem and computed the mass of a one-dimensional rod with a mass-density distribution. There, as here, the key idea is that if density is constant, mass is the product of density and volume.

Exploration3.4.1

Suppose that we have a flat, thin object (called a lamina) whose density varies across the object. We can think of the density on a lamina as a measure of mass per unit area. As an example, consider a circular plate \(D\) of radius 1 cm centered at the origin whose density \(\delta\) varies depending on the distance from its center so that the density in grams per square centimeter at point \((x, y)\) is

\begin{equation*} \delta(x,y) = 10-2(x^2+y^2). \end{equation*}

  1. Suppose that we partition the plate into subrectangles \(R_{ij}\text{,}\) where \(1 \leq i \leq m\) and \(1 \leq j \leq n\text{,}\) of equal area \(\Delta A\text{,}\) and select a point \((x_{ij}^*,y_{ij}^*)\) in \(R_{ij}\) for each \(i\) and \(j\text{.}\) What is the meaning of the quantity \(\delta(x_{ij}^*,y_{ij}^*) \Delta A\text{?}\)

  2. State a double Riemann sum that provides an approximation of the mass of the plate.

  3. Explain why the double integral

    \begin{equation*} \iint_D \delta(x,y) \, dA \end{equation*}

    tells us the exact mass of the plate.

  4. Determine an iterated integral which, if evaluated, would give the exact mass of the plate. Do not actually evaluate the integral. (This integral is considerably easier to evaluate in polar coordinates, which we will learn more about in Section3.5.)

Subsection3.4.1Mass

Density is a measure of some quantity per unit area or volume. For example, we can measure the human population density of some region as the number of humans in that region divided by the area of that region. In physics, the mass density of an object is the mass of the object per unit area or volume. As suggested by Preview Activity3.4.1, the following holds in general.

The mass of a lamina

If \(\delta(x, y)\) describes the density of a lamina defined by a planar region \(D\text{,}\) then the mass of \(D\) is given by the double integral \(\iint_D \delta(x,y) \, dA\text{.}\)

Example3.4.1

Let \(D\) be a half-disk lamina of radius 3 in quadrants IV and I, centered at the origin as shown in Figure3.4.2. Assume the density at point \((x,y)\) is given by \(\delta(x,y) = x\text{.}\) Find the exact mass of the lamina.

Figure3.4.2A half disk lamina.
Solution

We integrate the density function over the region \(D\text{.}\) Notice that \(D\) can be described by the inequalities \(-3 \le y \le 3\text{,}\) \(0 \le x \le \sqrt{9 - y^2}\text{.}\)

\begin{align*} \iint_D x\,dA \amp = \int_{y=-3}^3 \int_{x=0}^{\sqrt{9-y^2}} x\,dydx \\ \amp = \int_{y=-3}^3 \left[ \frac{(\sqrt{9-y^2})^2}{2} \right]\,dy\\ \amp = \frac12 \int_{y=-3}^3 9 - y^2\,dy\\ \amp = 28 - 18 = 10 \end{align*}

Subsection3.4.2Area

If we consider the situation where the mass-density distribution is constant, we can also see how a double integral may be used to determine the area of a region. Assuming that \(\delta(x,y) = 1\) over a closed bounded region \(D\text{,}\) where the units of \(\delta\) are mass per unit of area, it follows that \(\iint_D 1 \, dA\) is the mass of the lamina. But since the density is constant, the numerical value of the integral is simply the area.

As the following activity demonstrates, we can also see this fact by considering a three-dimensional solid whose height is always 1.

Example3.4.3

Suppose we want to find the area of the bounded region \(D\) between the curves

\begin{equation*} y = 1-x^2 \ \ \ \ \ \text{ and } \ \ \ \ \ y=x-1. \end{equation*}

A picture of this region is shown in Figure3.4.4.

  1. The volume of a solid with constant height is given by the area of the base times the height. Hence, we may interpret the area of the region \(D\) as the volume of a solid with base \(D\) and of uniform height 1. That is, the area of \(D\) is given by \(\iint_D 1 \, dA\text{.}\) Write an iterated integral whose value is \(\iint_D 1 \, dA\text{.}\) (Hint: Which order of integration might be more efficient? Why?)

    Figure3.4.4The graphs of \(y = 1-x^2\) and \(y=x-1\text{.}\)

  2. Evaluate the iterated integral from (a). What does the result tell you?

Solution

  1. It will be more efficient to express \(D\) as a y-simple region, since as an x-simple region its right boundary would have to be divided into two pieces.

    \begin{equation*} \int_{x=-2}^{1}\int_{y=x-1}^{1-x^2} 1\,dydx \end{equation*}
  2. \begin{align*} \int_{x=-2}^{1}\int_{y=x-1}^{1-x^2} 1\,dydx \amp = \int_{x=-2}^{1} 1-x^2 - (x-1) \,dydx \\ \amp = \int_{x=-2}^{1} 2-x^2 - x \,dydx \\ \amp = \left. 2x - \frac{x^3}{3} - \frac{x^2}{2} \right|_{-2}^1 \\ \amp = (2 - \frac{1^3}{3} - \frac{1}{2}) - (-4 - \frac{(-2)^3}{3} - \frac{(-2)^2}{2})\\ \amp = \frac{29}{6} \end{align*}

    The area of the region is \(29/6\text{.}\)

We now formally state the conclusion from our earlier discussion and Example3.4.3.

The double integral and area

Given a closed, bounded region \(D\) in the plane, the area of \(D\text{,}\) denoted \(\operatorname{area}(D)\text{,}\) is given by the double integral

\begin{equation*} \operatorname{area}(D) = \iint_D 1 \, dA. \end{equation*}

Subsection3.4.3Center of Mass

The center of mass of an object is a point at which the object will balance perfectly. For example, the center of mass of a circular disk of uniform density is located at its center. For any object, if we throw it through the air, it will spin around its center of mass and behave as if all the mass is located at the center of mass.

In order to understand the role that integrals play in determining the center of mass of an object with a nonuniform mass distribution, we start by finding the center of mass of a collection of \(N\) distinct point-masses in the plane.

Let \(m_1\text{,}\) \(m_2\text{,}\) \(\ldots\text{,}\) \(m_N\) be \(N\) masses located in the plane. Think of these masses as connected by rigid rods of negligible weight from some central point \((x,y)\text{.}\) A picture with four masses is shown in Figure3.4.5. Now imagine balancing this system by placing it on a thin pole at the point \((x,y)\) perpendicular to the plane containing the masses. Unless the masses are perfectly balanced, the system will fall off the pole. The point \((\overline{x}, \overline{y})\) at which the system will balance perfectly is called the center of mass of the system. Our goal is to determine the center of mass of a system of discrete masses, then extend this to a continuous lamina.

Figure3.4.5A center of mass \((\overline{x}, \overline{y})\) of four masses.

Each mass exerts a force (called a moment) around the lines \(x=\overline{x}\) and \(y=\overline{y}\) that causes the system to tilt in the direction of the mass. These moments are dependent on the mass and the distance from the given line. Let \((x_1,y_1)\) be the location of mass \(m_1\text{,}\) \((x_2,y_2)\) the location of mass \(m_2\text{,}\) etc. In order to balance perfectly, the moments in the \(x\) direction and in the \(y\) direction must be in equilibrium. We determine these moments and solve the resulting system to find the equilibrium point \((\overline{x}, \overline{y})\) at the center of mass.

The force that mass \(m_1\) exerts to tilt the system from the line \(y=\overline{y}\) is

\begin{equation*} m_1g(\overline{y} - y_1), \end{equation*}

where \(g\) is the gravitational constant. Similarly, the force mass \(m_2\) exerts to tilt the system from the line \(y= \overline{y}\) is

\begin{equation*} m_2g(\overline{y}-y_2). \end{equation*}

In general, the force that mass \(m_k\) exerts to tilt the system from the line \(y= \overline{y}\) is

\begin{equation*} m_kg(\overline{y}-y_k). \end{equation*}

For the system to balance, we need the forces to sum to 0, so that

\begin{equation*} \sum_{k=1}^N m_kg(\overline{y}-y_k) = 0. \end{equation*}

Solving for \(\overline{y}\text{,}\) we find that

\begin{equation*} \overline{y} = \frac{\sum_{k=1}^N m_ky_k}{\sum_{k=1}^N m_k}. \end{equation*}

A similar argument shows that

\begin{equation*} \overline{x} = \frac{\sum_{k=1}^N m_kx_k}{\sum_{k=1}^N m_k}. \end{equation*}

The value \(M_x~=~\sum_{k=1}^N m_ky_k\) is called the total moment with respect to the \(x\)-axis; \(M_y~=~\sum_{k=1}^N m_kx_k\) is the total moment with respect to the \(y\)-axis. Hence, the respective quotients of the moments to the total mass, \(M\text{,}\) determines the center of mass of a point-mass system:

\begin{equation*} (\overline{x}, \overline{y}) = \left(\frac{M_y}{M}, \frac{M_x}{M}\right). \end{equation*}

Now, suppose that rather than a point-mass system, we have a continuous lamina with a variable mass-density \(\delta(x, y)\text{.}\) We may estimate its center of mass by partitioning the lamina into \(mn\) subrectangles of equal area \(\Delta A\text{,}\) and treating the resulting partitioned lamina as a point-mass system. In particular, we select a point \((x_{ij}^*,y_{ij}^*)\) in the \(ij\)th subrectangle, and observe that the quanity

\begin{equation*} \delta(x_{ij}^*,y_{ij}^*) \Delta A \end{equation*}

is density times area, so \(\delta(x_{ij}^*,y_{ij}^*) \Delta A\) approximates the mass of the small portion of the lamina determined by the subrectangle \(R_{ij}\text{.}\)

We now treat \(\delta(x_{ij}^*,y_{ij}^*) \Delta A\) as a point mass at the point \((x_{ij}^*,y_{ij}^*)\text{.}\) The coordinates \((\overline{x}, \overline{y})\) of the center of mass of these \(mn\) point masses are thus given by

\begin{equation*} \overline{x} = \frac{\sum_{j=1}^n \sum_{i=1}^m x_{ij}^*\delta(x_{ij}^*,y_{ij}^*) \Delta A}{\sum_{j=1}^n\sum_{i=1}^m \delta(x_{ij}^*,y_{ij}^*) \Delta A} \ \ \ \ \ \text{ and } \ \ \ \ \ \overline{y} = \frac{\sum_{j=1}^n \sum_{i=1}^m y_{ij}^*\delta(x_{ij}^*,y_{ij}^*) \Delta A}{\sum_{j=1}^n\sum_{i=1}^m \delta(x_{ij}^*,y_{ij}^*) \Delta A}. \end{equation*}

If we take the limit as \(m\) and \(n\) go to infinity, we obtain the exact center of mass \((\overline{x}, \overline{y})\) of the continuous lamina.

The center of mass of a lamina

The coordinates \((\overline{x}, \overline{y})\) of the center of mass of a lamina \(D\) with density \(\delta = \delta(x,y)\) are given by

\begin{equation*} \overline{x} = \frac{\iint_D x\delta(x,y) \, dA}{\iint_D \delta(x,y) \, dA} \ \ \ \ \ \text{ and } \ \ \ \ \ \overline{y} = \frac{\iint_D y\delta(x,y) \, dA}{\iint_D \delta(x,y) \, dA}. \end{equation*}

The center of mass of a lamina can then be thought of as a weighted average of all of the points in the lamina with the weights given by the density at each point. The centroid of a lamina is the just the average of all of the points in the lamina, or the center of mass if the density at each point is 1.

The numerators of \(\overline{x}\) and \(\overline{y}\) are called the respective moments of the lamina about the coordinate axes. Thus, the moment of a lamina \(D\) with density \(\delta = \delta(x,y)\) about the \(y\)-axis is

\begin{equation*} M_y = \iint_D x\delta(x,y) \, dA \end{equation*}

and the moment of \(D\) about the \(x\)-axis is

\begin{equation*} M_x = \iint_D y\delta(x,y) \, dA. \end{equation*}

If \(M\) is the mass of the lamina, it follows that the center of mass is

\begin{equation*} (\overline{x}, \overline{y}) = \left(\frac{M_y}{M}, \frac{M_x}{M}\right)\text{.} \end{equation*}
Example3.4.6

A triangular lamina \(D\) with vertices \((0,0)\text{,}\) \((0,3)\text{,}\) and \((3,0)\) has density function \(\delta(x,y) = xy\text{.}\) Find its center of mass.

Solution
\begin{align*} \iint_D x\delta(x,y) \,dA\amp = \int_{x=0}^3 \int_{y=0}^{3-x} x^2y\,dydx = \frac{81}{20} \\ \iint_D y\delta(x,y) \, dA\amp = \int_{x=0}^3 \int_{y=0}^{3-x} xy^2\,dydx = \frac{81}{20} \end{align*}

Next we find the mass:

\begin{align*} \iint_{D} \delta(x,y)\,dA \amp = \int_{x=0}^3 \int_{y=0}^{3-x} xy\,dydx \\ \amp = \int_{x=0}^3 x\cdot \frac{(3-x)^2}{2}\,dx\\ \amp = \frac{27}{8} \end{align*}

Now we are in a position to find the coordinates \(\bar x\) and \(\bar y\) of the center of mass:

\begin{equation*} \bar x = \bar y = \frac{81/20}{27/8} = \frac{6}{5} \end{equation*}

So the center of mass is \((\frac{6}{5},\frac{6}{5})\text{.}\) Contrast this with the centroid of the triangle at \((1,1)\text{.}\)

Subsection3.4.4Probability

Calculating probabilities is a very important application of integration in the physical, social, and life sciences. To understand the basics, consider the game of darts in which a player throws a dart at a board and tries to hit a particular target. Let us suppose that a dart board is in the form of a disk \(D\) with radius 10 inches. If we assume that a player throws a dart at random, and is not aiming at any particular point, then it is equally probable that the dart will strike any single point on the board. For instance, the probability that the dart will strike a particular 1 square inch region is \(\frac{1}{100 \pi}\text{,}\) or the ratio of the area of the desired target to the total area of \(D\) (assuming that the dart thrower always hits the board itself at some point). Similarly, the probability that the dart strikes a point in the disk \(D_3\) of radius 3 inches is given by the area of \(D_3\) divided by the area of \(D\text{.}\) In other words, the probability that the dart strikes the disk \(D_3\) is

\begin{equation*} \frac{9 \pi}{100\pi} = \iint_{D_3} \frac{1}{100 \pi} \, dA. \end{equation*}

The integrand, \(\frac{1}{100\pi}\text{,}\) may be thought of as a distribution function, describing how the dart strikes are distributed across the board. In this case the distribution function is constant since we are assuming a uniform distribution, but we can easily envision situations where the distribution function varies. For example, if the player is fairly good and is aiming for the bulls eye (the center of \(D\)), then the distribution function \(f\) could be skewed toward the center, say

\begin{equation*} f(x,y) = K e^{-(x^2+y^2)} \end{equation*}

for some constant positive \(K\text{.}\) If we assume that the player is consistent enough so that the dart always strikes the board, then the probability that the dart strikes the board somewhere is 1, and the distribution function \(f\) will have to satisfy1This makes \(K = \frac{1}{\pi\left(1-e^{-100}\right)}\text{,}\) which you can check.

\begin{equation*} \iint_D f(x,y) \, dA = 1. \end{equation*}

For such a function \(f\text{,}\) the probability that the dart strikes in the disk \(D_1\) of radius 1 would be

\begin{equation*} \iint_{D_1} f(x,y) \, dA. \end{equation*}

Indeed, the probability that the dart strikes in any region \(R\) that lies within \(D\) is given by

\begin{equation*} \iint_R f(x,y) \, dA. \end{equation*}

The preceding discussion highlights the general idea behind calculating probabilities. We assume we have a joint probability density function \(f\text{,}\) a function of two independent variables \(x\) and \(y\) defined on a domain \(D\) that satisfies the conditions

  • \(f(x,y) \geq 0\) for all \(x\) and \(y\) in \(D\text{,}\)

  • the probability that \(x\) is between some values \(a\) and \(b\) while \(y\) is between some values \(c\) and \(d\) is given by

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

  • The probability that the point \((x,y)\) is in \(D\) is 1, that is

    \begin{equation} \iint_D f(x,y) \, dA = 1.\label{eq_11_4_pdf}\tag{3.4.1} \end{equation}

Note that it is possible that \(D\) could be an infinite region and the limits on the integral in Equation(3.4.1) could be infinite. When we have such a probability density function \(f=f(x,y)\text{,}\) the probability that the point \((x,y)\) is in some region \(R\) contained in the domain \(D\) (the notation we use here is \(P((x,y)\in R)\)) is determined by

\begin{equation*} P((x,y)\in R) = \iint_R f(x,y) \, dA. \end{equation*}
Activity3.4.2

A firm manufactures smoke detectors. Two components for the detectors come from different suppliers one in Michigan and one in Ohio. The company studies these components for their reliability and their data suggests that if \(x\) is the life span (in years) of a randomly chosen component from the Michigan supplier and \(y\) the life span (in years) of a randomly chosen component from the Ohio supplier, then the joint probability density function \(f\) might be given by

\begin{equation*} f(x,y) = e^{-x} e^{-y}. \end{equation*}

  1. Theoretically, the components might last forever, so the domain \(D\) of the function \(f\) is the set \(D\) of all \((x,y)\) such that \(x \ge 0\) and \(y \ge 0\text{.}\) To show that \(f\) is a probability density function on \(D\) we need to demonstrate that

    \begin{equation*} \int \int_D f(x,y) \, dA = 1, \end{equation*}

    or that

    \begin{equation*} \int_0^{\infty} \int_0^{\infty} f(x,y) \, dy \, dx = 1. \end{equation*}

    Use your knowledge of improper integrals to verify that \(f\) is indeed a probability density function.

  2. Assume that the smoke detector fails only if both of the supplied components fail. To determine the probability that a randomly selected detector will fail within one year, we will need to determine the probability that the life span of each component is between 0 and 1 years. Set up an appropriate iterated integral, and evaluate the integral to determine the probability.

  3. What is the probability that a randomly chosen smoke detector will fail between years 3 and 7?

  4. Suppose that the manufacturer determines that one of the components is more likely to fail than the other, and hence conjectures that the probability density function is instead \(f(x,y) = K e^{-x} e^{-2y}.\) What is the value of \(K\text{?}\)

Subsection3.4.5Summary

  • The mass of a lamina \(D\) with a mass density function \(\delta = \delta(x,y)\) is \(\iint_D \delta(x,y) \, dA.\)

  • The area of a region \(D\) in the plane has the same numerical value as the volume of a solid of uniform height 1 and base \(D\text{,}\) so the area of \(D\) is given by \(\iint_D 1 \, dA.\)

  • The center of mass, \((\overline{x},\overline{y})\text{,}\) of a continuous lamina with a variable density \(\delta(x,y)\) is given by

    \begin{equation*} \overline{x} = \frac{\iint_D x\delta(x,y) \, dA}{\iint_D \delta(x,y) \, dA} \ \ \ \ \ \text{ and } \ \ \ \ \ \overline{y} = \frac{\iint_D y\delta(x,y) \, dA}{\iint_D \delta(x,y) \, dA}. \end{equation*}

  • Given a joint probability density function \(f\) is a function of two independent variables \(x\) and \(y\) defined on a domain \(D\text{,}\) if \(R\) is some subregion of \(D\text{,}\) then the probability that \((x,y)\) is in \(R\) is given by

    \begin{equation*} \iint_R f(x,y) \, dA. \end{equation*}

Subsection3.4.6Exercises

A triangular plate is bounded by the graphs of the equations \(y = 2x\text{,}\) \(y = 4x\text{,}\) and \(y = 4\text{.}\) The plate's density at \((x,y)\) is given by \(\delta(x,y) = 4xy^2 + 1\text{,}\) measured in grams per square centimeter (and \(x\) and \(y\) are measured in centimeters).

  1. Set up an iterated integral whose value is the mass of the plate. Include a labeled sketch of the region of integration. Why did you choose the order of integration you did?

  2. Determine the mass of the plate.

  3. Determine the exact center of mass of the plate. Draw and label the point you find on your sketch from (a).

  4. What is the average density of the plate? Include units on your answer.

Let \(D\) be a half-disk lamina of radius 3 in quadrants IV and I, centered at the origin as in Activity3.4.1. Assume the density at point \((x,y)\) is equal to \(x\text{.}\)

  1. Before doing any calculations, what do you expect the \(y\)-coordinate of the center of mass to be? Why?

  2. Set up iterated integral expressions which, if evaluated, will determine the exact center of mass of the lamina.

  3. Use appropriate technology to evaluate the integrals to find the center of mass numerically.

Let \(x\) denote the time (in minutes) that a person spends waiting in a checkout line at a grocery store and \(y\) the time (in minutes) that it takes to check out. Suppose the joint probability density for \(x\) and \(y\) is

\begin{equation*} f(x,y) = \frac{1}{8} e^{-x/4-y/2}. \end{equation*}

  1. What is the exact probability that a person spends between 0 to 5 minutes waiting in line, and then 0 to 5 minutes waiting to check out?

  2. Set up, but do not evaluate, an iterated integral whose value determines the exact probability that a person spends at most 10 minutes total both waiting in line and checking out at this grocery store.

  3. Set up, but do not evaluate, an iterated integral expression whose value determines the exact probability that a person spends at least 10 minutes total both waiting in line and checking out, but not more than 20 minutes.