Using Newton’s Law of Cooling

Exponential decay can also be applied to temperature. When a hot object is left in surrounding air that is at a lower temperature, the object’s temperature will decrease exponentially, leveling off as it approaches the surrounding air temperature. On a graph of the temperature function, the leveling off will correspond to a horizontal asymptote at the temperature of the surrounding air. Unless the room temperature is zero, this will correspond to a vertical shift of the generic exponential decay function. This translation leads to Newton’s Law of Cooling, the scientific formula for temperature as a function of time as an object’s temperature is equalized with the ambient temperature.

The formula is derived as follows:

[latex]\begin{array}{l}T\left(t\right)=A{b}^{ct}+{T}_{s}\hfill & \hfill \\ T\left(t\right)=A{e}^{\mathrm{ln}\left({b}^{ct}\right)}+{T}_{s}\hfill & \text{Properties of logarithms}.\hfill \\ T\left(t\right)=A{e}^{ct\mathrm{ln}b}+{T}_{s}\hfill & \text{Properties of logarithms}.\hfill \\ T\left(t\right)=A{e}^{kt}+{T}_{s}\hfill & \text{Rename the constant }c \mathrm{ln} b,\text{ calling it }k.\hfill \end{array}[/latex]

Exponential growth cannot continue forever. Exponential models, while they may be useful in the short term, tend to fall apart the longer they continue. Consider an aspiring writer who writes a single line on day one and plans to double the number of lines she writes each day for a month. By the end of the month, she must write over 17 billion lines or one-half-billion pages. It is impractical, if not impossible, for anyone to write that much in such a short period of time. Eventually an exponential model must begin to approach some limiting value and then the growth is forced to slow. For this reason, it is often better to use a model with an upper bound instead of an exponential growth model although the exponential growth model is still useful over a short term before approaching the limiting value.

The logistic growth model is approximately exponential at first, but it has a reduced rate of growth as the output approaches the model’s upper bound called the carrying capacity. For constants a, b, and c, the logistic growth of a population over time x is represented by the model

[latex]f\left(x\right)=\dfrac{c}{1+a{e}^{-bx}}[/latex]

The graph below shows how the growth rate changes over time. The graph increases from left to right, but the growth rate only increases until it reaches its point of maximum growth rate at which the rate of increase decreases.

Example: Using the Logistic-Growth Model

An influenza epidemic spreads through a population rapidly at a rate that depends on two factors. The more people who have the flu, the more rapidly it spreads, and also the more uninfected people there are, the more rapidly it spreads. These two factors make the logistic model good for studying the spread of communicable diseases. And, clearly, there is a maximum value for the number of people infected: the entire population.

For example, at time t = 0 there is one person in a community of 1,000 people who has the flu. So, in that community, at most 1,000 people can have the flu. Researchers find that for this particular strain of the flu, the logistic growth constant is b = 0.6030. Estimate the number of people in this community who will have had this flu after ten days. Predict how many people in this community will have had this flu after a long period of time has passed.

Show Solution

We substitute the given data into the logistic growth model

[latex]f\left(x\right)=\frac{c}{1+a{e}^{-bx}}[/latex]

Because at most 1,000 people, the entire population of the community, can get the flu, we know the limiting value is c = 1000. To find a, we use the formula that the number of cases at time t = 0 is [latex]\frac{c}{1+a}=1[/latex], from which it follows that a = 999. This model predicts that, after ten days, the number of people who have had the flu is [latex]f\left(x\right)=\frac{1000}{1+999{e}^{-0.6030x}}\approx 293.8[/latex]. Because the actual number must be a whole number (a person has either had the flu or not) we round to 294. In the long term, the number of people who will contract the flu is the limiting value, c = 1000.

Analysis of the Solution

Remember that because we are dealing with a virus, we cannot predict with certainty the number of people infected. The model only approximates the number of people infected and will not give us exact or actual values.

The graph below gives a good picture of how this model fits the data.

The graph of [latex]f\left(x\right)=\frac{1000}{1+999{e}^{-0.6030x}}[/latex].

Choosing an Appropriate Model

Now that we have discussed various mathematical models, we need to learn how to choose the appropriate model for the raw data we have. Many factors influence the choice of a mathematical model among which are experience, scientific laws, and patterns in the data itself. Not all data can be described by elementary functions. Sometimes a function is chosen that approximates the data over a given interval. For instance, suppose data were gathered on the number of homes bought in the United States from the years 1960 to 2013. After plotting these data in a scatter plot, we notice that the shape of the data from the years 2000 to 2013 follow a logarithmic curve. We could restrict the interval from 2000 to 2010, apply regression analysis using a logarithmic model, and use it to predict the number of home buyers for the year 2015.

