Skip to main content

Modelling with differential equations (SACE Stage 2 Specialist Mathematics)

Syllabus dot point

“Formulate and solve differential equations modelling real situations such as growth, cooling and limited growth, and interpret the solutions.”

9 min read

Quick answer

Modelling turns a worded rate description into a differential equation. "Rate proportional to amount" gives dydt=ky\dfrac{dy}{dt} = ky (exponential growth or decay). "Rate proportional to the difference from a fixed value" gives Newton's law of cooling dTdt=−k(T−Ts)\dfrac{dT}{dt} = -k(T - T_s). "Limited growth toward a capacity" gives a logistic equation. Solve by separation, apply conditions, then interpret the long-term behaviour.

Jump to a section
  1. What this dot point is asking
  2. Translating words into a differential equation
  3. Exponential growth and decay
  4. Newton's law of cooling
  5. Limited (logistic) growth
  6. Interpreting the logistic curve
  7. Setting up from scratch
  8. Exam-style questions

What this dot point is asking

You need to formulate differential equations from real situations, solve them, and interpret the resulting solutions including their long-term behaviour.

Translating words into a differential equation

The phrase describing the rate dictates the equation:

Reading the rate phrase
  • "rate of change proportional to the amount present" gives dydt=ky\dfrac{dy}{dt} = ky.
  • "rate proportional to the difference between TT and the surroundings TsT_s" gives dTdt=−k(T−Ts)\dfrac{dT}{dt} = -k(T - T_s).
  • "growth limited by a carrying capacity MM" gives a logistic form such as dPdt=kP(M−P)\dfrac{dP}{dt} = kP(M - P).
    The proportionality constant kk is positive; signs are set by whether the quantity rises or falls.

Exponential growth and decay

The simplest model, dydt=ky\dfrac{dy}{dt} = ky, solves to y=Aekty = A e^{kt} (see the separable-equations method). Use two pieces of data: one to find AA, another to find kk. Interpretation: as t→∞t \to \infty the quantity grows without bound if k>0k > 0 or decays to 00 if k<0k < 0.

Newton's law of cooling

A body at temperature TT in surroundings at TsT_s cools at a rate proportional to the temperature difference:

dTdt=−k(T−Ts).\frac{dT}{dt} = -k(T - T_s).

Newton's law of cooling

A drink at 80∘80^{\circ}C cools in a 20∘20^{\circ}C room. After 55 minutes it is 60∘60^{\circ}C. Find T(t)T(t)

Step 1: Write the differential equation, separate variables and integrate

Newton's law of cooling says the rate of change of temperature is proportional to the gap between the object's temperature and the surrounding temperature. With room temperature 20∘20^{\circ}C the model is dTdt=−k(T−20)\dfrac{dT}{dt} = -k(T - 20). Separating variables and integrating both sides gives:

∫1T−20 dT=∫−k dt  ⟹  ln⁡∣T−20∣=−kt+C.\int \frac{1}{T - 20}\,dT = \int -k\,dt \implies \ln|T - 20| = -kt + C.

Step 2: Exponentiate and write the general solution

Exponentiating both sides removes the logarithm. Writing A=eCA = e^{C} absorbs the integration constant into a single parameter.

T−20=Ae−ktT - 20 = A e^{-kt}, so T=20+Ae−ktT = 20 + A e^{-kt}.

Step 3: Apply the initial condition T=80T = 80 at t=0t = 0 to find AA

Substituting t=0t = 0 into the general solution immediately gives AA, because e0=1e^{0} = 1.

80=20+A80 = 20 + A, so A=60A = 60. Thus T=20+60e−ktT = 20 + 60 e^{-kt}.

Step 4: Apply the second condition T=60T = 60 at t=5t = 5 to find kk

The second piece of data gives an equation in kk alone, which we solve by taking a logarithm.

60=20+60e−5k  ⟹  40=60e−5k  ⟹  e−5k=23.60 = 20 + 60 e^{-5k} \implies 40 = 60 e^{-5k} \implies e^{-5k} = \tfrac{2}{3}.

