$\newcommand{\dollar}{\} \DeclareMathOperator{\erf}{erf} \DeclareMathOperator{\arctanh}{arctanh} \newcommand{\lt}{<} \newcommand{\gt}{>} \newcommand{\amp}{&}$

## Section6.4Density, Mass, and Center of Mass

###### Motivating Questions
• How are mass, density, and volume related?

• How is the mass of an object with varying density computed?

• What is is the center of mass of an object, and how are definite integrals used to compute it?

Studying the units of the integrand and variable of integration helps us understand the meaning of a definite integral. For instance, if $v(t)$ is the velocity of an object moving along an axis, measured in feet per second, and $t$ measures time in seconds, then both the definite integral and its Riemann sum approximation,

\begin{equation*} \int_a^b v(t) \, dt \approx \sum_{i=1}^n v(t_i) \Delta t\text{,} \end{equation*}

have units given by the product of the units of $v(t)$ and $t\text{:}$

\begin{equation*} \text{(feet/sec)} \cdot \text{(sec)} = \text{feet}\text{.} \end{equation*}

Thus, $\int_a^b v(t) \, dt$ measures the total change in position of the moving object in feet.

Unit analysis will be particularly helpful to us in what follows.

###### Example6.48

In each of the following scenarios, we consider the distribution of a quantity along an axis.

1. Suppose that the function $c(x) = 200 + 100 e^{-0.1x}$ models the density of traffic on a straight road, measured in cars per mile, where $x$ is the number of miles east of a major interchange, and consider the definite integral

\begin{equation*} \int_0^2 (200 + 100 e^{-0.1x}) \, dx\text{.} \end{equation*}
1. What are the units of the product $c(x) \cdot \Delta x\text{?}$

2. What are the units of the definite integral and its Riemann sum approximation given by

\begin{equation*} \int_0^2 c(x) \, dx \approx \sum_{i=1}^n c(x_i) \Delta x? \end{equation*}
3. Evaluate the definite integral $\int_0^2 c(x) \, dx = \int_0^2 \left(200 + 100 e^{-0.1x}\right) \, dx$ and write one sentence to explain the meaning of the value you find.

2. On a 6 foot long shelf filled with books, the function $B$ models the distribution of the weight of the books, in pounds per inch, where $x$ is the number of inches from the left end of the bookshelf. Let $B(x)$ be given by the rule $B(x) = 0.5 + \frac{1}{(x+1)^2}\text{.}$

1. What are the units of the product $B(x) \cdot \Delta x\text{?}$

2. What are the units of the definite integral and its Riemann sum approximation given by

\begin{equation*} \int_{12}^{36} B(x) \, dx \approx \sum_{i=1}^n B(x_i) \Delta x? \end{equation*}
3. Evaluate the definite integral $\int_{0}^{72} B(x) \, dx = \int_0^{72} \left(0.5 + \frac{1}{(x+1)^2}\right) \, dx$ and write one sentence to explain the meaning of the value you find.

Solution
1. Since $c(x)$ has units of cars per mile, and $\Delta x$ has units of miles, then $c(x) \cdot \Delta x$ has units of cars.

2. From the answer to the previous part, $c(x) \cdot \Delta x$ has units of cars, so the summation and integral will be the sum of things with units of cars, and therefore also have units of cars.

3. Using the exponent rule for integrating:

\begin{equation*} \int_0^2 c(x) \, dx = \int_0^2 \left(200 + 100 e^{-0.1x}\right) \, dx = \left. 200x -1000e^{-0.1x} \right |_0^2\\ = (200(2) - 1000e^{-0.1(2)}) - (200(0)-1000e^{-0.1(0)})=581 \end{equation*}

This represents that, on the two mile stretch of road, there are 581 cars.

1. Since $B(x)$ has units of pounds per inch, and $\Delta x$ has units of inches, then $B(x) \cdot \Delta x$ has units of pounds.

2. From the answer to the previous part, $B(x) \cdot \Delta x$ has units of pounds, so the summation and integral will be the sum of things with units of pounds, and therefore also have units of pounds.

3. Using a power rule, we get that

\begin{equation*} \int_{0}^{72} B(x) \, dx = \int_0^{72} \left(0.5 + \frac{1}{(x+1)^2}\right) \, dx\\ = \left. 0.5 x -\frac{1}{x+1} \right |_0^{72} \\ = \left(0.5(72) - \frac{1}{72+1} \right)-\left(0.5(0)- \frac1{0+1} \right)=36.986 \end{equation*}

This represents that, on the 72 inches of book shelf, there are 36.986 pounds of books.

### SubsectionDensity

The mass of an object, typically measured in metric units such as grams or kilograms, is a measure of the amount of material in the object. The density of an object measures the distribution of mass per unit volume. For instance, if a brick has mass 3 kg and volume 0.002 m$^3\text{,}$ then the density of the brick is

\begin{equation*} \frac{3 \mbox{kg} }{0.002 \mbox{m} ^3} = 1500 \frac{\mbox{kg} }{\mbox{m} ^3}\text{.} \end{equation*}

As another example, the mass density of water is 1000 kg/m$^3\text{.}$ Each of these relationships demonstrate the following general principle.

###### Density Formula

For an object of constant density $d\text{,}$ with mass $m$ and volume $V\text{,}$

\begin{equation*} d = \frac{m}{V}, \ \text{or} \ m = d \cdot V\text{.} \end{equation*}

But what happens when the density is not constant?

The formula $m = d \cdot V$ is reminiscent of two other equations that we have used in our work: for a body moving in a fixed direction, distance = rate $\cdot$ time, and, for a rectangle, its area is given by $A = l \cdot w\text{.}$ These formulas hold when the principal quantities involved, such as the rate the body moves and the height of the rectangle, are constant. When these quantities are not constant, we have turned to the definite integral for assistance. By working with small slices on which the quantity of interest (such as velocity) is approximately constant, we can use a definite integral to add up the values on the pieces.

