# 6.3: Density, Mass, and Center of Mass

- Page ID
- 4329

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Learning Objectives

In this section, we strive to understand the ideas generated by the following important questions:

- How are mass, density, and volume related?
- How is the mass of an object with varying density computed?
- What is the center of mass of an object, and how are definite integrals used to compute it?

We have seen in several different circumstances how studying the units on the integrand and variable of integration enables us to better 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, while \( t\) measures time in seconds, then both the definite integral and its Riemann sum approximation,

\( \int^b_a v(t) \,dt ≈ \sum^n_{i=1} v(t_i)\Delta t\),

have their overall units given by the product of the units of \(v(t)\) and \(t\):

(feet/sec)·(sec) = feet.

Thus,

\(\int^b_a v(t) \,dt\)

measures the total change in position (in feet) of the moving object. This type of unit analysis will be particularly helpful to us in what follows. To begin, in the following preview activity we consider two different definite integrals where the integrand is a function that measures how a particular quantity is distributed over a region and think about how the units on the integrand and the variable of integration indicate the meaning of the integral.

Preview Activity \(\PageIndex{1}\):

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

- Suppose that the function \(c(x) = 200 + 100e^{−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^2_0 (200 + 100e^{−0.1x} )\, dx\).
- What are the units on the product \( c(x) \cdot \Delta x\)?
- What are the units on the definite integral and its Riemann sum approximation given by

\( \int_0^2 c(x)dx \approx\sum_{i=1}^{n}c(x_{i})\Delta x\) ? - Evaluate the definite integral \( \int^2_0 c(x)\ dx = \int^2_0 (200 + 100e^{−0.1x})\ dx\) and write one sentence to explain the meaning of the value you find.

- 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 \( B(x) = 0.5 + \dfrac{1}{(x+1)^2}\).
- What are the units on the product \( B(x) \cdot \Delta x\)?
- What are the units on the definite integral and its Riemann sum approximation given by

\( \int^{36}_{12} B(x)\ dx ≈ \sum^n_{i=1}B(x_i) \Delta x\)? - Evaluate the definite integral \( \int^{72}_0 B(x) dx = \int^{72}_0 \left(0.5 + \dfrac{1}{(x+1)^2}\right) dx\) and write one sentence to explain the meaning of the value you find.

## Density

The *mass* of a quantity, typically measured in metric units such as grams or kilograms, is a measure of the amount of a 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 \( \text{m}^3\), then the density of the brick is

\[ \dfrac{3\text{kg}}{0.002 \text{m}^3} = 1500 \dfrac{\text{kg}} {\text{m}^3}.\]

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

For an object of constant density \( d\), with mass \( m\) and volume \( V\),

\[ d = \dfrac{m}{V}\]

or

\[ m = d \cdot V.\]

But what happens when the density is not constant?

If we consider the formula \( m = d \cdot V\), it is reminiscent of two other equations that we have used frequently in recent work: for a body moving in a fixed direction, distance = rate · time, and, for a rectangle, its area is given by \( A\ =\ l \cdot w\). 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. The main idea in each situation is that by working with small slices of the quantity that is varying, we can use a definite integral to add up the values of small pieces on which the quantity of interest (such as the velocity of a moving object) are approximately constant.

For example, in the setting where 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\), since \( v(t)\) is almost constant on a small interval, and for a constant rate, distance = rate · time. 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\): \( f(x)\Delta x\). Both of these principles are represented visually in Figure \(\PageIndex{1}\).

**Figure \(\PageIndex{1}\):** At left, estimating a small amount of distance traveled, \( v(t)\Delta t\), and at right, a small amount of area under the curve, \( f(x)\Delta x\).

In a similar way, if we consider the setting where the density of some quantity is not constant, the definite integral enables us to still compute the overall mass of the quantity. Throughout, we will focus on problems where the density varies in only one dimension, say along a single axis, and think about how mass is distributed relative to location along the axis. Let’s consider a thin bar of length b that is situated so its left end is at the origin, where x = 0, and assume that the bar has constant cross-sectional area of 1 \( \text{cm}^2\). We let the function \( \rho(x)\) represent the mass density function of the bar, measured in grams per cubic centimeter. That is, given a location \( x\), \( \rho(x)\) tells us approximately how much mass will be found in a one-centimeter wide slice of the bar at \( x\).