So −5k=ln⁡23-5k = \ln\tfrac{2}{3}, giving k=−15ln⁡23=15ln⁡32≈0.0811k = -\tfrac{1}{5}\ln\tfrac{2}{3} = \tfrac{1}{5}\ln\tfrac{3}{2} \approx 0.0811.

Step 5: Write the final model and interpret the long-term behaviour

Substituting the numerical value of kk gives the complete model. Examining what happens as time grows large confirms the physical meaning of the solution.

The model is T=20+60e−0.0811tT = 20 + 60 e^{-0.0811 t}. As t→∞t \to \infty, e−0.0811t→0e^{-0.0811t} \to 0, so T→20∘T \to 20^{\circ}C. The drink approaches room temperature, as expected.

Limited (logistic) growth

When growth slows as a population nears a carrying capacity MM, the logistic equation dPdt=kP(M−P)\dfrac{dP}{dt} = kP(M - P) applies. Its solution is an S-shaped curve that rises from a small initial value and levels off at MM. The key interpretation: dPdt→0\dfrac{dP}{dt} \to 0 as P→MP \to M, so the population stabilises at the carrying capacity, and the steepest growth occurs at the halfway point P=M2P = \tfrac{M}{2}.

Always interpret the long-term behaviour

After solving, examine t→∞t \to \infty. Decay and cooling models approach a limiting value; exponential growth is unbounded; logistic growth levels off at the capacity. SACE questions almost always ask you to state and interpret this limit.

Interpreting the logistic curve

The logistic solution P=M1+Ae−ktP = \dfrac{M}{1 + A e^{-kt}} has features SACE frequently asks you to read off. The carrying capacity MM is the horizontal asymptote the curve approaches as t→∞t \to \infty. The constant AA is fixed by the initial population: substituting t=0t = 0 gives P(0)=M1+AP(0) = \dfrac{M}{1 + A}, so A=M−P(0)P(0)A = \dfrac{M - P(0)}{P(0)}. The growth rate is slow at first when PP is small, fastest at the inflection point P=M2P = \tfrac{M}{2}, then slows again as PP nears MM. Being able to identify the capacity, the initial value, and the point of fastest growth from the formula is exactly what the interpretation marks reward.

Setting up from scratch

A worked setup: a tank gains salt at a constant rate and loses it proportionally to the amount present, giving dSdt=a−bS\dfrac{dS}{dt} = a - bS. This is separable and solves to an exponential approaching the equilibrium S=abS = \tfrac{a}{b}, where the rate of change is zero. Recognising the equilibrium (where dSdt=0\dfrac{dS}{dt} = 0) often answers the long-term question without fully solving.

Common errors
Wrong sign on kk
Cooling and decay need a negative rate, so the term is −k(T−Ts)-k(T - T_s) or −kN-kN with k>0k > 0. A sign error makes the model grow instead of decay.
Using only one data point
Most models have two unknowns (AA and kk); you need two conditions to find both.
Forgetting the surroundings term
Newton's cooling solution is T=Ts+Ae−ktT = T_s + A e^{-kt}, not Ae−ktA e^{-kt} alone; the TsT_s offset is essential and sets the limiting temperature.
Skipping the interpretation
Stating the limiting value and what it means (room temperature, carrying capacity, equilibrium) is required for full marks.

Exam-style questions

Questions in the style of SACE Board exam questions on this dot point, each with a worked answer. They are written by ExamExplained unless tagged "Past paper"; the year shows the paper a question is modelled on.

SACE 2023-style4 marks
Calculator-assumed. A butterfly population BB grows by dBdt=0.1B(1−BK)\dfrac{dB}{dt} = 0.1B\left(1 - \dfrac{B}{K}\right), where KK is a positive constant. You may assume KB(K−B)=1B+1K−B\dfrac{K}{B(K-B)} = \dfrac{1}{B} + \dfrac{1}{K-B}. Show that B=K1+Ae−0.1tB = \dfrac{K}{1 + A e^{-0.1t}} for some constant AA.
Show worked answer →