Three kinds of functions that are often useful in mathematical models are linear functions, exponential functions, and logarithmic functions. If the data lies on a straight line or seems to lie approximately along a straight line, a linear model may be best. If the data is non-linear, we often consider an exponential or logarithmic model although other models, such as quadratic models, may also be considered.

In choosing between an exponential model and a logarithmic model, we look at the way the data curves. This is called the concavity. If we draw a line between two data points, and all (or most) of the data between those two points lies above that line, we say the curve is concave down. We can think of it as a bowl that bends downward and therefore cannot hold water. If all (or most) of the data between those two points lies below the line, we say the curve is concave up. In this case, we can think of a bowl that bends upward and can therefore hold water. An exponential curve, whether rising or falling, whether representing growth or decay, is always concave up away from its horizontal asymptote. A logarithmic curve is always concave down away from its vertical asymptote. In the case of positive data, which is the most common case, an exponential curve is always concave up and a logarithmic curve always concave down.

A logistic curve changes concavity. It starts out concave up and then changes to concave down beyond a certain point, called a point of inflection.

After using the graph to help us choose a type of function to use as a model, we substitute points, and solve to find the parameters. We reduce round-off error by choosing points as far apart as possible.

Example: Choosing a Mathematical Model

Does a linear, exponential, logarithmic, or logistic model best fit the values listed below? Find the model, and use a graph to check your choice.

x	1	2	3	4	5	6	7	8	9
y	0	1.386	2.197	2.773	3.219	3.584	3.892	4.159	4.394

Show Solution

First, plot the data on a graph as in the graph below. For the purpose of graphing, round the data to two significant digits.

Clearly, the points do not lie on a straight line, so we reject a linear model. If we draw a line between any two of the points, most or all of the points between those two points lie above the line, so the graph is concave down, suggesting a logarithmic model. We can try [latex]y=a\mathrm{ln}\left(bx\right)[/latex]. Plugging in the first point, [latex]\left(\text{1,0}\right)[/latex], gives [latex]0=a\mathrm{ln}b[/latex]. We reject the case that a = 0 (if it were, all outputs would be 0), so we know

[latex]\mathrm{ln}\left(b\right)=0[/latex]. Thus b = 1 and [latex]y=a\mathrm{ln}\left(\text{x}\right)[/latex]. Next we can use the point [latex]\left(\text{9,4}\text{.394}\right)[/latex] to solve for a:
[latex]\begin{array}{l}y=a\mathrm{ln}\left(x\right)\hfill \\ 4.394=a\mathrm{ln}\left(9\right)\hfill \\ a=\frac{4.394}{\mathrm{ln}\left(9\right)}\hfill \end{array}[/latex]

Because [latex]a=\frac{4.394}{\mathrm{ln}\left(9\right)}\approx 2[/latex], an appropriate model for the data is [latex]y=2\mathrm{ln}\left(x\right)[/latex].

To check the accuracy of the model, we graph the function together with the given points.

The graph of [latex]y=2\mathrm{ln}x[/latex].

We can conclude that the model is a good fit to the data. Compare the figure above to the graph of [latex]y=\mathrm{ln}\left({x}^{2}\right)[/latex] shown below.

The graph of [latex]y=\mathrm{ln}\left({x}^{2}\right)[/latex]

The graphs appear to be identical when x > 0. A quick check confirms this conclusion: [latex]y=\mathrm{ln}\left({x}^{2}\right)=2\mathrm{ln}\left(x\right)[/latex] for x > 0. However, if x < 0, the graph of [latex]y=\mathrm{ln}\left({x}^{2}\right)[/latex] includes an “extra” branch as shown below. This occurs because while [latex]y=2\mathrm{ln}\left(x\right)[/latex] cannot have negative values in the domain (as such values would force the argument to be negative), the function [latex]y=\mathrm{ln}\left({x}^{2}\right)[/latex] can have negative domain values.

Expressing an Exponential Model in Base e

While powers and logarithms of any base can be used in modeling, the two most common bases are [latex]10[/latex] and [latex]e[/latex]. In science and mathematics, the base e is often preferred. We can use properties of exponents and properties of logarithms to change any base to base e.