Skip to main content

Section 12.7 Functions Defined by Their Rates

When we learned about definite integrals, we learned that the definite integral \(\displaystyle \int_a^b f(x) \, dx\) computes the total change in a quantity that depends on \(x\) when \(x\) changes from \(a\) to \(b\) and where \(f(x)\) represents the rate of change of that quantity with respect to \(x\text{.}\) We have worked from an intuitive idea of rate of change with respect to time using concepts like velocity being the rate of change of position or flow rates (as in gallons per minute) as being the rate of change of volume with respect to time. We are now preparing to learn more specifically what we mean by rate of change, namely introducing the concept of the derivative.

Subsection 12.7.1 Describing Function Behavior

One of the consequences of the properties of definite integrals is that we can describe certain behaviors of a quantity in terms of the properties of the rate of change. For sequences, we learned that the properties of a sequence are determined from the increments. That is, a sequence is

  • increasing when its increments are positive,

  • decreasing when its increments are negative,

  • concave up when its increments are increasing,

  • concave down when its increments are decreasing.

(See Theorem 13.1.7 and Definition 13.1.8.) We learned analogous properties of functions defined by an accumulation of changes defined by a rate of change. That is, for a quantity \(Q\) with a rate of change \(R\text{:}\)

  • \(Q\) is increasing when its rate of change is \(R\) is positive,

  • \(Q\) is decreasing when its rate of change \(R\) is negative.

(See Theorem 3.4.5.) We define concavity for functions analogously.

Definition 12.7.1. Concavity of Functions.

A quantity \(Q\) that is a function of an independent variable \(x\) with a corresponding rate of change \(R\) that is itself a function of \(x\) has concavity that is determined by whether the rate is increasing or decreasing. Suppose that \(I=(a,b)\) is an interval.

  • \(Q\) is concave up on \(I\) if \(R\) is increasing on \(I\text{.}\)

  • \(Q\) is concave down on \(I\) if \(R\) is decreasing on \(I\text{.}\)

Concavity is closely related to the concept of acceleration (by which we also include the idea of deceleration). A constant rate of change leads to a linear relation. On a graph, this is a straight line. Concavity refers to a rate of change that is itself changing. This is acceleration. In physics, acceleration is caused by a force, so we can think of concavity as the effect on an object in the presence of a force.

Suppose that a quantity has a positive rate of change (increasing) and is also concave up (an increasing rate of change). This would be like a car moving forward (positive rate) with a rocket pushing it forward (positive acceleration). The result would be that the car continues to go faster (increasing rate), covering ever increasing distances per unit time. A graph of the position would be rising (increasing) and bending up (concave up).

Figure 12.7.2. \(Q\) is an increasing and concave up function of \(x\text{.}\) The rate is positive and increasing.

Next, suppose that a quantity has a positive rate of change (increasing) but is concave down (a decreasing rate of change). This would be like a car moving forward (positive rate) but with a rocket in reverse (negative acceleration). The car would still be moving forward, but the rocket is slowing it down. A graph of position in this case would be rising (increasing) but bending down (concave down).

Figure 12.7.3. \(Q\) is an increasing and concave down function of \(x\text{.}\) The rate is positive but decreasing.

If the rocket continues to exert a negative force, there will be a moment when all of the forward momentum is gone and then the car begins to go backwards. Consequently, we learn that a quantity that is concave down can switch from increasing to decreasing. Graphically, this is exactly what a parabola that opens down does. However, if the rocket is gradually reduced, we might be able to slow the car down without ever changing direction.

Figure 12.7.4. \(Q\) is a concave down function of \(x\) that is increasing for \(x \lt a\) and decreasing for \(x \gt a\text{.}\)

Similar behaviors might be described for negative rates of change. This would correspond to a car that is already going backwards. Being concave up (increasing rate of change) corresponds to a positive acceleration (rocket force), which in this case is opposite the motion and would serve to slow the car down maybe to the point of reversing direction. Graphically, this corresponds to a dropping graph that is bending up (moving toward flat). Being concave down (decreasing rate of change) corresponds to negative acceleration (rocket force) which now is the same direction as the motion. This would cause the car to speed up (in the negative direction). Graphically, being concave down corresponds to a graph that is bending down and growing ever steeper.

Example 12.7.5.

Suppose \(V\) measures the volume of water (liters) in a container and that \(V\) is a function of time \(t\) (minutes) such that the rate of change (liters per minute) is also a function of time defined by

\begin{equation*} \frac{dV}{dt} = R(t) = 10-0.25t. \end{equation*}