Rewrite as dBdt=0.1B(K−B)K\dfrac{dB}{dt} = \dfrac{0.1B(K-B)}{K} and separate: KB(K−B) dB=0.1 dt\dfrac{K}{B(K-B)}\,dB = 0.1\,dt. [1 mark]

Use the given partial fractions: ∫(1B+1K−B)dB=∫0.1 dt\int \left(\dfrac{1}{B} + \dfrac{1}{K-B}\right)dB = \int 0.1\,dt, giving ln⁡∣B∣−ln⁡∣K−B∣=0.1t+c\ln|B| - \ln|K-B| = 0.1t + c, i.e. ln⁡BK−B=0.1t+c\ln\dfrac{B}{K-B} = 0.1t + c. [1 mark]

Exponentiate: BK−B=Ce0.1t\dfrac{B}{K-B} = Ce^{0.1t}. Take reciprocals: K−BB=1Ce−0.1t\dfrac{K-B}{B} = \dfrac{1}{C}e^{-0.1t}, so KB−1=Ae−0.1t\dfrac{K}{B} - 1 = A e^{-0.1t} with A=1CA = \tfrac{1}{C}. [1 mark]

Then KB=1+Ae−0.1t\dfrac{K}{B} = 1 + A e^{-0.1t}, so B=K1+Ae−0.1tB = \dfrac{K}{1 + A e^{-0.1t}}, as required. [1 mark]

SACE 2022-style5 marks
Calculator-assumed. The area covered by bacteria satisfies dAdt=12A(50−A)50\dfrac{dA}{dt} = \tfrac{1}{2}A\dfrac{(50 - A)}{50} with A(0)=1A(0) = 1. Solve to show that A=501+49e−0.5tA = \dfrac{50}{1 + 49 e^{-0.5t}}.
Show worked answer →

Separate: 50A(50−A) dA=12 dt\dfrac{50}{A(50-A)}\,dA = \tfrac{1}{2}\,dt. Partial fractions: 50A(50−A)=1A+150−A\dfrac{50}{A(50-A)} = \dfrac{1}{A} + \dfrac{1}{50-A}. [1 mark]

Integrate: ln⁡∣A∣−ln⁡∣50−A∣=12t+c\ln|A| - \ln|50-A| = \tfrac{1}{2}t + c, so ln⁡A50−A=0.5t+c\ln\dfrac{A}{50-A} = 0.5t + c. [1 mark]

Exponentiate: A50−A=Ce0.5t\dfrac{A}{50-A} = Ce^{0.5t}. Apply A(0)=1A(0) = 1: 149=C\dfrac{1}{49} = C. [1 mark]

So A50−A=149e0.5t\dfrac{A}{50-A} = \tfrac{1}{49}e^{0.5t}, giving 49A=(50−A)e0.5t49A = (50-A)e^{0.5t}, hence A=50e0.5t49+e0.5tA = \dfrac{50e^{0.5t}}{49 + e^{0.5t}}. [1 mark]

Divide through by e0.5te^{0.5t}: A=501+49e−0.5tA = \dfrac{50}{1 + 49 e^{-0.5t}}, as required. [1 mark]

SACE 2021-style1 mark
Calculator-assumed. Given A=501+49e−0.5tA = \dfrac{50}{1 + 49 e^{-0.5t}} models the area covered by bacteria, state the maximum area available for growth.
Show worked answer →

As t→∞t \to \infty, e−0.5t→0e^{-0.5t} \to 0, so the denominator 1+49e−0.5t→11 + 49e^{-0.5t} \to 1.

Therefore A→50A \to 50. The maximum (limiting) area is 5050, the carrying capacity of the logistic model. [1 mark]

Practise this

Sources & how we know this

ExamExplained