Skip to content

Logistic growth and phase lines

For an equation dydt=f(y)\dfrac{dy}{dt} = f(y), the steady states are the zeros of ff, and the sign of ff between them says which way every solution moves. Drawn as arrows on the yy-axis, this state-space diagram (phase line) tells you what a population does in the long run without solving anything. The logistic equation is the standard case: growth that slows as the population nears its carrying capacity.

In this topic

Taught
Before the final, flavour B
Questions
11
Total marks
32
Suggested time
57 min
Difficulty

Introductory to Capstone

Show the 3 hardest

Key results

Steady states
f(y∗)=0  ⟺  y(t)=y∗ is a constant solutionf(y^{*}) = 0 \iff y(t) = y^{*} \text{ is a constant solution}

For an equation dydt=f(y)\dfrac{dy}{dt} = f(y) whose right side depends on yy alone.

Stability
f ⁣:+ below, − above  ⟹  y∗ stablef ⁣:− below, + above  ⟹  y∗ unstable\begin{aligned} f \colon + \text{ below},\ - \text{ above} &\implies y^{*} \text{ stable} \\ f \colon - \text{ below},\ + \text{ above} &\implies y^{*} \text{ unstable} \end{aligned}

For a steady state y∗y^{*} of dydt=f(y)\dfrac{dy}{dt} = f(y), with the signs of ff taken on intervals just below and just above y∗y^{*}.

The logistic equation
dNdt=rN(1−NK)\frac{dN}{dt} = rN\Bigl(1 - \frac{N}{K}\Bigr)

For constants r>0r > 0 and K>0K > 0: the steady states are 0 (unstable) and KK (stable), and NN grows fastest at N=K/2N = K/2, at the rate rK/4rK/4.

Solution of the logistic equation
N(t)=K1+Ae−rt,A=K−N(0)N(0)N(t) = \frac{K}{1 + A e^{-rt}}, \qquad A = \frac{K - N(0)}{N(0)}

For N(0)>0N(0) > 0 and t≥0t \ge 0.

Harvesting
dNdt=rN(1−NK)−h has a steady state  ⟺  h≤rK4\frac{dN}{dt} = rN\Bigl(1 - \frac{N}{K}\Bigr) - h \text{ has a steady state} \iff h \le \frac{rK}{4}

For constants r>0r > 0, K>0K > 0 and h≥0h \ge 0; rK/4rK/4 is the maximum sustainable yield.

Read the theoryUBC Math 100 §12.1 Linear and nonlinear differential equationsUBC Math 100 §13.1 State-space diagramsUBC Math 100 §13.2 Logistic growth

Worked example

Standard3 marks · 5 min

A fish stock of N(t)N(t) tonnes grows logistically and is harvested at a constant 32 tonnes per year, so that the equation below holds. Find the steady states. If the stock is 150 tonnes when harvesting begins, what does it approach in the long run?

dNdt=0.4N(1−N500)−32\frac{dN}{dt} = 0.4N\left(1 - \frac{N}{500}\right) - 32

Give the long-run stock in tonnes.

Strategy

Find the steady states as the zeros of the right-hand side, then use its sign between them to see which way a stock of 150 tonnes moves and where it stops.

Working

  1. 0.4N−0.0008N2−32=0  ⟺  N2−500N+40 000=00.4N - 0.0008N^{2} - 32 = 0 \iff N^{2} - 500N + 40\,000 = 0

    Multiply through by −1250 to clear the decimals.

  2. (N−100)(N−400)=0⟹N=100 or N=400(N - 100)(N - 400) = 0 \quad\Longrightarrow\quad N = 100 \text{ or } N = 400

    Harvesting moves the steady states of the unfished stock, 0 and 500, to 100 and 400.

  3. dNdt=−0.0008 (N−100)(N−400)\frac{dN}{dt} = -0.0008\,(N - 100)(N - 400)

    Negative below 100, positive between 100 and 400, negative above 400: 400 is stable and 100 is unstable.

  4. N=150:dNdt=0.4(150)(0.7)−32=10>0N = 150:\quad \frac{dN}{dt} = 0.4(150)(0.7) - 32 = 10 > 0

    The stock starts between the steady states, so it rises and approaches 400.