For example, if we have a nonnegative velocity function that is not constant, over a short time interval $\Delta t\text{,}$ we know that the distance traveled is approximately $v(t) \Delta t$ since $v(t)$ is almost constant on a small interval. Similarly, if we are thinking about the area under a nonnegative function $f$ whose value is changing, on a short interval $\Delta x\text{,}$ the area under the curve is approximately the area of the rectangle whose height is $f(x)$ and whose width is $\Delta x\text{:}$ $f(x) \Delta x\text{.}$ Both of these principles are represented visually in Figure6.49. Figure6.49At left, estimating a small amount of distance traveled, $v(t) \Delta t\text{,}$ and at right, a small amount of area under the curve, $f(x) \Delta x\text{.}$

In a similar way, if the density of some object is not constant, we can use a definite integral to compute the overall mass of the object. We will focus on problems where the density varies in only one dimension, say along a single axis.

###### Example6.50

Let's consider a thin bar of length $b$ whose left end is at the origin, where $x = 0\text{,}$ and assume that the bar has constant cross-sectional area of 1 cm$^2\text{.}$ We let $\rho(x)$ represent the mass density function of the bar, measured in grams per cubic centimeter. That is, given a location $x\text{,}$ $\rho(x)$ tells us approximately how much mass will be found in a one-centimeter wide slice of the bar at $x\text{.}$ Figure6.51A thin bar of constant cross-sectional area 1 cm$^2$ with density function $\rho(x)$ g/cm$^3\text{.}$

The volume of a thin slice of the bar of width $\Delta x\text{,}$ as pictured in Figure6.51, is the cross-sectional area times $\Delta x\text{.}$ Since the cross-sections each have constant area 1 cm$^2\text{,}$ it follows that the volume of the slice is $1 \Delta x$ cm$^3\text{.}$ And because mass is the product of density and volume, we see that the mass of this slice is approximately

\begin{equation*} \text{mass} _{\text{slice} } \approx \rho(x) \ \frac{\mbox{g} }{\mbox{cm} ^3} \cdot 1 \Delta x \ \mbox{cm} ^3 = \rho(x) \cdot \Delta x \ \mbox{g}\text{.} \end{equation*}

Therefore, the total mass of the bar is measured by the corresponding Riemann sum (and the integral that it approximates)

\begin{equation*} \sum_{i=1}^n \rho(x_i) \Delta x \approx \int_0^b \rho(x) \, dx=\text{mass}_{\text{total}}\text{.} \end{equation*}

(As usual, the Riemann sum is an approximation, while the integral will be the exact mass.)

In the example above, computing the mass of the rod was made easier by the fact that the rod had constant cross-sectional area. This meant that the volume of any thin slice was approximately

\begin{equation*} 1\text{cm}^{2} \cdot \Delta x\text{ cm} = \Delta x \text{ cm}^{3} \end{equation*}

indepdently of the value of $x \text{,}$ which allowed us to conclude that the mass of a thin slice at $x$ was approximately $\rho(x) \cdot \Delta x \text{ g} \text{.}$ In general, it may be the case that the volume of our thin slices also varies with $x \text{.}$ In these situations, if $A(x)$ gives the area of a cross section at $x \text{,}$ then the volume of a thin slice at $x$ is approximately $A(x) \Delta x$ and thus the mass of the slice will be approximately $\rho(x) \cdot A(x) \Delta x \text{.}$ From this expression, we pass to a definite integral in the usual way to reach the following.

###### Calculating Mass using a Definite Integral

For an object whose mass is distributed along a single axis according to the function $\rho(x)$ and whose cross-sectional area at $x$ is given by the function $A(x) \text{,}$ the total mass $M$ of the object between $x = a$ and $x = b$ is

\begin{equation*} M = \int_{a}^{b} \rho(x) A(x) \, dx \end{equation*}

Before embarking on an example, we require a short aside on units. The formula above implicitly assumes that the object in question is three-dimensional, since its cross sections have area $A(x)$ and therefore its thin slices have volume $A(x) \Delta x \approx A(x)\, dx \text{.}$ However, we can and often will consider mass and density in one- and two-dimensional settings. In these cases, cross sections don't have area, but rather they are points (in the one-dimensional case) or lengths (in the two-dimensional one). To avoid clashing with the usual convention that density is mass per unit volume, or $\frac{m}{V} \text{,}$ we leave the formula above in terms of three-dimensional objects, albeit with the caveat that for some problems terms like "area" and "volume" must be interpreted in context.

###### Example6.52

Consider the triangle below, which has density given by $\rho(x) = 2x\text{ g/cm}^{2}$ for $0 \leq x \leq 4 \text{.}$ That is, the density of the triangle along a vertical slice at $x$ is $2x \text{ g/cm}^{2} \text{.}$ Figure6.53A triangle with density function $\rho(x) = 2x\text{ g/cm}^{2} \text{.}$

We'll compute the total mass $M$ of the triangle using a definite integral. To accomplish this, we first find the height of a vertical cross section at $x$ for $0 \leq x \leq 4 \text{.}$ The equations of the lines representing the top and bottom edges of the triangle are $y = -\frac{x}{2}+2$ and $y = \frac{x}{2}-2 \text{;}$ therefore, the height of a vertical cross section at $x$ is

\begin{equation*} -\frac{x}{2}+2 - \left(\frac{x}{2}-2\right) = -x+4. \end{equation*}

Thus the mass of a thin slice at $x$ is

\begin{equation*} 2x\text{ g/cm}^{2} \cdot (-x+4) \, dx\text { cm}^{2} = 2x(-x+4) \, dx \text{ g} \end{equation*}

It follows that

\begin{equation*} M = \int_{0}^{4}2x(-x+4) \,dx \end{equation*}

which evaluates to $\frac{64}{3}\text{ g} \text{.}$

###### Calculating Mass using a Definite Integral

For an object of constant cross-sectional area whose mass is distributed along a single axis according to the function $\rho(x)$ (whose units are units of mass per unit of length), the total mass, $M\text{,}$ of the object between $x = a$ and $x = b$ is given by

\begin{equation*} M = \int_a^b \rho(x) \, dx\text{.} \end{equation*}
###### Example6.54

Consider the following situations in which mass is distributed in a non-constant manner.

1. Suppose that a thin rod with constant cross-sectional area of 1 cm$^2$ has its mass distributed according to the density function $\rho(x) = 2e^{-0.2x}\text{,}$ where $x$ is the distance in cm from the left end of the rod, and the units of $\rho(x)$ are g/cm. If the rod is 10 cm long, determine the exact mass of the rod.