Describe the behavior of \(V\) and sketch a representative graph.

Solution

The rate of change of the volume in the container, \(R(t)\) determines the behavior of the volume. Because \(R\) has a negative slope \(m=-0.25\text{,}\) the rate is decreasing. This tells us that the volume is a concave down function. Solving the inequalities \(R(t) \gt 0\) and \(R(t) \lt 0\) will allow us to see when the rate is positive or negative, which will imply when the volume is increasing or decreasing, respectively.

The inequalities are solved by solving the equation \(R(t)=0\) and then testing the inequalities in the resulting intervals.

\begin{gather*} 10-0.25t = 0 \\ -0.25t = -10\\ t = \frac{-10}{-0.25} = 40 \end{gather*}

Testing the sign of \(R(t)\) when \(t \lt 40\text{,}\) we find, for example \(R(20) = 10-0.25(20) = 5\text{,}\) that the rate is positive. Consequently, \(V\) is increasing when \(t \lt 40\text{.}\) On the other hand, testing the sign of \(R(t)\) when \(t \gt 40\text{,}\) such as \(R(60) = 0-0.25(60)=-5\text{,}\) we find that the rate is negative so that the volume is decreasing when \(t \gt 40\text{.}\)

The graph of \(R(t)\text{,}\) shown below, is consistent with these analyses. The graph is decreasing (corresponding to the negative slope), above the axis for \(t \lt 40\) and below the axis for \(t \gt 40\text{.}\)

Figure 12.7.6. A graph of the rate of change of volume of water as a function of time.

