Exponential growth and decay (Part 16): Logistic growth model

In this series of posts, I provide a deeper look at common applications of exponential functions that arise in an Algebra II or Precalculus class. In the previous posts in this series, I considered financial applications, radioactive decay, and Newton’s Law of Cooling.

Today, I discuss the logistic growth model, which describes how an infection (like a disease, a rumor, or advertise) spreads in a population. In yesterday’s post, I described an in-class demonstration that engages students while also making the following formula believable:

A(t) = \displaystyle \frac{Ly_0}{y_0+ (L-y_0)e^{-rt}}.

I’d like to discuss some observations about this somewhat complicated function that will make producing its graph easier. The first two observations are within reach of Precalculus students.

1. Let’s figure out the y-intercept:

A(0) = \displaystyle \frac{Ly_0}{y_0+ (L-y_0)e^{-r \cdot 0}} = \displaystyle \frac{Ly_0}{y_0+ L-y_0} = y_0.

In other words, the number y_0 represents the initial number of people who have the infection.

2. Let’s figure out the limiting value as t gets large:

\displaystyle \lim_{t \to \infty} A(t) = \displaystyle \frac{Ly_0}{y_0+ (L-y_0) \cdot 0} = \displaystyle \frac{Ly_0}{y_0} = L.

As expected, all L people will get the infection eventually. (Of course, Precalculus students won’t be familiar with the $\displaystyle \lim$ notation, but they should understand that e^{-rt} decays to zero as t gets large.

3. Let’s now figure out the point of inflection. Ordinarily, points of inflection are found by setting the second derivative equal to zero. Though this can be done for the function A(t) above, it would be a somewhat daunting exercise!

The good news is that the points of inflection can be found quite simply using the governing differential equation, which is

A' = r A [ L - A] = r L A - r A^2

Let’s take the derivative of both sides, remembering that r and L are constants:

A'' = r L A' - 2 r A A'

A'' = A' (r L - 2 r A)

So the second derivative is equal to zero when either A' = 0 or else r L - 2 r A = 0. The first case corresponds to the trivial cases A(t) \equiv 0 and A(t) \equiv L; these constants are called the equilibrium solutions. The second case is the more interesting one:

r L - 2 r A = 0

r L = 2 r A

\displaystyle \frac{L}{2} = A

This suggests that, as the infection spreads throughout a population, the curve changes concavity at the time that half of the population becomes infected. In other words, the infection spreads fastest throughout the population at the time when half of the population has been infected.

The time at which the point of inflection occurs can be found by setting A(t) = \displaystyle \frac{L}{2} and solving for t:

\displaystyle \frac{L}{2} = \displaystyle \frac{Ly_0}{y_0+ (L-y_0)e^{-rt}}.

\displaystyle \frac{1}{2} = \displaystyle \frac{y_0}{y_0+ (L-y_0)e^{-rt}}.

y_0 + (L-y_0) e^{-rt} = 2y_0

(L-y_0) e^{-rt} = y_0

e^{-rt} = \displaystyle \frac{y_0}{L-y_0}

-rt = \displaystyle \ln \left( \frac{y_0}{L-y_0} \right)

t = \displaystyle - \frac{1}{r} \ln \left( \frac{y_0}{L-y_0} \right)

This technique for finding the points of inflection directly from the differential equation is possible whenever the differential equation is autonomous, which loosely means that the independent variable does not appear on the right-hand side.

 

Leave a comment

1 Comment

  1. Exponential growth and decay: Index | Mean Green Math

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

This site uses Akismet to reduce spam. Learn how your comment data is processed.

%d bloggers like this: