Skip to main content
Logo image

Active Calculus

Matthew Boelkins

Section 6.3 Density, Mass, and Center of Mass

Studying the units on 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.

Preview Activity 6.3.1.

In each of the following scenarios, we consider the distribution of a quantity along an axis.
Part a
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 number of miles east of a major interchange, and consider the definite integral \(\int_0^2 (200 + 100 e^{-0.1x}) \, dx\text{.}\)
(i) What are the units on the product \(c(x) \cdot \triangle x\text{?}\)
Units:
  • cars
  • miles
  • cars/mile
  • miles/car
  • cars*cars
  • miles*miles
(ii) What are the units on 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) \triangle x? \end{equation*}
Units:
  • cars
  • miles
  • cars/mile
  • miles/car
  • cars*cars
  • miles*miles
(iii) Evaluate the definite integral \(\int_0^2 c(x) \, dx = \int_0^2 (200 + 100 e^{-0.1x}) \, dx =\)
What does the answer to the definite integral \(\int_0^2 c(x) \, dx\) mean?
  • The number of cars in the first two miles.
  • The number of cars at mile 2.
  • The density of cars in the first two miles.
  • The density of cars at mile 2.
  • The area the cars occupy in the first two miles.
  • The area the cars occupy at mile 2.
Part b
On a 6 foot long shelf filled with books, the function \(B\) models the distribution of the weight of the books, measured 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 \(\displaystyle B(x) = 0.5 + \frac{1}{(x+1)^2}\text{.}\)
(i) What are the units on the product \(B(x) \cdot \triangle x\text{?}\)
Units:
  • pounds
  • inches
  • pound/inch
  • inch/pound
  • pound*pound
  • inch*inch
(ii) What are the units on 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) \triangle x? \end{equation*}
Units:
  • pounds
  • inches
  • pound/inch
  • inch/pound
  • pound*pound
  • inch*inch
(iii) Evaluate the definite integral \(\displaystyle \int_{0}^{72} B(x) \, dx = \int_0^{72} (0.5 + \frac{1}{(x+1)^2}) \, dx=\)
What does the answer to the definite integral \(\int_0^{72} B(x) \, dx\) mean?
  • The total weight of the books on the entire shelf.
  • The total weight of the books at the end of the shelf.
  • The distribution of the weight of the books on the entire shelf.
  • The distribution of the weight of the books at the end of the shelf.
  • The area the books occupy on the entire shelf.
  • The area the books occupy at the end of the shelf.

Subsection 6.3.1 Density

The mass of a quantity, typically measured in metric units such as grams or kilograms, is a measure of the amount of the quantity. In a corresponding way, 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.
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\) we know that the distance traveled is approximately \(v(t) \Delta t\text{,}\) 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\) 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 Figure 6.3.1.
Figure 6.3.1. At 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.
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{.}\)
Figure 6.3.2. A 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 Figure 6.3.2, 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*}
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{,} \end{equation*}
therefore measure the mass of the bar between \(0\) and \(b\text{.}\) (The Riemann sum is an approximation, while the integral will be the exact mass.)
For objects whose cross-sectional area is constant and whose mass is distributed relative to horizontal location, \(x\text{,}\) it makes sense to think of the density function \(\rho(x)\) with units “mass per unit length,” such as g/cm. Thus, when we compute \(\rho(x) \cdot \Delta x\) on a small slice \(\Delta x\text{,}\) the resulting units are g/cm \(\cdot\) cm = g, which thus measures the mass of the slice. The general principle follows.
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*}

Activity 6.3.2.

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 on \(\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 lying horizontally 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 \(z\) at which to cut the bar so that the two pieces will each have identical mass.

Subsection 6.3.2 Weighted 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 Table 6.3.3.
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
Table 6.3.3. A college student’s semester grades.
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 activity, we explore further how weighted averages can be used to find the balancing point of a physical system.

Activity 6.3.3.

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.

Subsection 6.3.3 Center of Mass

In Activity 6.3.3, we saw that the balancing point of a system of point-masses
 1 
In the activity, 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 activity, 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{.}\tag{6.3.1} \end{equation}
The greater the number of slices, the more accurate our estimate of the balancing point will be. The sums in Equation (6.3.1) 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.

Activity 6.3.4.

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

Subsection 6.3.4 Summary

  • For an object of constant density \(D\text{,}\) with volume \(V\) and mass \(m\text{,}\) we know that \(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*}
  • 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.

Exercises 6.3.5 Exercises

1. Center of mass for a linear density function.

A rod has length 4 meters. At a distance \(x\) meters from its left end, the density of the rod is given by \(\delta(x)=4 + 5 x\) g/m.
(a) Complete the Riemann sum for the total mass of the rod (use \(Dx\) in place of \(\Delta x\)):
mass = \(\Sigma\)
(b) Convert the Riemann sum to an integral and find the exact mass.
mass =
(include help (units)
 2 
/pg_files/helpFiles/Units.html
)

2. Center of mass for a nonlinear density function.

A rod with uniform density (mass/unit length) \(\delta(x) = 7+\sin\mathopen{}\left(x\right)\) lies on the \(x\)-axis between \(x= 0\) and \(x= \pi.\) Find the mass and center of mass of the rod.
mass =
center of mass =

3. Interpreting the density of cars on a road.

Suppose that the density of cars (in cars per mile) down a 20-mile stretch of the Pennsylvania Turnpike is approximated by \(\delta(x) = 350 \left(2+\sin\left(4\sqrt{x + 0.2}\right)\right),\) at a distance \(x\) miles from the Breezewood toll plaza. Sketch a graph of this function for \(0\leq x\leq 20\text{.}\)
(a) Complete the Riemann sum that approximates the total number of cars on this 20-mile stretch (use \(Dx\) instead of \(\Delta x\)):
Number = \(\Sigma\)
(b) Find the total number of cars on the 20-mile stretch.
Number =

4. Center of mass in a point-mass system.

A point mass of 5 grams located 7 centimeters to the left of the origin and a point mass of 7 grams located 9 centimeters to the right of the origin are connected by a thin, light rod. Find the center of mass of the system.
Center of Mass =
  • to the left of the origin
  • to the right of the origin
  • (at the origin)
(include help (units)
 3 
/pg_files/helpFiles/Units.html
in your center of mass)

5.

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?

6.

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.

7.

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.
You have attempted of activities on this page.
PrevTopNext
Feedback