BETAInsomnia StudyroomLogistic Models with Differential Equations
Logistic Models with Differential Equations

AP Course

AP Calculus AB/BC

Study the complete College Board sequence for AP Calculus AB and BC, from limits through infinite series.

Choose an official unit and topic to open its lecture, concept check, or focused practice.

AP Calculus AB/BC · Unit 7 · Topic 7.9 · BC Only

Logistic Models with Differential Equations

Analyze growth limited by a carrying capacity and identify equilibria and maximum growth.

1. Topic Focus

Model rates with differential equations, read slope fields, approximate solutions, solve separable equations, and interpret exponential or logistic models.

This topic: Analyze growth limited by a carrying capacity and identify equilibria and maximum growth.

2. Key Relationship

\(\frac{dP}{dt}=kP\left(1-\frac PL\right)\)

Read every symbol with its domain, direction, units, and hypotheses before applying the relationship.

3. Visual Connection

carrying capacity LL/2
Logistic growthGrowth is fastest at half the carrying capacity, then slows as the solution approaches the upper equilibrium.

4. Worked Example

Equilibria occur at P=0 and P=L; growth is fastest at P=L/2.

Write the governing relationship first, carry out the algebra cleanly, and finish with a sentence that answers the mathematical question.

5. Concept Development

AP Calculus BC topic: Logistic differential equations extend exponential models by making the relative growth rate decrease as the quantity approaches a limiting capacity.

The logistic model

A standard logistic initial-value problem is

\(\frac{dP}{dt}=rP\left(1-\frac{P}{K}\right),\qquad r>0,\quad K>0,\quad P(0)=P_0.\)
SymbolMeaning
\(P(t)\)Quantity or population at time \(t\)
\(K\)Carrying capacity, the long-run sustainable level
\(r\)Intrinsic relative growth constant, with units of reciprocal time
\(P_0\)Initial quantity

The factor \(rP\) produces exponential-like growth. The crowding factor \(1-P/K\) reduces that growth as \(P\) increases:

\(\frac{P'}{P}=r\left(1-\frac{P}{K}\right).\)

Thus the relative growth rate decreases linearly with population. When \(P\ll K\), the crowding factor is near \(1\), so the model initially behaves approximately like \(P'=rP\).

Equilibria and qualitative behavior

Setting \(P'=0\) gives equilibrium solutions \(P=0\) and \(P=K\). A sign analysis determines the direction of every nonnegative solution.

PopulationSign of \(P'\)Behavior
\(P=0\)\(0\)Constant equilibrium
\(0PositiveIncreases toward \(K\)
\(P=K\)\(0\)Constant equilibrium
\(P>K\)NegativeDecreases toward \(K\)

In the nonnegative model, \(K\) is a stable equilibrium: nearby solutions move toward it. The equilibrium \(0\) is unstable because a small positive population grows away from it. Carrying capacity is therefore a limiting level, not necessarily a hard upper bound; a population that begins above \(K\) decreases toward it.

Solving by separation of variables

For a non-equilibrium positive solution, separate and use partial fractions:

\(\frac{dP}{P(1-P/K)}=r\,dt,\qquad \frac{K}{P(K-P)}=\frac1P+\frac1{K-P}.\)

Integration gives

\(\ln|P|-\ln|K-P|=rt+C.\)

Solving for \(P\) and applying \(P(0)=P_0>0\) produces

\(P(t)=\frac{K}{1+Ae^{-rt}},\qquad A=\frac{K-P_0}{P_0}.\)

The solution \(P=0\) must be listed separately because the separation step divides by \(P\). If \(P_0=K\), then \(A=0\) and the formula correctly gives the constant solution \(P=K\).

Concavity and fastest growth

Differentiate the differential equation with respect to time:

\(P''=rP'\left(1-\frac{2P}{K}\right).\)

For \(0

\(G(P)=rP\left(1-\frac{P}{K}\right)=rP-\frac rK P^2.\)

This downward-opening parabola has its vertex at \(P=K/2\), so

\(\text{maximum growth rate}=G\left(\frac K2\right)=\frac{rK}{4}.\)

If \(0

\(t_*=\frac1r\ln\left(\frac{K-P_0}{P_0}\right).\)

If \(P_0>K/2\), that inflection occurred before \(t=0\), so the future solution is already concave down. A solution starting above \(K\) decreases and is concave up as it approaches \(K\).

Estimating a parameter from data

If \(K\), \(P_0\), and a later value \(P(t_1)=P_1\) are known with \(0

\(r=\frac1{t_1}\ln\left(\frac{P_1(K-P_0)}{P_0(K-P_1)}\right).\)

Always check units and assumptions. The basic model treats \(r\) and \(K\) as constant and ignores migration, delays, seasonal effects, and sudden environmental changes. A good AP response interprets conclusions within the model rather than claiming that a real population must follow the curve forever.

6. Detailed Worked Example and Error Check

Example 1: Read the model

For

\(P'=0.6P\left(1-\frac{P}{800}\right),\)

the intrinsic growth constant is \(r=0.6\), the carrying capacity is \(K=800\), and the equilibria are \(0\) and \(800\). Growth is fastest at \(P=400\), with rate

\(\frac{rK}{4}=\frac{(0.6)(800)}4=120\)

individuals per time unit.

Example 2: Solve an initial-value problem

Suppose \(P'=0.4P(1-P/1000)\) and \(P(0)=100\). Then

\(A=\frac{1000-100}{100}=9,\qquad P(t)=\frac{1000}{1+9e^{-0.4t}}.\)

After \(10\) time units,

\(P(10)=\frac{1000}{1+9e^{-4}}\approx858.49.\)

The value remains below \(1000\) and approaches it over time.

Example 3: Verify the explicit solution

Let \(P=K(1+Ae^{-rt})^{-1}\). Differentiating gives

\(P'=\frac{KrAe^{-rt}}{(1+Ae^{-rt})^2}.\)

Also,

\(rP\left(1-\frac PK\right)=r\frac{K}{1+Ae^{-rt}}\left(\frac{Ae^{-rt}}{1+Ae^{-rt}}\right)=\frac{KrAe^{-rt}}{(1+Ae^{-rt})^2}.\)

The two expressions agree, so the formula satisfies the logistic differential equation wherever it is defined.

Example 4: Begin above carrying capacity

For \(P'=0.2P(1-P/500)\) and \(P(0)=700\),

\(A=\frac{500-700}{700}=-\frac27,\qquad P(t)=\frac{500}{1-(2/7)e^{-0.2t}}.\)

Because \(P>500\), \(P'<0\). Also \(P''=0.2P'(1-2P/500)>0\), so the population decreases, is concave up, and approaches \(500\).

Example 5: Locate the inflection point

For the model in Example 2, fastest growth occurs at \(P=500\). Solve

\(500=\frac{1000}{1+9e^{-0.4t}}\)

to obtain

\(t=\frac{\ln9}{0.4}\approx5.49.\)

At that time the growth rate is \(rK/4=(0.4)(1000)/4=100\). The population is not largest there; its rate of increase is largest there.

Example 6: Estimate the growth constant

A population has \(K=1000\), \(P(0)=100\), and \(P(5)=300\). Since \(A=9\),

\(300=\frac{1000}{1+9e^{-5r}}\quad\Longrightarrow\quad e^{-5r}=\frac7{27}.\)

Therefore

\(r=\frac15\ln\left(\frac{27}{7}\right)\approx0.270\)

per time unit.

Example 7: Compare relative and absolute growth

For \(P'=0.5P(1-P/10000)\), at \(P=100\) the relative growth rate is

\(\frac{P'}P=0.5(0.99)=0.495,\)

which is close to the exponential rate \(0.5\). At \(P=5000\), the relative rate has fallen to \(0.25\), although the absolute rate is maximal there.

Example 8: Read the rate as a function of population

For \(P'=0.4P(1-P/1000)\), the rates at \(P=0,250,500,750,1000\) are respectively

\(0,\quad75,\quad100,\quad75,\quad0.\)

The equal rates at \(250\) and \(750\) reflect the symmetry of the quadratic \(G(P)\), and its vertex confirms the maximum at \(P=500\).

Common errors

  • Calling \(K/2\) the maximum population instead of the population at maximum growth.
  • Treating \(r\) as the actual relative rate at every population; the relative rate is \(r(1-P/K)\).
  • Forgetting the equilibrium solution \(P=0\) after dividing by \(P\).
  • Assuming every logistic curve is S-shaped for \(t\ge0\); the visible concavity depends on the initial value.
  • Using exponential growth indefinitely even when a carrying capacity is part of the context.

7. AP Reasoning Routine

Translate the context into a rate equation, verify candidate solutions by substitution, carry constants through integration, and apply initial conditions last.

  • Identify the representation and requested quantity.
  • State the rule or theorem and verify its conditions.
  • Keep exact values until the final requested approximation.
  • Interpret sign, units, interval, and context.
AP Checkpoint

Analyze each logistic model. Give exact values before decimal approximations.
(a) For \(P'=0.3P(1-P/1200)\), identify \(r\), \(K\), the equilibria, and the sign of \(P'\) on the nonnegative intervals they determine.
(b) In the model from part (a), at what population is growth fastest, and what is the maximum growth rate?
(c) Solve \(P'=0.2P(1-P/500)\), \(P(0)=50\).
(d) Use the solution from part (c) to find \(P(10)\).
(e) When does the solution from part (c) reach its inflection point?
(f) Describe the direction, concavity, and long-run behavior of a solution to \(P'=0.1P(1-P/500)\) with \(P(0)=700\).
(g) For \(P'=0.4P(1-P/1000)\), find both the relative and absolute growth rates when \(P=200\).
(h) A logistic population has \(K=800\), \(P(0)=80\), and \(P(4)=200\). Find \(r\).
(i) Explain one mathematical difference between exponential and logistic relative growth rates, and state when the logistic model resembles exponential growth.
(j) State the solutions for the initial conditions \(P(0)=0\) and \(P(0)=K\). Explain why one requires special care during separation of variables.

Check the solution

(a) \(r=0.3\), \(K=1200\), and the equilibria are \(P=0,1200\). For \(00\); for \(P>1200\), \(P'<0\).
(b) Growth is fastest at \(P=K/2=600\), and the maximum rate is \(rK/4=(0.3)(1200)/4=90\).
(c) \(A=(500-50)/50=9\), so \(P(t)=500/(1+9e^{-0.2t})\).
(d) \(P(10)=500/(1+9e^{-2})\approx225.43\).
(e) The inflection occurs at \(P=250\), so \(t=(\ln9)/0.2=5\ln9\approx10.99\).
(f) Since \(P>500\), the solution decreases and is concave up. It approaches the stable equilibrium \(P=500\) from above.
(g) \(P'/P=0.4(1-200/1000)=0.32\) per time unit, and \(P'=0.32(200)=64\) individuals per time unit.
(h) Here \(A=(800-80)/80=9\). Since \(200=800/(1+9e^{-4r})\), \(e^{-4r}=1/3\), so \(r=(\ln3)/4\approx0.275\).
(i) Exponential growth has constant relative rate \(P'/P=r\). Logistic growth has decreasing relative rate \(P'/P=r(1-P/K)\). When \(P\) is small relative to \(K\), the latter is approximately \(r\), so the models initially look similar.
(j) The solutions are \(P(t)=0\) and \(P(t)=K\), respectively. Separation divides by \(P\), so it loses the zero solution unless that equilibrium is recorded first.