Math Core

Lesson 3.3 · Differential Equations

Logistic growth

Exponential growth dPdt=kP\dfrac{dP}{dt} = kP predicts that a population grows forever, faster and faster. Real populations run out of food, space or customers, and their growth levels off. The logistic differential equation models that leveling off, and AP Calculus BC expects you to read a lot from it without ever solving it.

The logistic differential equation

In the logistic model, the growth rate is proportional both to the population PP and to how far PP is from a ceiling LL.

Definition

Logistic differential equation

dPdt=kP(1−PL),k>0, L>0.\frac{dP}{dt} = kP\left(1 - \frac{P}{L}\right), \qquad k \gt 0,\ L \gt 0.

The constant kk is the relative growth rate when PP is small, and LL is the carrying capacity: the population level the environment can sustain.

The factor (1−PL)\left(1 - \dfrac{P}{L}\right) is what makes this model different from exponential growth:

  • When PP is small compared with LL, the factor is close to 11 and dPdt≈kP\dfrac{dP}{dt} \approx kP. Growth looks exponential.
  • As PP approaches LL, the factor approaches 00 and growth slows to a stop.
  • If PP is larger than LL, the factor is negative, so PP decreases back toward LL.

You will also see the equivalent form dPdt=cP(L−P)\dfrac{dP}{dt} = cP(L - P). Here LL is still the carrying capacity, and c=kLc = \dfrac{k}{L}. To find LL from any logistic equation, factor out PP and find the value of PP that makes the remaining factor zero.

Equilibrium solutions and long-run behavior

Setting dPdt=0\dfrac{dP}{dt} = 0 gives P=0P = 0 or P=LP = L. These constant functions are the equilibrium solutions. Every other solution with P(0)>0P(0) \gt 0 heads toward LL:

starting valuesign of dPdt\dfrac{dP}{dt}behaviorlim⁡t→∞P(t)\displaystyle\lim_{t \to \infty} P(t)
P(0)=0P(0) = 000stays at 0000
0<P(0)<L0 \lt P(0) \lt Lpositiveincreases toward LLLL
P(0)=LP(0) = L00stays at LLLL
P(0)>LP(0) \gt Lnegativedecreases toward LLLL

Where the population grows fastest

The rate dPdt=kP−kLP2\dfrac{dP}{dt} = kP - \dfrac{k}{L}P^2 is a downward-opening quadratic in PP with zeros at P=0P = 0 and P=LP = L. A parabola reaches its maximum halfway between its zeros, so the growth rate is largest when P=L2P = \dfrac{L}{2}.

You can also see this through the second derivative. Differentiating with respect to tt (using the chain rule, since PP depends on tt):

d2Pdt2=(k−2kLP)dPdt.\frac{d^2P}{dt^2} = \left(k - \frac{2k}{L}P\right)\frac{dP}{dt}.

For 0<P<L0 \lt P \lt L, the factor dPdt\dfrac{dP}{dt} is positive, so the sign of d2Pdt2\dfrac{d^2P}{dt^2} matches the sign of k−2kLPk - \dfrac{2k}{L}P: positive when P<L2P \lt \dfrac{L}{2} and negative when P>L2P \gt \dfrac{L}{2}.

Reading a logistic equation

For dPdt=kP(1−PL)\dfrac{dP}{dt} = kP\left(1 - \dfrac{P}{L}\right) with 0<P(0)<L0 \lt P(0) \lt L:

  • PP is increasing, and lim⁡t→∞P(t)=L\displaystyle\lim_{t \to \infty} P(t) = L.
  • PP is growing fastest when P=L2P = \dfrac{L}{2}. This is the inflection point of the solution curve.
  • The graph is concave up while P<L2P \lt \dfrac{L}{2} and concave down while P>L2P \gt \dfrac{L}{2}, which gives the familiar S-shape.

Worked example: Reading the equation

A fish population satisfies dPdt=0.4P(1−P2000)\dfrac{dP}{dt} = 0.4P\left(1 - \dfrac{P}{2000}\right) with P(0)=300P(0) = 300. Find the carrying capacity, the population when it is growing fastest, and the growth rate when P=500P = 500.

Solution. The carrying capacity is L=2000L = 2000 fish. Growth is fastest at P=20002=1000P = \dfrac{2000}{2} = 1000. When P=500P = 500,

dPdt=0.4(500)(1−5002000)=200(0.75)=150 fish per unit of time.\frac{dP}{dt} = 0.4(500)\left(1 - \frac{500}{2000}\right) = 200(0.75) = 150 \text{ fish per unit of time}.

Worked example: The other form

A rumor spreads through a school of 800 students according to dydt=0.003 y (800−y)\dfrac{dy}{dt} = 0.003\,y\,(800 - y), where yy is the number of students who have heard it after tt hours. What is the maximum rate at which the rumor spreads?

Solution. The factor 800−y800 - y is zero at y=800y = 800, so L=800L = 800. The rate is greatest when y=400y = 400:

dydt=0.003(400)(400)=480 students per hour.\frac{dy}{dt} = 0.003(400)(400) = 480 \text{ students per hour}.

Rewriting as 0.003(800) y(1−y800)0.003(800)\,y\left(1 - \dfrac{y}{800}\right) shows k=2.4k = 2.4.

The logistic solution

The logistic equation is separable (with partial fractions, a BC technique), and its solution is

P(t)=L1+Ae−kt,where A=L−P(0)P(0).P(t) = \frac{L}{1 + Ae^{-kt}}, \qquad \text{where } A = \frac{L - P(0)}{P(0)}.

You can check the value of AA: at t=0t = 0, P(0)=L1+AP(0) = \dfrac{L}{1 + A}, so 1+A=LP(0)1 + A = \dfrac{L}{P(0)}. As t→∞t \to \infty, e−kt→0e^{-kt} \to 0 and P(t)→LP(t) \to L, matching everything above.

P(t) = 1000/(1 + 9e^(−0.5t)) rises in an S-shape toward the carrying capacity 1000. It grows fastest at P = 500, when t = 2 ln 9 ≈ 4.39.Open in grapher →

Worked example: Using the solution

A population satisfies dPdt=0.5P(1−P1000)\dfrac{dP}{dt} = 0.5P\left(1 - \dfrac{P}{1000}\right) with P(0)=100P(0) = 100. Write P(t)P(t) and find when the population is growing fastest.

Solution. Here L=1000L = 1000, k=0.5k = 0.5, and A=1000−100100=9A = \dfrac{1000 - 100}{100} = 9, so

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

Growth is fastest when P=500P = 500: 1+9e−0.5t=21 + 9e^{-0.5t} = 2, so e−0.5t=19e^{-0.5t} = \dfrac{1}{9}, and t=2ln⁡9≈4.394t = 2\ln 9 \approx 4.394.

Common mistake

The fastest growth happens at P=L2P = \dfrac{L}{2}, which is a population value, not a time. If a question asks when growth is fastest, you need the solution formula (or given data) to turn P=L2P = \dfrac{L}{2} into a value of tt. And if P(0)>L2P(0) \gt \dfrac{L}{2}, the population never passes through L2\dfrac{L}{2}, so the growth rate is largest at the start.

Tip

A quick way to spot a logistic equation: dPdt\dfrac{dP}{dt} is a quadratic in PP with no constant term and a negative P2P^2 coefficient, like dPdt=0.2P−0.0005P2\dfrac{dP}{dt} = 0.2P - 0.0005P^2. Factor out PP to read off the carrying capacity: 0.2P(1−P400)0.2P\left(1 - \dfrac{P}{400}\right) gives L=400L = 400.

Practice

Practice 1

A population satisfies dPdt=0.05P(1−P1200)\dfrac{dP}{dt} = 0.05P\left(1 - \dfrac{P}{1200}\right). What is the carrying capacity?

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 2

A population satisfies dPdt=0.002P(900−P)\dfrac{dP}{dt} = 0.002P(900 - P) with P(0)=50P(0) = 50. For what value of PP is the population growing fastest?

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 3

For the population in the previous problem, dPdt=0.002P(900−P)\dfrac{dP}{dt} = 0.002P(900 - P), what is the maximum value of dPdt\dfrac{dP}{dt}?

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 4

A population satisfies dPdt=0.3P(1−P1000)\dfrac{dP}{dt} = 0.3P\left(1 - \dfrac{P}{1000}\right) with P(0)=1500P(0) = 1500. Find lim⁡t→∞P(t)\displaystyle\lim_{t \to \infty} P(t).

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 5

Let P(t)P(t) satisfy dPdt=0.2P(1−P1000)\dfrac{dP}{dt} = 0.2P\left(1 - \dfrac{P}{1000}\right) with P(0)=300P(0) = 300. Which describes the graph of PP for t≥0t \ge 0?

Practice 6

Find the particular solution P(t)P(t) of dPdt=0.2P(1−P600)\dfrac{dP}{dt} = 0.2P\left(1 - \dfrac{P}{600}\right) with P(0)=100P(0) = 100.

Enter an expression, e.g. 3x^2 - 2x + 1

Practice 7

For P(t)=6001+5e−0.2tP(t) = \dfrac{600}{1 + 5e^{-0.2t}} from the previous problem, at what time tt is the population growing fastest? Give an exact answer or a decimal to the nearest hundredth.

Enter a number. Fractions like 3/4 and sqrt(2) are OK.

Practice 8

A population satisfies dPdt=P5−P24000\dfrac{dP}{dt} = \dfrac{P}{5} - \dfrac{P^2}{4000} with P(0)=20P(0) = 20. Find lim⁡t→∞P(t)\displaystyle\lim_{t \to \infty} P(t).

Enter a number. Fractions like 3/4 and sqrt(2) are OK.