2. Consider the cone that has a base of radius 4 m and a height of 5 m. Picture the cone with the center of its base at the origin and think of the cone as a solid of revolution.

1. Write and evaluate a definite integral whose value is the volume of the cone.

2. Next, suppose that the cone has uniform density of 800 kg/m$^3\text{.}$ What is the mass of the solid cone?

3. Now suppose that the cone's density is not uniform, but rather that the cone is most dense at its base. In particular, assume that the density of the cone is uniform across cross sections parallel to its base, but that in each such cross section that is a distance $x$ units from the origin, the density of the cross section is given by the function $\rho(x) = 400 + \frac{200}{1+x^2}\text{,}$ measured in kg/m$^3\text{.}$ Determine and evaluate a definite integral whose value is the mass of this cone of non-uniform density. Do so by first thinking about the mass of a given slice of the cone $x$ units away from the base; remember that in such a slice, the density will be essentially constant.

3. Let a thin rod of constant cross-sectional area 1 cm$^2$ and length 12 cm have its mass be distributed according to the density function $\rho(x) = \frac{1}{25}(x-15)^2\text{,}$ measured in g/cm. Find the exact location $x$ at which to cut the bar so that the two pieces will each have identical mass.

Hint
1. Remember that $M = \int_a^b \rho(x) \, dx\text{.}$

2. Think about slices: what's the volume of a slice? what's the mass of a slice?

3. Consider $\int_0^b \rho(x) \, dx\text{.}$

1. $M = 0.1 - 0.1e^{-4} \approx 0.0981684$ grams.

1. $V = \int_{0}^{5} \pi (4 - \frac{4}{5}x)^2 \, dx = \frac{80\pi}{3} \approx 83.7758 \mbox{m}^3\text{.}$

