Unit 3, lesson 46

Logistic Growth and Model Limits

How the logistic function models growth that starts out nearly exponential and levels off at a carrying capacity, what each parameter means, and where the model stops matching reality.

  • Theory part 91
  • Practice part 92
  • Models

A few fish released into an empty lake multiply quickly at first, because food and space are plentiful. As the lake fills, each new fish has less to eat, and growth slows until the population levels off. An exponential model gets the beginning right and the end badly wrong. The logistic model describes both stages with three parameters, each with a clear meaning.

Key terms

Logistic function

Definition A logistic function is a function of the form

P(t)=L1+Ae−ktP(t) = \frac{L}{1 + Ae^{-kt}}

with constants L>0L > 0, A>0A > 0 and k>0k > 0. Its domain is all real numbers and its range is the open interval (0,L)(0, L). It is increasing, and its graph has the horizontal asymptote y=0y = 0 as t→−∞t \to -\infty and y=Ly = L as t→∞t \to \infty. The graph has the shape of a stretched letter S. Books use different letters for these constants, and the calculus lesson on logistic growth writes KK in place of LL and rr in place of kk.

In plain words A logistic function starts near zero, climbs, and then flattens out under a ceiling it never touches.

Example A lake is stocked with 1010 fish of a new species, and the fish population tt years later is modeled by

F(t)=10001+99e−0.5t.F(t) = \frac{1000}{1 + 99e^{-0.5t}}.

Rounded to whole fish, the model gives these values.

tt in years004488121216162020
F(t)F(t)10106969355355803803968968996996

Positive parameters

Definition In P(t)=L1+Ae−ktP(t) = \frac{L}{1 + Ae^{-kt}} all three parameters must be positive. The starting value is

P(0)=L1+A,soA=LP(0)−1.\begin{aligned} P(0) &= \frac{L}{1 + A}, \\[0.8em] \text{so} \quad A &= \frac{L}{P(0)} - 1. \end{aligned}

With L>0L > 0, the condition A>0A > 0 holds exactly when 0<P(0)<L0 < P(0) < L, and k>0k > 0 makes PP increasing. Under these conditions the denominator 1+Ae−kt1 + Ae^{-kt} is greater than 11 for every tt, so PP is defined everywhere and stays between 00 and LL. If instead A<0A < 0 while k>0k > 0, the denominator is zero at t=ln⁡(−A)kt = \frac{\ln(-A)}{k} and the graph has a vertical asymptote there. If A>0A > 0 but k<0k < 0, the curve falls toward 00 as tt grows.

In plain words Positive parameters guarantee the S shape. The ceiling is positive, the start is between zero and the ceiling, and the curve goes up.

Example For the lake, L=1000L = 1000 and F(0)=10F(0) = 10, so

A=100010−1=99.A = \frac{1000}{10} - 1 = 99.

If a survey had instead counted 12001200 fish at the start, more than the lake can support, the formula would give

A=10001200−1=−16.A = \frac{1000}{1200} - 1 = -\frac{1}{6}.

With that negative AA and k=0.5k = 0.5, the curve falls from 12001200 toward 10001000 for t≥0t \ge 0 instead of rising in an S. A start above the capacity is outside the definition above, which describes growth from below.

Carrying capacity

Definition In a logistic model the parameter LL is the carrying capacity, the size the population approaches over time. Because 1+Ae−kt>11 + Ae^{-kt} > 1, every value satisfies 0<P(t)<L0 < P(t) < L, and because Ae−kt→0Ae^{-kt} \to 0 as t→∞t \to \infty,

P(t)→Last→∞.P(t) \to L \quad \text{as} \quad t \to \infty.

So the line y=Ly = L is a horizontal asymptote of the graph. It is approached but never reached.

In plain words The carrying capacity is the most the environment can hold in the long run. The model gets as close to it as you like but never gets there.

Example For the lake the carrying capacity is 10001000 fish. After 2020 years the model gives about 995.5995.5 fish, and after 3030 years

F(30)=10001+99e−15≈999.97,F(30) = \frac{1000}{1 + 99e^{-15}} \approx 999.97,

still below 10001000.

Early growth

Definition Multiplying P(t)P(t) by 1+Ae−kt1 + Ae^{-kt} and rearranging gives, for every tt, the exact identity

P(t)=LAekt(1−P(t)L).P(t) = \frac{L}{A}e^{kt}\left(1 - \frac{P(t)}{L}\right).

While P(t)P(t) is small compared with LL, the factor in parentheses is close to 11, so

P(t)≈LAekt.P(t) \approx \frac{L}{A}e^{kt}.

Early logistic growth is therefore nearly exponential with continuous growth rate kk. The exponential formula is too large by the factor 11−P(t)/L\frac{1}{1 - P(t)/L}, which grows as P(t)P(t) approaches LL.

In plain words At the start the ceiling is far away and barely matters, so the population grows like an exponential function. The approximation only fails once the population is a noticeable fraction of the capacity.

Example For the lake, LA=100099\frac{L}{A} = \frac{1000}{99}.

ttF(t)F(t)100099e0.5t\frac{1000}{99}e^{0.5t}
0010.0010.0010.1010.10
2226.7226.7227.4627.46
4469.4569.4574.6474.64
66168.66168.66202.88202.88
88355.46355.46551.50551.50

The two columns stay within about 88 percent for four years. By year 88 the exponential guess is more than 5050 percent too large.

Saturation

Definition Saturation is the late stage of logistic growth, when P(t)P(t) is close to LL and growth has nearly stopped. The remaining gap is

L−P(t)=LAe−kt1+Ae−kt,L - P(t) = \frac{LAe^{-kt}}{1 + Ae^{-kt}},