**Figure \(\PageIndex{2}\):** A thin bar of constant cross-sectional area 1 \( \text{cm}^2\) with density function \( \rho(x) \dfrac{\text{g}}{\text{cm}^{3}}\).

If we now consider a thin slice of the bar of width \( \Delta x\), as pictured in Figure \(\PageIndex{2}\), the volume of such a slice is the cross-sectional area times \( \Delta x\). Since the cross-sections each have constant area 1 \( \text{cm}^2\), it follows that the volume of the slice is \( 1\Delta x\ \text{cm}^2\). Moreover, since mass is the product of density and volume (when density is constant), we see that the mass of this given slice is approximately

\( \text{mass}\ _{\text{slice}}\approx\rho(x)\dfrac{\text{g}}{\text{cm}^{3}}\cdot1\Delta x\ \text{cm}^{3}=\rho(x)\cdot\Delta x \cdot \text{g}\)

Hence, for the corresponding Riemann sum (and thus for the integral that it approximates),

\[ \sum^n_{i=1}\rho(x_i)\Delta x ≈ \int^b_0 \rho(x) dx,\]

we see that these quantities measure the mass of the bar between 0 and \( b\). (The Riemann sum is an approximation, while the integral will be the exact mass.)

At this point, we note that we will be focused primarily on situations where mass is distributed relative to horizontal location, \( x\), for objects whose cross-sectional area is constant. In that setting, 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\), the resulting units are \( \dfrac{\text{g}}{\text{cm}}\cdot\text{cm}=\text{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\) of the object between \( x = a\) and \( x = b\) is given by

\[M = \int^b_a ρ(x) dx.\]

Exercise \(\PageIndex{1}\):

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

- Suppose that a thin rod with constant cross-sectional area of 1 \( \text{cm}^{2}\) has its mass distributed according to the density function \( \rho(x) = 2e^{−0.2x}\), 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.
- 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.
- Write and evaluate a definite integral whose value is the volume of the cone.
- Next, suppose that the cone has uniform density of 800 kg/m3 . What is the mass of the solid cone?
- 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 + \dfrac{200}{1+x^2}\), measured in kg/\( \text{m}^{3}\). 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.

- Let a thin rod of constant cross-sectional area 1 \( \text{cm}^{2}\) and length 12 cm have its mass be distributed according to the density function \( \rho(x) = \dfrac{1}{25} (x − 15)^2\), 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.

## Weighted Averages

class grade grade points credits chemistry B+ 3.3

*Table 6.1: A college student’s semester grades.*

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,\ a_2, . . .,\ a_n\), we know that their average is given by

\[ \dfrac{a_1,\ a_2, . . .,\ a_n}{n},\]

and for a quantity being measured by a function \( f\) on an interval \( [a, b]\), the average value of the quantity on \( [a, b]\) is

\p \dfrac{1}{b − a} \int^b_a f (x) dx.\[

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

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

\[ \dfrac{3.3 + 3.7 + 2.7 + 2.7}{4} = 3.1.\]

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

\[ \dfrac{3.3 \cdot 5 + 3.7 \cdot 4 + 2.7 \cdot 3 + 2.7 \cdot 3}{5 + 4 + 3 + 3} =3.\overline{16}.\]

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 \(\PageIndex{1}\):

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.

- Suppose that a shelf is 6 feet long, with its left end situated at \( x = 0\). If one book of weight 1 lb is placed at \( x_1 = 0\), and another book of weight 1 lb is placed at \( x_2 = 6\), what is the location of \( \overline{x}\), the point at which the shelf would (theoretically) balance on a fulcrum?
- Now, say that we place four books on the shelf, each weighing 1 lb: at \( x_1 = 0\), at \( x_2 = 2\), at \( x_3 = 4\), and at \( x_4 = 6\). Find \( \overline{x}\), the balancing point of the shelf.
- 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\), \( x_2 = 2\), \( x_3 = 3\), and \( x_4 = 6\).
- 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, and \( 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}\), the balancing point of the shelf. How does the balancing point in this scenario compare to that found in (b)?
- 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\). How does \( \overline{x}\) change?
- 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?
- Experiment with a couple of different scenarios of your choosing where you move the location of one of the books to the left, or you decrease the weight of one of the books.
- 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.