2. $M = \int_{0}^{5} (400 + \frac{200}{1+x^2}) \cdot \pi (4-\frac{4}{5}x)^2 \, dx$ $= \frac{32}{5}\pi(25015 - 5 \ln(17576) + 72\tan^{-1}(5) \approx 33597.4664 \mbox{kg}$

2. $b \approx 3.0652\text{.}$

Solution
1. Since the mass, $M\text{,}$ is given by $M = \int_a^b \rho(x) \, dx\text{,}$ it follows that

\begin{equation*} M = \int_{0}^{10} 2e^{-0.2x} = \left. \frac{2}{-0.2}e^{-0.2x} \right|_{0}^{20} = -0.1e^{-4}+0.1e^{0}\text{,} \end{equation*}

and hence $M = 0.1 - 0.1e^{-4} \approx 0.0981684$ grams.

2. Consider the cone that has a base of radius 4 m and a height of 5 m.

1. With the cone having its base centered at the origin and being a solid of revolution about the $x$-axis, the cone is generated by the line $y = 4 - \frac{4}{5}x\text{.}$ Using the disk method, it follows that the cone's volume is

\begin{equation*} V = \int_{0}^{5} \pi (4 - \frac{4}{5}x)^2 \, dx = \frac{80\pi}{3} \approx 83.7758 \mbox{m}^3\text{.} \end{equation*}
2. If the cone has uniform density 800 kg/m$^3\text{,}$ we know that the cone's mass is the product of its density and volume, and thus

\begin{equation*} M = \rho \cdot V = 800 \cdot \frac{80\pi}{3} = \frac{6400\pi}{3} \approx 67020.6433 \mbox{kg}\text{.} \end{equation*}
3. If the cone has non-uniform density, but is uniform across cross sections parallel to its base such that in each cross section a distance $x$ units from the origin has its density given by $\rho(x) = 400 + \frac{200}{1+x^2}$ kg/m$^3\text{,}$ we naturally consider slices of the cone that are perpendicular to the $x$-axis. For each such slice, the radius is $r(x) = 4 - \frac{4}{5}x\text{,}$ which makes the volume of the slice

\begin{equation*} V_{\text{slice}} \approx \pi (4-\frac{4}{5}x)^2 \Delta x\text{.} \end{equation*}

Since the density is constant on the slice with value $\rho(x) = 400 + \frac{200}{1+x^2}\text{,}$ it follows that the mass of the slice is

\begin{equation*} M_{\text{slice}} \approx (400 + \frac{200}{1+x^2}) \cdot \pi (4-\frac{4}{5}x)^2 \Delta x\text{.} \end{equation*}

Using a definite integral to let $\Delta x \to 0$ and add up the values of all of the slices, we thus find the mass of the cone with this density distribution to be

\begin{align*} M \amp = \int_{0}^{5} (400 + \frac{200}{1+x^2}) \cdot \pi (4-\frac{4}{5}x)^2 \, dx\\ \amp = \frac{32}{5}\pi(25015 - 5 \ln(17576) + 72\tan^{-1}(5) \approx 33597.4664 \mbox{kg}\\ \amp \approx 33597.4664 \mbox{kg}\text{.} \end{align*}
3. The total mass of the bar is given by $M = \int_{0}^{12} \frac{1}{25}(x-15)^2 \, dx = \frac{1116}{25}\text{.}$ To find where to cut the bar in two so that the two pieces have equal mass, we have to determine the value of $b$ such that

\begin{equation*} \int_{0}^{b} \frac{1}{25}(x-15)^2 \, dx = \frac{1}{2} \cdot \frac{1116}{25}\text{.} \end{equation*}

Since $\int_{0}^{b} \frac{1}{25}(x-15)^2 \, dx = \frac{1}{75}( (b-15)^3 - (-15)^3 )\text{,}$ we need to have $b$ satisfy the equation

\begin{equation*} \frac{1}{75}( (b-15)^3 - (-15)^3 ) = \frac{1}{2} \cdot \frac{1116}{25}\text{.} \end{equation*}

Thus, $(b-15)^3 - (-15)^3 = 75 \cdot \frac{558}{25} = 3 \cdot 558 = 1674\text{,}$ so $(b-15)^3 = 1675 - 3375 = -1700\text{.}$ It follows that $b-15 = \sqrt{-1700} \approx -11.9348\text{,}$ so $b \approx 3.0652\text{.}$

In some 2D or 3D objects, density depends on distance from a center point. In this case, we will apply the principle from the previous section (that mass is the integral of density), but must make a slight adjustment to account for the fact that layers of the same thickness will have different volumes depending on their distance from the center.

###### Example6.55

Consider a 1 meter long metal pipe covered with a layer of insulation. How would we find the total mass of the pipe? We would need to add up the masses of the metal and the insulation, and these masses can be found by taking density times volume:

\begin{equation*} \text{mass}_{\text{layer}}=\rho_{\text{layer}}V_{\text{layer}} \end{equation*}

But note that the density of the layer depends on its distance from the center: close to the center of the pipe, the hollow core of the pipe has 0 density, at a radius of $r_1$ the metal pipe has a high density, and at a radius of $r_2$ the insulation has a low density. Let us use $\rho(r)$ to denote the density as a function of distance $r$ from the center.

Next, to get $V_{\text{layer}}$ we can use approximate volume as circumference times thickness times length, so $V \approx (2\pi r_i)(\Delta r)(1 \text{ meter})=(2\pi r_i)(\Delta r)\text{.}$ Thus the mass of a single layer is approximately

\begin{equation*} \text{mass}_i \approx \rho(r_i)(2\pi r_i)(\Delta r). \end{equation*}

Summing both layers, we get that

\begin{equation*} \text{mass} \approx \displaystyle \sum_{i=1}^2 \rho(r_i)(2\pi r_i)(\Delta r) \end{equation*}

is the total mass of the pipe.

We can generalize the formula we found in the previous example to any object for which density depends on radius. That more general formula is

\begin{equation*} \text{mass} \approx \displaystyle \sum_{i=1}^n \rho(r_i)(2\pi r_i)(\Delta r) \end{equation*}

where $\rho(r)$ is the density of the object a distance of $r_i$ from the center.

Similar to the earlier applications in this chapter: we will take the limit as $\Delta r \rightarrow 0 \text{,}$ and doing so turns the sum into an integral.

###### Calculating Mass using a Definite Integral for Radially Symmetric Objects

For a 2D disk (respectively, 3D cylinder), if the density of the object depends only on the distance from the center point (respectively, center axis) according to the function $\rho(r)\text{,}$ then the mass per unit length (respectively, mass of the object) is given by the formula

\begin{equation*} \text{mass} = \displaystyle \int_{0}^R \rho(r)(2\pi r)\, dr= \displaystyle \int_{0}^R 2 \pi r\rho(r)\, dr \end{equation*}

where $R$ is the outer radius of the object.

Note that this formula is almost identical to the one we use for calculating mass of an object with constant cross-sectional area, except that this formula incorporates a factor of $2 \pi r$ (the circumference of a circle) inside the integral.

###### Example6.56

After an oil spill, scientists measure that the concentration of oil depends on the distance from the source according to the function $\rho(r)=150 e^{-r^2}$ where $r$ is in miles, and $\rho(r)$ is in barrels per square mile. If the spill has reached a maximum radius of 5 miles, how much oil was spilled?

Hint

Set up the integral like for the mass of a disk, then integrate. To evaluate the integral use the $u$-substitution $u=-r^2\text{.}$

$471.23$ barrels.

Solution

Since the density depends on the radius, we will use the mass of a disk formula $\displaystyle \int_{0}^R \rho(r)(2\pi r)\, dr \text{.}$

For us, the maximum radius is $R=5$ miles, so the total amount of oil is

\begin{equation*} \displaystyle \int_{0}^5 (150e^{-r^2})(2\pi r)\, dr = \left.-150 \pi e^{-r^2} \right|_0^5 \\ (-150\pi e^{-(5)^2}) - (-150\pi e^{-(0)^2})=471.23 \text{ barrels} \end{equation*}
###### Example6.57

The 2D region bounded by $y=6-x\text{,}$ $x=0\text{,}$ $x=3\text{,}$ and $y=0$ (pictured below), is rotated around the $x$-axis to make a solid of revolution. Figure6.58The 2D region before it is rotated around the $x$-axis.

When this object is carved out of material, its density depends on the distance from the $x$-axis according to $\rho_{radial}(r)=9-r\text{.}$

1. Fix a value of $x$ in the interval $[0,3]\text{.}$ Set up and evaluate the integral to find $\rho_{axis}(x)\text{,}$ the mass per unit length of the circular cross-section at that $x$-value. (Your answer will be in terms of $x\text{.}$)

2. Use your answer to the previous question and the techniques in the Calculating Mass using a Definite Integral section to find the total mass of the object.

Hint
1. Each cross-section is a disk with outer radius $y=6-x\text{.}$

2. Use your formula for $\rho_{axis}(x)$ and the formula from "Calculating Mass using a Definite Integral":
\begin{equation*} M=\int_a^b \rho(x)dx \end{equation*}
1. $\rho_{axis}(x)=\pi(\frac{2}{3} x^3 - 3x^2 -36x+180)$

2. $M=364.5 \pi \approx 1145.11$
Solution
1. Since we have circular cross-sections and their density only depends on the distance to the center of the circle (which is the $x$-axis), we will integrate as we do for a disk. Note that our outer radius depends on the $x$-value, and in fact, the outer radius is $y=6-x\text{.}$ Therefore,

\begin{equation*} \rho_{axis}(x)= \int_0^{6-x}(2 \pi r) \rho_{radial}(r) \, dr \\ = \int_0^{6-x}(2 \pi r) (9-r) \, dr \\ = \left. (9 \pi r^2- \frac{2\pi}{3} r^3) \right|_{r=0}^{r=6-x}\\ = (9 \pi (6-x)^2- \frac{2\pi}{3} (6-x)^3)-(0)\\ = \pi(\frac{2}{3} x^3 - 3x^2 -36x+180) \end{equation*}
2. What we found in the previous part is the mass per unit length, so just as in "Calculating Mass using a Definite Integral", we simply need to integrate $\rho_{axis}(x)$ from the endpoints $a=0$ to $b=3\text{:}$

\begin{equation*} M= \int_0^3 \rho_{axis}(x) \, dx =\int_0^3 \pi(\frac{2}{3} x^3 - 3x -36x+180) \, dx \\ = \left. \pi(\frac{1}{6} x^4 - x^3 -18x^2+180x) \right|_{0}^3\\ =\pi(\frac{1}{6} (3)^4 - (3)^3 -18(3)^2+180(3))-(0)\\ = 364.5 \pi \approx 1145.11 \end{equation*}
###### Example6.59

Suppose that a circular city's population density is given by $\rho(r)=1000-100r\text{,}$ where $r$ is the distance from the city's center in miles and $\rho(r)$ is in units of people per square mile. If the city has a radius of 6 miles, find the population of the city.

Hint

Set up the integral like for the mass of a disk, then integrate.

$67858$ people.

Solution

Since the density depends on the radius, we will use the mass of a disk formula $\displaystyle \int_{0}^R \rho(r)(2\pi r)\, dr \text{.}$

For us, the maximum radius is $R=6$ miles, so the total population is

\begin{equation*} \displaystyle \int_{0}^6 ((1000-100r)(2\pi r)\, dr = \left. 1000 \pi r^2-\frac{200\pi}{3}r^3 \right|_0^6 \\ (1000 \pi (6)^2 - \frac{200 \pi}{3} (6)^3) - (0) \approx 67858 \text{ people} \end{equation*}
###### Example6.60

When the density of a spherical object only depends on the distance from the center, the formula for the object's mass is

\begin{equation*} \text{mass}= \int_0^R 4 \pi r^2 \rho(r)\, dr \end{equation*}

which is identical to the formula for a circular object, except that the circumference of a circle ($C=2\pi r$) is replaced by the surface area of a sphere ($A=4\pi r^2$).

The earth is denser near its core. Suppose that the density of the earth was given by $\rho(r)=(-0.0016)r+13$ where $r$ has units of kilometers and $\rho$ is in trillions of kilograms per cubic kilometer (in more normal units, $1 \text{ trillion }\frac{kg}{km^3}=1 \frac{g}{cm^3}$).

1. The earth has a radius of approximately 6300 kilometers. Find the total mass of the earth.

2. The volume of a sphere is $V=\frac 43 \pi R^3\text{.}$ Find the volume of the earth, and calculate the average density of the earth. (Average density is the total mass divided by total volume.)

3. Goofus sees that at the center of the earth, the density is $\rho(0)=13\text{,}$ and at the surface the density is $\rho(6300)=2.92\text{,}$ and therefore concludes that the average density of the earth is $\frac{13+2.92}{2}=7.96\text{.}$ Is Goofus correct? If he is not, explain why his calculation is incorrect.

4. The earth's inner core has a radius of 1200 kilometers. Find the average density of the earth's inner core.

5. What percent of the earth's mass is in the inner core? What percent of the earth's volume is in the inner core?

Hint
1. Set up and evaluate the integral given by the equation

\begin{equation*} \text{mass}= \int_0^R 4 \pi r^2 \rho(r)\, dr \end{equation*}
2. Use the volume formula to get the volume, then divide mass by volume to get the average density.

3. Goofus is incorrect.

4. Perform the same calculations as parts (a.) and (b.) but for the smaller radius.

1. The mass of the earth is roughly $5.7 \times 10^{12}$ trillions of kilograms.

2. The volume of the earth is $V=1.05 \times 10^{12}$ cubic kilometers. Its average density is $5.44$ trillions of kilograms per cubic kilometer.

3. Goofus is incorrect.

4. The density is $11.56$ trillion kilograms per cubic kilometer.

5. Approximately $1.47 \%$ of the earth's mass is in the inner core.

Approximately $0.689 \%$ of the earth's volume is in the inner core.

Solution
1. We must take the integral:

\begin{equation*} \text{mass}= \int_0^R 4 \pi r^2 \rho(r)\, dr = \int_0^{6300} 4 \pi r^2(-0.0016r+13)\, dr \\ = \left. 4\pi(-0.0004 r^4+\frac{13}{3}r^3) \right |_0^{6300}\\ \approx 5.7 \times 10^{12} \end{equation*}

so the mass of the earth is roughly $5.7 \times 10^{12}$ trillions of kilograms.

2. The volume formula for a sphere says that the volume of the earth is $V=\frac{4}{3}\pi (6300)^3=1.05 \times 10^{12}$ cubic kilometers. Therefore, its average density is $\displaystyle \frac{5.7 \times 10^{12}}{1.05 \times 10^{12}}=5.44$ trillions of kilograms per cubic kilometer.

3. Goofus is incorrect. His calculation is wrong because he is not taking into account the factor of $r^2$ inside the integral, which counts the lower-density parts on the outside of the earth more than the higher-density parts near the center. For instance, his average treated innermost 1 kilometer of the earth as having the same volume as the outermost 1 kilometer of earth, even though the latter is much larger in terms of volume.

4. We will essentially repeat our calculations from parts (a.) and (b.), except with 1200 as the maximum radius instead of 6300. Therefore, the volume of the inner core is $V=\frac{4}{3}\pi (1200)^3=7.24 \times 10^{9}$ cubic killometers, and the mass is

\begin{equation*} \text{mass}= \int_0^R 4 \pi r^2 \rho(r)\, dr = \int_0^{1200} 4 \pi r^2(-0.0016r+13)\, dr \\ = \left. 4\pi(-0.0004 r^4+\frac{13}{3}r^3) \right |_0^{1200}\\ \approx 8.37 \times 10^{10} \end{equation*}

so the density is $(8.37 \times 10^{10})/(7.24 \times 10^{9})=11.56$ trillion kilograms per cubic kilometer.

5. In the previous parts, we saw that the mass of the inner core was $8.37 \times 10^{10}$ and the mass of the whole earth was $5.7 \times 10^{12}\text{,}$ so dividing gives that approximately $1.47 \%$ of the earth's mass is in the inner core.

In the previous parts, we saw that the volume of the inner core was $7.24 \times 10^{9}$ and the volume of the whole earth was $1.05 \times 10^{12}\text{,}$ so dividing gives that approximately $0.689 \%$ of the earth's volume is in the inner core.

### SubsectionWeighted Averages

The concept of an average is a natural one, and one that we have used repeatedly as part of our understanding of the meaning of the definite integral. If we have $n$ values $a_1\text{,}$ $a_2\text{,}$ $\ldots\text{,}$ $a_n\text{,}$ we know that their average is given by

\begin{equation*} \frac{a_1 + a_2 + \cdots + a_n}{n}\text{,} \end{equation*}

and for a quantity being measured by a function $f$ on an interval $[a,b]\text{,}$ the average value of the quantity on $[a,b]$ is

\begin{equation*} \frac{1}{b-a} \int_a^b f(x) \, dx\text{.} \end{equation*}

As we continue to think about problems involving the distribution of mass, it is natural to consider the idea of a weighted average, where certain quantities involved are counted more in the average.

A common use of weighted averages is in the computation of a student's GPA, where grades are weighted according to credit hours. Let's consider the scenario in Table6.61.

 class grade grade points credits chemistry B+ 3.3 5 calculus A- 3.7 4 history B- 2.7 3 psychology B- 2.7 3

If all of the classes were of the same weight (i.e., the same number of credits), the student's GPA would simply be calculated by taking the average

\begin{equation*} \frac{3.3 + 3.7 + 2.7 + 2.7}{4} = 3.1\text{.} \end{equation*}

But since the chemistry and calculus courses have higher weights (of 5 and 4 credits respectively), we actually compute the GPA according to the weighted average

\begin{equation*} \frac{3.3 \cdot 5 + 3.7 \cdot 4 + 2.7 \cdot 3 + 2.7 \cdot 3}{5 + 4 + 3 + 3} = 3.1\overline{6}\text{.} \end{equation*}

The weighted average reflects the fact that chemistry and calculus, as courses with higher credits, have a greater impact on the students' grade point average. Note particularly that in the weighted average, each grade gets multiplied by its weight, and we divide by the sum of the weights.

In the following example, we explore further how weighted averages can be used to find the balancing point of a physical system.

###### Example6.62

For quantities of equal weight, such as two children on a teeter-totter, the balancing point is found by taking the average of their locations. When the weights of the quantities differ, we use a weighted average of their respective locations to find the balancing point.

1. Suppose that a shelf is 6 feet long, with its left end situated at $x = 0\text{.}$ If one book of weight 1 lb is placed at $x_1 = 0\text{,}$ and another book of weight 1 lb is placed at $x_2 = 6\text{,}$ what is the location of $\overline{x}\text{,}$ the point at which the shelf would (theoretically) balance on a fulcrum?

2. Now, say that we place four books on the shelf, each weighing 1 lb: at $x_1 = 0\text{,}$ at $x_2 = 2\text{,}$ at $x_3 = 4\text{,}$ and at $x_4 = 6\text{.}$ Find $\overline{x}\text{,}$ the balancing point of the shelf.

3. How does $\overline{x}$ change if we change the location of the third book? Say the locations of the 1-lb books are $x_1 = 0\text{,}$ $x_2 = 2\text{,}$ $x_3 = 3\text{,}$ and $x_4 = 6\text{.}$

4. Next, suppose that we place four books on the shelf, but of varying weights: at $x_1 = 0$ a 2-lb book, at $x_2 = 2$ a 3-lb book, at $x_3 = 4$ a 1-lb book, and at $x_4 = 6$ a 1-lb book. Use a weighted average of the locations to find $\overline{x}\text{,}$ the balancing point of the shelf. How does the balancing point in this scenario compare to that found in (b)?

5. What happens if we change the location of one of the books? Say that we keep everything the same in (d), except that $x_3 = 5\text{.}$ How does $\overline{x}$ change?

6. What happens if we change the weight of one of the books? Say that we keep everything the same in (d), except that the book at $x_3 = 4$ now weighs 2 lbs. How does $\overline{x}$ change?

7. Experiment with a couple of different scenarios of your choosing where you move one of the books to the left, or you decrease the weight of one of the books.

8. Write a couple of sentences to explain how adjusting the location of one of the books or the weight of one of the books affects the location of the balancing point of the shelf. Think carefully here about how your changes should be considered relative to the location of the balancing point $\overline{x}$ of the current scenario.

Hint
1. Find the average location.

2. Note that you are averaging 4 locations.

3. Compare (b).

4. If the book at location $x_1$ weighs 2 pounds, it's like there are two books at $x_1\text{,}$ so $2x_1$ must play a role in the average.

5. Use a weighted average.

6. Compare (e).

7. Experiment.

8. Compare the effects of moving locations left and adding weight to the left of the balancing point.

1. $\overline{x} = \frac{x_1 + x_2}{2} = 3\text{.}$

2. $\overline{x} = \frac{x_1 + x_2 + x_3 + x_4}{4} = 3\text{.}$

3. $\overline{x} = \frac{x_1 + x_2 + x_3 + x_4}{4} = 2.75\text{.}$

4. $\overline{x} = \frac{2x_1 + 3x_2 + 1x_3 + 1x_4}{7} = \frac{16}{7}\text{.}$

5. $\overline{x} = \frac{2x_1 + 3x_2 + 1x_3 + 1x_4}{7} = \frac{17}{7}\text{.}$

6. $\overline{x} = \frac{2x_1 + 3x_2 + 1x_3 + 2x_4}{7} = \frac{22}{7}\text{.}$

8. If we have an existing arrangement and balancing point, moving one of the locations to the left will move the balancing point to the left; similarly, moving one of the locations to the right will move the balancing point to the right. If instead we add weight to an existing location, if that location is left of the balancing point, the balancing point will move left; the behavior is similar if on the right.

Solution
1. The shelf will balance at $\overline{x} = \frac{x_1 + x_2}{2} = \frac{0 + 6}{2} = 3\text{.}$

2. The shelf will balance at the average of the locations: $\overline{x} = \frac{x_1 + x_2 + x_3 + x_4}{4} = \frac{0 + 2 + 4 + 6}{4} = 3\text{.}$

3. Changing the location of $x_3$ to the left will move the average to the left, so the new location for $\overline{x}$ is $\overline{x} = \frac{x_1 + x_2 + x_3 + x_4}{4} = \frac{0 + 2 + 3 + 6}{4} = \frac{11}{4} = 2.75\text{.}$

4. We can think of the weighted books almost as if there were that number of 1-lb books at each location. The balancing point is given by the weighted average

\begin{equation*} \overline{x} = \frac{2x_1 + 3x_2 + 1x_3 + 1x_4}{7} = \frac{2 \cdot 0 + 3 \cdot 2 + 1 \cdot 4 + 1 \cdot 6}{7} = \frac{16}{7}\text{.} \end{equation*}

Note particularly that we divide by the total weight of the books, not the total number.

5. Here we find

\begin{equation*} \overline{x} = \frac{2x_1 + 3x_2 + 1x_3 + 1x_4}{7} = \frac{2 \cdot 0 + 3 \cdot 2 + 1 \cdot 5 + 1 \cdot 6}{7} = \frac{17}{7}\text{,} \end{equation*}

which moves the balancing point slightly to the right.

6. In this situation, we find the balancing point moves significantly to the right:

\begin{equation*} \overline{x} = \frac{2x_1 + 3x_2 + 1x_3 + 2x_4}{7} = \frac{2 \cdot 0 + 3 \cdot 2 + 1 \cdot 4 + 2 \cdot 6}{7} = \frac{22}{7}\text{.} \end{equation*}
7. Try, for instance, having $x_2 = 1$ and $x_3 = 2$ and experiment with different weights.

8. If we have an existing arrangement and balancing point, moving one of the locations to the left will move the balancing point to the left; similarly, moving one of the locations to the right will move the balancing point to the right. If instead we add weight to an existing location, if that location is left of the balancing point, the balancing point will move left; the behavior is similar if on the right.

### SubsectionCenter of Mass

In Example6.62, we saw that the balancing point of a system of point-masses2In the example, we actually used weight rather than mass. Since weight is proportional to mass, the computations for the balancing point result in the same location regardless of whether we use weight or mass. The gravitational constant is present in both the numerator and denominator of the weighted average.(such as books on a shelf) is found by taking a weighted average of their respective locations. In the example, we were computing the center of mass of a system of masses distributed along an axis, which is the balancing point of the axis on which the masses rest.

###### Center of Mass (point-masses)

For a collection of $n$ masses $m_1\text{,}$ $\ldots\text{,}$ $m_n$ that are distributed along a single axis at the locations $x_1\text{,}$ $\ldots\text{,}$ $x_n\text{,}$ the center of mass is given by

\begin{equation*} \overline{x} = \frac{x_1 m_1 + x_2 m_2 + \cdots + x_n m_n}{m_1 + m_2 + \cdots + m_n}\text{.} \end{equation*}

Now consider a thin bar over which density is distributed continuously. If the density is constant, it is obvious that the balancing point of the bar is its midpoint. But if density is not constant, we must compute a weighted average. Let's say that the function $\rho(x)$ tells us the density distribution along the bar, measured in g/cm. If we slice the bar into small sections, we can think of the bar as holding a collection of adjacent point-masses. The mass $m_i$ of a slice of thickness $\Delta x$ at location $x_i\text{,}$ is $m_i \approx \rho(x_i) \Delta x\text{.}$

If we slice the bar into $n$ pieces, we can approximate its center of mass by

\begin{equation*} \overline{x} \approx \frac{x_1 \cdot \rho(x_1) \Delta x + x_2 \cdot \rho(x_2) \Delta x + \cdots + x_n \cdot \rho(x_n) \Delta x }{\rho(x_1) \Delta x + \rho(x_2) \Delta x + \cdots + \rho(x_n) \Delta x}\text{.} \end{equation*}

Rewriting the sums in sigma notation, we have

\begin{equation} \overline{x} \approx \frac{\sum_{i = 1}^{n} x_i \cdot \rho(x_i) \Delta x}{\sum_{i = 1}^{n} \rho(x_i) \Delta x}\text{.}\label{E-CtrOfM}\tag{6.4} \end{equation}

The greater the number of slices, the more accurate our estimate of the balancing point will be. The sums in Equation(6.4) can be viewed as Riemann sums, so in the limit as $n \to \infty\text{,}$ we find that the center of mass is given by the quotient of two integrals.

###### Center of Mass (continuous mass distribution)

For a thin rod of density $\rho(x)$ distributed along an axis from $x = a$ to $x = b\text{,}$ the center of mass of the rod is given by

\begin{equation*} \overline{x} = \frac{\int_a^b x \rho(x) \, dx}{\int_a^b \rho(x) \, dx}\text{.} \end{equation*}

Note that the denominator of $\overline{x}$ is the mass of the bar, and that this quotient of integrals is simply the continuous version of the weighted average of locations, $x\text{,}$ along the bar.

###### Example6.63

Consider a thin bar of length 20 cm whose density is distributed according to the function $\rho(x) = 4 + 0.1x\text{,}$ where $x = 0$ represents the left end of the bar. Assume that $\rho$ is measured in g/cm and $x$ is measured in cm.

1. Find the total mass, $M\text{,}$ of the bar.

2. Without doing any calculations, do you expect the center of mass of the bar to be equal to 10, less than 10, or greater than 10? Why?

3. Compute $\overline{x}\text{,}$ the exact center of mass of the bar.

4. What is the average density of the bar?

5. Now consider a different density function, given by $p(x) = 4e^{0.020732x}\text{,}$ also for a bar of length 20 cm whose left end is at $x = 0\text{.}$ Plot both $\rho(x)$ and $p(x)$ on the same axes. Without doing any calculations, which bar do you expect to have the greater center of mass? Why?

6. Compute the exact center of mass of the bar described in (e) whose density function is $p(x) = 4e^{0.020732x}\text{.}$ Check the result against the prediction you made in (e).

Hint
1. Recall $M = \int_a^b \rho(x) \, dx\text{.}$

2. Think about whether there's more mass to the left of 10 or to the right of 10.

3. Use the formula just developed in the preceding subsection of the text.

4. Recall that when density, $\rho$ is constant, $M = \rho \cdot V\text{.}$

5. Think about how the two density functions distribute mass. Note that they appear to distribute the same total mass, as the areas under the two curves appear to be equal.

6. Use $p(x)$ in place of $\rho(x)$ in standard calculations.

1. $M = \int_{0}^{20} 4 + 0.1x \, dx = 100$ g.

2. Greater than 10.

3. $\overline{x} = \frac{\int_{0}^{20} x (4 + 0.1x)) \, dx}{\int_{0}^{20} 4 + 0.1x \, dx} = \frac{32}{3}\text{.}$

4. 5 g/cm.

5. Slightly to the right of the center of mass for $\rho(x)\text{.}$

6. $\overline{x} = \frac{\int_{0}^{20} x 4e^{0.020732x} \, dx}{\int_{0}^{20} 4e^{0.020732x} \, dx} \approx 10.6891\text{,}$

Solution
1. $M = \int_{0}^{20} 4 + 0.1x \, dx = 100$ g.

2. The center of mass will be greater than 10 since the mass is more concentrated as $x$ increases.

3. Using the formula we developed, $\overline{x} = \frac{\int_{0}^{20} x \rho(x) \, dx}{\int_{0}^{20} \rho(x) \, dx}\text{,}$ so

\begin{equation*} \overline{x} = \frac{\int_{0}^{20} x (4 + 0.1x)) \, dx}{\int_{0}^{20} 4 + 0.1x \, dx} = \frac{\frac{3200}{3}}{100} = \frac{32}{3}\text{,} \end{equation*}

which is indeed just greater than 10.

4. Recall that average density is given by mass divided by length, so since the bar is 20 cm long and its mass is 100 g, the average density is 100/20 = 5 g/cm. This can also be computed using the formula for average value of a function:

\begin{equation*} \rho_{\text{AVG} [0,20]} = \frac{1}{20-0} \int_{0}^{20} \rho(x) \,dx = \frac{100}{20}\text{.} \end{equation*}
5. Since $p(x) = 4e^{0.020732x}$ lies below $\rho(x) = 4 + 0.1x$ up until where they intersect at about $17.53\text{,}$ after which $p(x) \gt \rho(x)\text{,}$ we see that $p(x)$ distributes mass more to the right. It also appears as if the two curves may bound the same total area, and thus the generate the same total mass. If that is the case, we expect $\overline{x}$ for $p(x)$ to lie slightly to the right of the center of mass for $\rho(x)\text{.}$

6. We see that

\begin{equation*} \overline{x} = \frac{\int_{0}^{20} x 4e^{0.020732x} \, dx}{\int_{0}^{20} 4e^{0.020732x} \, dx} \approx 10.6891\text{,} \end{equation*}

which is indeed slightly right of $\overline{x}$ for $\rho(x)\text{.}$

### SubsectionSummary

• For an object of constant density $D\text{,}$ with volume $V$ and mass $m\text{,}$ we have $m = D \cdot V\text{.}$

• If an object with constant cross-sectional area (such as a thin bar) has its density distributed along an axis according to the function $\rho(x)\text{,}$ then we can find the mass of the object between $x = a$ and $x = b$ by

\begin{equation*} m = \int_a^b \rho(x) \, dx\text{.} \end{equation*}
• If the density of an object (such as a disk) depends on the distance from one central point according to the function $\rho(r)\text{,}$ then we can find the mass of the object by

\begin{equation*} m = \int_0^R 2\pi r \rho(x) \, dx\text{.} \end{equation*}
• For a system of point-masses distributed along an axis, say $m_1, \ldots, m_n$ at locations $x_1, \ldots, x_n\text{,}$ the center of mass, $\overline{x}\text{,}$ is given by the weighted average

\begin{equation*} \overline{x} = \frac{\sum_{i=1}^n x_i m_i}{\sum_{i=1}^n m_i}\text{.} \end{equation*}

If instead we have mass continuously distributed along an axis, such as by a density function $\rho(x)$ for a thin bar of constant cross-sectional area, the center of mass of the portion of the bar between $x = a$ and $x = b$ is given by

\begin{equation*} \overline{x} = \frac{\int_a^b x \rho(x) \, dx}{\int_a^b \rho(x) \, dx}\text{.} \end{equation*}

In each situation, $\overline{x}$ represents the balancing point of the system of masses or of the portion of the bar.

### SubsectionExercises

Let a thin rod of length $a$ have density distribution function $\rho(x) = 10e^{-0.1x}\text{,}$ where $x$ is measured in cm and $\rho$ in grams per centimeter.

1. If the mass of the rod is 30 g, what is the value of $a\text{?}$

2. For the 30g rod, will the center of mass lie at its midpoint, to the left of the midpoint, or to the right of the midpoint? Why?

3. For the 30g rod, find the center of mass, and compare your prediction in (b).

4. At what value of $x$ should the 30g rod be cut in order to form two pieces of equal mass?

Consider two thin bars of constant cross-sectional area, each of length 10 cm, with respective mass density functions $\rho(x) = \frac{1}{1+x^2}$ and $p(x) = e^{-0.1x}\text{.}$

1. Find the mass of each bar.

2. Find the center of mass of each bar.

3. Now consider a new 10 cm bar whose mass density function is $f(x) = \rho(x) + p(x)\text{.}$

1. Explain how you can easily find the mass of this new bar with little to no additional work.

2. Similarly, compute $\int_0^{10} xf(x) \, dx$ as simply as possible, in light of earlier computations.

3. True or false: the center of mass of this new bar is the average of the centers of mass of the two earlier bars. Write at least one sentence to say why your conclusion makes sense.

Consider the curve given by $y = f(x) = 2xe^{-1.25x} + (30-x) e^{-0.25(30-x)}\text{.}$

1. Plot this curve in the window $x = 0 \ldots 30\text{,}$ $y = 0 \ldots 3$ (with constrained scaling so the units on the $x$ and $y$ axis are equal), and use it to generate a solid of revolution about the $x$-axis. Explain why this curve could generate a reasonable model of a baseball bat.

2. Let $x$ and $y$ be measured in inches. Find the total volume of the baseball bat generated by revolving the given curve about the $x$-axis. Include units on your answer.

3. Suppose that the baseball bat has constant weight density, and that the weight density is $0.6$ ounces per cubic inch. Find the total weight of the bat whose volume you found in (b).

4. Because the baseball bat does not have constant cross-sectional area, we see that the amount of weight concentrated at a location $x$ along the bat is determined by the volume of a slice at location $x\text{.}$ Explain why we can think about the function $\rho(x) = 0.6 \pi f(x)^2$ (where $f$ is the function given at the start of the problem) as being the weight density function for how the weight of the baseball bat is distributed from $x = 0$ to $x = 30\text{.}$

5. Compute the center of mass of the baseball bat.

The Math Olympics gives out prize medals, which are disks with varying density.

1. The first place medal has a radius of 2 inches, and the density of the disk is given by $\rho(r)=9-2r$ where $r$ is the distance from the center of the disk. Find the mass of the first place medal.

2. To make the second place medal, the Math Olympics takes a first-place medal and removes material from the outside until the radius is 1.8 inches. Find the mass of the second place medal.