which is positive for every tt and approaches 00 as t→∞t \to \infty. Once Ae−ktAe^{-kt} is small the gap is approximately LAe−ktLAe^{-kt}, so it shrinks by about the factor e−ke^{-k} with each unit of time.

In plain words Near the ceiling the population creeps upward, and each year closes about the same fraction of the remaining gap.

Example Here is the gap for the lake every two years. Over two years the factor is e−0.5⋅2=e−1≈0.37e^{-0.5 \cdot 2} = e^{-1} \approx 0.37.

tt1414161618182020
1000−F(t)1000 - F(t)82.882.832.132.112.112.14.54.5

Each entry is close to 0.370.37 times the one before it.

Symmetry about the midpoint

Definition The midpoint of a logistic curve is the point (t0,L2)(t_0, \frac{L}{2}) with

t0=ln⁡Ak.t_0 = \frac{\ln A}{k}.

Since Ae−kt0=1Ae^{-kt_0} = 1, the function can be written as

P(t)=L1+e−k(t−t0),P(t) = \frac{L}{1 + e^{-k(t - t_0)}},

and P(t0)=L2P(t_0) = \frac{L}{2}. The graph is symmetric about the midpoint, which means

P(t0+s)+P(t0−s)=Lfor every real s.\begin{gathered} P(t_0 + s) + P(t_0 - s) = L \\[0.4em] \text{for every real } s. \end{gathered}

To see this, put u=e−ksu = e^{-ks}. Then

P(t0+s)=L1+uandP(t0−s)=L1+1/u=Lu1+u,\begin{gathered} P(t_0 + s) = \frac{L}{1 + u} \\[0.8em] \text{and} \\[0.4em] P(t_0 - s) = \frac{L}{1 + 1/u} = \frac{Lu}{1 + u}, \end{gathered}

which add to LL. The curve is also steepest at the midpoint, where it changes from bending upward to bending downward. Calculus proves that fact, and here it is taken as known.

In plain words Turn the graph half a turn about its midpoint and it lands on itself. The growth that happens just before the midpoint is mirrored by the growth just after it.

Example For the lake, t0=ln⁡990.5≈9.19t_0 = \frac{\ln 99}{0.5} \approx 9.19 years, when the population is 500500. Four years on either side,

F(t0+4)≈880.8andF(t0−4)≈119.2,\begin{aligned} F(t_0 + 4) &\approx 880.8 \\[0.4em] \text{and} \quad F(t_0 - 4) &\approx 119.2, \end{aligned}

which add to 10001000. The population gains about 380.8380.8 fish in the four years before the midpoint and the same number in the four years after.

Model interpretation without differential equations

Definition In calculus the logistic function comes from an equation about its rate of change, a differential equation. Without that equation, a logistic model can still be read entirely from its formula. The carrying capacity is LL, the starting value is L1+A\frac{L}{1 + A}, the early continuous growth rate is kk, and t0=ln⁡Akt_0 = \frac{\ln A}{k} is the time of half capacity and fastest growth. For any level yy with 0<y<L0 < y < L, solving y=L1+Ae−kty = \frac{L}{1 + Ae^{-kt}} with logarithms gives the unique time at which the model reaches yy,

t=1kln⁡(AyL−y).t = \frac{1}{k}\ln\left(\frac{Ay}{L - y}\right).

The levels 00 and LL are never reached. The model also rests on assumptions that real populations can break. The capacity stays fixed, growth is smooth with no seasons, harvests or disasters, and the rise is exactly symmetric about the midpoint.

In plain words The formula and logarithms answer every question the model can answer. Questions outside its assumptions, like what happens if the lake shrinks in a drought, it cannot answer at all.

Example The lake model reaches 900900 fish when

t=10.5ln⁡(99⋅9001000−900)=2ln⁡891≈13.58 years.\begin{aligned} t &= \frac{1}{0.5}\ln\left(\frac{99 \cdot 900}{1000 - 900}\right) \\[0.8em] &= 2\ln 891 \approx 13.58 \text{ years.} \end{aligned}

Asking when it reaches 10001000 fish has no answer, since the model only approaches that level. A real lake is also unlikely to settle at exactly 10001000. Its population usually rises and falls around the capacity from year to year, which this smooth curve cannot show.

Check yourself

  1. For the model P(t)=6001+5e−0.4tP(t) = \frac{600}{1 + 5e^{-0.4t}}, find the carrying capacity, the starting value and the time of the midpoint.
Show the answer

The carrying capacity is 600600. The starting value is P(0)=6001+5=100P(0) = \frac{600}{1 + 5} = 100. The midpoint time is

t0=ln⁡50.4≈4.02,t_0 = \frac{\ln 5}{0.4} \approx 4.02,

when the value is 300300.

  1. For the same model, when does P(t)P(t) reach 450450?
Show the answer

Setting P(t)=450P(t) = 450 gives

1+5e−0.4t=600450=43.1 + 5e^{-0.4t} = \frac{600}{450} = \frac{4}{3}.

So 5e−0.4t=135e^{-0.4t} = \frac{1}{3} and e−0.4t=115e^{-0.4t} = \frac{1}{15}. Taking natural logarithms gives −0.4t=−ln⁡15-0.4t = -\ln 15, so

t=ln⁡150.4≈6.77.t = \frac{\ln 15}{0.4} \approx 6.77.
  1. A logistic model has carrying capacity 800800 and midpoint time t0t_0. Three time units before the midpoint its value is 140140. What is its value three time units after the midpoint?
Show the answer

By symmetry about the midpoint, P(t0+3)+P(t0−3)=800P(t_0 + 3) + P(t_0 - 3) = 800, so P(t0+3)=800−140=660P(t_0 + 3) = 800 - 140 = 660.

How this lesson connects

Where this leads in Calculus