Answer

Steady states 100 and 400;N(t)→400 tonnes\text{Steady states } 100 \text{ and } 400;\quad N(t) \to 400 \text{ tonnes}

The answer is steady states at 100 and 400 tonnes; the stock approaches 400 tonnes.

Check

At 400 tonnes growth 0.4 · 400 · 0.2 = 32 a year replaces the catch, and the stock gains 10 tonnes a year at 350 but loses 14 at 450.

Where this goes wrong

Choosing the nearer steady state, 100. It is unstable: just above 100 the stock grows, so from 150 tonnes it climbs away from 100 to the stable steady state at 400.

On the paper

Find the equilibria by setting the rate to zero, draw the phase line with the sign of the rate between them, and classify each equilibrium from the arrows.

Leans on

Solving a quadratic equation, from Functions, graphs and inverses

Your turn

Challenging4 marks · 8 min

A fishery works with a fixed effort EE, so its catch is proportional to the stock: ENEN tonnes per year. The stock N(t)N(t), in tonnes, then satisfies the equation below, with 0<E<0.80 < E < 0.8. Find the positive steady state N∗N^{*} in terms of EE, and the catch EN∗EN^{*} taken each year once the stock is there. Which effort EE makes that sustained catch as large as possible?

dNdt=0.8N(1−N1000)−EN\frac{dN}{dt} = 0.8N\left(1 - \frac{N}{1000}\right) - EN

Give the effort E, per year.

Answer this question

Practice questions

5 of the 11, across every level. Open a solution to read the working, or answer the question yourself and it counts towards your progress.

  • Introductory1 mark · 2 min

    How many steady states does this differential equation have?

    dydt=y4−4y2\frac{dy}{dt} = y^{4} - 4y^{2}
    Answer this question
  • Introductory2 marks · 3 min

    The differential equation below has steady states y=1y = 1 and y=4y = 4. Which statement describes their stability?

    dydt=(y−1)(y−4)\frac{dy}{dt} = (y - 1)(y - 4)
    Answer this question
  • Standard2 marks · 4 min

    A patient receives a drug through an intravenous drip at 10 mg per hour. The liver removes 4 % of the drug in the body each hour, and the kidneys remove a further 1 % each hour. Write a differential equation for the amount A(t)A(t) of drug in the body, in mg, and find the amount it approaches in the long run.

    Give the amount in milligrams.

    Answer this question
  • Challenging4 marks · 7 min

    Left alone, a fish stock of N(t)N(t) tonnes would grow logistically with r=0.3r = 0.3 per year and carrying capacity 4000 tonnes. A fishery removes a constant hh tonnes per year, so the stock satisfies the equation below. Find the maximum sustainable yield: the largest hh for which the stock has a steady state. For any larger hh the stock collapses, whatever its starting size.

    dNdt=0.3N(1−N4000)−h\frac{dN}{dt} = 0.3N\left(1 - \frac{N}{4000}\right) - h

    Give h in tonnes per year.

    Answer this question
  • Capstone6 marks · 12 min

    A protein switches on its own production. Its concentration c(t)≥0c(t) \ge 0 satisfies the equation below, where the first term is production and −c-c is decay.

    dcdt=10c29+c2−c\frac{dc}{dt} = \frac{10c^{2}}{9 + c^{2}} - c

    (a) Find every steady state. (b) Draw the state-space diagram (phase line) for c ≥ 0 and classify each steady state as stable or unstable. (c) Describe what happens in the long run if c(0) = 0.5, and if c(0) = 3.

    Answer this question

Browse all 11 questions

  • Needs

    Slope fields and Euler’s method

    Checking that a function solves a differential equation, reading a slope field, and stepping a solution forward with Euler’s method.

    14 questions