The graph of the volume therefore needs to be concave down, increasing for \(t \lt 40\) and decreasing for \(t \gt 40\text{.}\) We do not know the starting volume (it wasn't given), so we might measure the change of volume from the initial value. That is, \(V=0\) on the graph will mean the volume is the same as the initial volume, rather than meaning there is no water. Using our knowledge of (((Unresolved xref, reference "thm-elementary-definite-integrals"; check spelling or use "provisional" attribute)))definite integrals of elementary algebraic formulas and the (((Unresolved xref, reference "thm-definite-integral-linearity"; check spelling or use "provisional" attribute)))linearity of definite integrals, we can find an explicit formula for our volume:

\begin{align*} V(t) &= \int_0^t R(z) \, dz = \int_0^t 10-0.25z \, dz\\ &= \int_0^t 10 \, dz - 0.25 \int_0^t z \, dz = 10t - 0.25 \frac{t^2}{2}. \end{align*}

The graph of this functions is shown below.

Figure 12.7.7. A graph of the change in volume of water as a function of time.

Definite integrals allow us to compute the change in a quantity when we know the rate of change. We are now turning our attention to the reverse question. If we know how to describe the quantity itself as a function of time, how do we find its corresponding rate of change? That rate of change is called the derivative.

Subsection 12.7.2 Introduction to Differential Equations

The rate of change or derivative of a quantity \(Q\) with respect to an independent variable \(x\) is itself another variable in the state of the system. It is often the case for physical and biological systems that there is a relationship between \(Q\) and its derivative \(\displaystyle \frac{dQ}{dx}\text{.}\) When expressed as an equation, such a relationship is called a differential equation. A function for which the differential equation is satisfied is called a solution.

Example 12.7.8.

A population \(P\) is a function of time \(t\text{.}\) The derivative \(\frac{dP}{dt}\) is the rate of change of the population, which consists of the birth rate (total births per unit time) and death rate (total deaths per unit time). When each of these rates is proportional to the population size, we have Malthusian growth and can write the differential equation

\begin{equation*} \frac{dP}{dt} = b \cdot P - d \cdot P \end{equation*}

where \(b\) is the per capita birth rate and \(d\) is the per capita death rate. If we write \(r=b-d\) to define a per capita net growth rate \(r\text{,}\) then the differential equation simplifies to

\begin{equation*} \frac{dP}{dt} = r P. \end{equation*}
Example 12.7.9.

In physics, Newton's second law of motion is usually written as its equation \(F=ma\text{,}\) which says that the total force acting on a body is always equal to the mass of that body times the acceleration of the body. Because acceleration is the rate of change of velocity,

\begin{equation*} a=\frac{dv}{dt}, \end{equation*}

Newton's law is actually a differential equation if we can compute the force.

Consider a falling object with mass \(m\text{.}\) Gravity is a downward force acting on the body with gravitational force \(F_g = -mg\) where \(g\) is the constant acceleration due to gravity. In the absence of other forces (no air resistance), the differential equation from Newton's law would be written

\begin{equation*} F_g = ma \quad \Leftrightarrow \quad -mg = m\frac{dv}{dt} \quad \Leftrightarrow \quad \frac{dv}{dt} = -g. \end{equation*}

If there is air resistance, the resulting force is itself usually a function of the velocity of the object passing through the air, \(F_a = f(v)\text{.}\) Experimentally, it has been found that many objects follow a square law, that the air resistance is proportional to the square of the velocity. Because air resistance is always in opposition to motion, \(F_a\) must have the opposite sign as \(v\text{.}\) So we have

\begin{equation*} F_a = - \gamma v |v| = \begin{cases} -\gamma v^2, & v \ge 0 \\ +\gamma v^2, & v \lt 0. \end{cases} \end{equation*}

Consequently, the differential equation of a falling object with air resistance is

\begin{equation*} F=ma \quad \Leftrightarrow \quad F_a+F_g=m \frac{dv}{dt} \quad \Leftrightarrow \quad -\gamma v|v| - m g = m \frac{dv}{dt}. \end{equation*}

Knowing a differential equation can often allow us to understand much of the behavior of the system of interest simply by determining when the rate is predicted to be positive and negative. Reasoning through concavity can be a little more difficult. We will focus on the case where the differential equation only involves the quantity \(y\) and its derivative \(\frac{dy}{dx}\) but does not involve the independent variable \(x\text{.}\) Such an equation is called an autonomous differential equation.

Values where an autonomous rate function \(f(y)=0\) are called equilibrium solutions because the value of \(y\) will never change if it has that value. That is, \(y\) will be a constant function that satisfies the differential equation.

We often summarize the behavior of an autonomous differential equation using a phase line. A phase line is a number line representing the dependent quantity \(y\text{.}\) We mark the equilibrium solutions, where \(f(y)=0\text{,}\) on the line. Between these points, we draw arrows representing whether \(y\) is increasing or decreasing between those points. If we also know where the rate \(f(y)\) reaches extreme values, we mark the location of those extreme rate points on the number line as well to represent inflection points.

The behavior of the dependent quantity then depends on where its initial value starts. If the initial value is at an equilibrium, then the dependent variable has the same constant value for all values of the independent variable. Otherwise, the function will be either increasing or decreasing, moving toward or away from equilibrium solutions.

Example 12.7.11.

Consider an autonomous differential equation

\begin{equation*} \frac{dy}{dx} = f(y) \end{equation*}

where the rate function \(f(y)\) is illustrated below. Create a phase line and use it to sketch expected shapes for solutions representing initial values in the different regions identified in the phase line.

Solution

The graph of \(\frac{dy}{dx}=f(y)\) has zeros at \(y=0\) and \(y=\pm 3\text{.}\) These correspond to equilibrium solutions of the differential equation. That is, \(y(x)=3\) (constant function) is one of the possible solutions. If the initial value is \(y(0)=3\text{,}\) then \(y(x)=3\) for all other values of \(x\) as well. Similarly for \(y(0)=0\) or \(y(0)=-3\text{.}\) These will be the points shown on the phase line.

The three equilibrium points divide the phase line into four intervals: \((-\infty, -3)\text{,}\) \((-3,0)\text{,}\) \((0,3)\) and \((3,\infty)\text{.}\) The graph of \(\frac{dy}{dx}=f(y)\) allows us to look at the sign of \(\frac{dy}{dx}\) for each of those regions. We see that \(\frac{dy}{dx} \gt 0\) for \(y \in (-\infty,-3)\text{,}\) so \(y(x)\) will be an increasing function if the initial value starts in this interval. We summarize the behavior in each interval with the following table.

Interval \((-\infty,-3)\) \((-3,0)\) \((0,3)\) \((3,\infty)\)
Sign of \(\frac{dy}{dx}=f(y)\) \(+\) (pos) \(-\) (neg) \(+\) (pos) \(-\) (neg)
Behavior of \(y\) increasing decreasing increasing decreasing

The phase line is a graphical summary of the table. Marked points on the line represent the equilibrium solutions. Arrows above the line represent the behavior as direction of motion. We will also include the locations of extreme rates, namely \(y=\pm \sqrt{3}\text{,}\) on the phase line as tick marks (not points) to indicate the locations of inflection points.

When we translate the information in the phase line to a sketch of the graph of solutions, we think of the phase line as the \(y\)-axis. The information about whether the function is increasing or decreasing is translated into the graph. Equilibrium solutions are also horizontal asymptotes for the solutions in the adjacent intervals. So we need to level off whenever the solution gets close to an equilibrium.