## Center of Mass

In Activity \(\PageIndex{2}\), we saw that the balancing point of a system of point-masses^{1} (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.

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

\( \overline{x} = \dfrac{x_1m_1 + x_2m_2 + · · · x_n m_n}{m_1 + m_2 + · · · + m_n} \).

What if we instead 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 1 In the activity, we actually used weight rather than mass. Since weight is computed by the gravitational constant times mass, the computations for the balancing point result in the same location regardless of whether we use weight or mass, since the gravitational constant is present in both the numerator and denominator of the weighted average.

function \( \rho(x)\) tells us the density distribution along the bar, measured in g/cm. If we slice the bar into small sections, this enables us to think of the bar as holding a collection of adjacent point-masses. For a slice of thickness \( \Delta x\) at location \( x_i\), note that the mass of the slice, \( m_i\), satisfies \( m_i \approx \rho(x_i)\Delta x\).

Taking \( n\) slices of the bar, we can approximate its center of mass by

\( \overline{x} \approx \dfrac{ x_1 · ρ(x_1)\Delta x + x_2 · ρ(x_2)\Delta x + · · · + x_n · ρ(x_n)\Delta x}{ρ(x_1)\Delta x + ρ(x_2)\Delta x + · · · + ρ(x_n)\Delta x}\).

. Rewriting the sums in sigma notation, it follows that

\[ \overline{x} \approx \dfrac{\sum^n_{i=1}x_i · \rho(x_i)\Delta x}{\sum^n_{i=1}\rho(x_i)\Delta x}. \label{6.4} \]

Moreover, it is apparent that the greater the number of slices, the more accurate our estimate of the balancing point will be, and that the sums in Equation (\( \ref{6.4}\)) can be viewed as Riemann sums. Hence, in the limit as \( n\rightarrow\infty\), we find that the center of mass is given by the quotient of two integrals.

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

\( \overline{x} = \dfrac{\int^b_a x \rho(x) dx}{\int^b_a \rho(x) dx}\).

Note particularly 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\), along the bar.

Activity \(\PageIndex{1}\):

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

- Find the total mass, \( M\), of the bar.
- 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?
- Compute \( \overline{x}\), the exact center of mass of the bar.
- What is the average density of the bar?
- Now consider a different density function, given by \( p(x) = 4e^{0.020732x}\), also for a bar of length 20 cm whose left end is at \( x = 0\). 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?
- Compute the exact center of mass of the bar described in (e) whose density function is \( p(x) = 4e^{0.020732x}\). Check the result against the prediction you made in (e).

## Summary

In this section, we encountered the following important ideas:

- For an object of constant density \( D\), with volume \( V\) and mass \( m\), we know that \( m = D·V\).
- 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)\), then we can find the mass of the object between \( x = a\) and \( x = b\) by \[ m = \int^b_a \rho(x) dx.\]
- For a system of point-masses distributed along an axis, say \( m_1, . . ., m_n\) at locations \( x_1, . . ., x_n\), the center of mass, \( \overline{x}\), is given by the weighted average \[ \overline{x} = \dfrac{\sum^n_{i=1}x_im_i}{\sum^n_{i=1}m_i}.\]
- 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 \[ \overline{x} = \dfrac{\int^b_a x\rho(x)dx}{\int^b_a \rho(x)dx}.\]
- In each situation, \( \overline{x}\) represents the balancing point of the system of masses or of the portion of the bar.

## Contributors and Attributions

Matt Boelkins (Grand Valley State University), David Austin (Grand Valley State University), Steve Schlicker (Grand Valley State University)