Skip to main content

First-order differential equations and separation of variables (QCE Specialist Mathematics Unit 4)

Syllabus dot point

“Formulate and solve first-order differential equations using separation of variables, including growth and decay and the logistic model, and interpret solutions in applied rates-of-change contexts”

8 min read

Quick answer

A first-order differential equation relates yy to its rate of change. If it factors as dydx=f(x)g(y)\frac{dy}{dx} = f(x)g(y), separate the variables and integrate both sides, adding one constant, then apply the initial condition for the particular solution. The model dydt=ky\frac{dy}{dt} = ky gives y=Aekty = Ae^{kt} (growth or decay); Newton's cooling gives T=Ts+(T0−Ts)e−ktT = T_s + (T_0 - T_s)e^{-kt}; the logistic model dPdt=rP(1−P/M)\frac{dP}{dt} = rP(1 - P/M) gives S-shaped growth toward carrying capacity MM.

Jump to a section
  1. What this dot point is asking
  2. The answer
  3. Exam-style questions

What this dot point is asking

QCAA wants you to set up a first-order differential equation from a worded rates-of-change situation, solve it by separation of variables, apply an initial condition to find the particular solution, and interpret the result. This dot point underpins much of IA1 modelling and is examined in IA3 and the external assessment, where growth, decay, cooling and logistic models recur.

The answer

What a differential equation says

A differential equation relates a quantity to its rate of change. The statement "the rate of growth is proportional to the current amount" becomes

dydt=ky.\frac{dy}{dt} = ky.

Solving means finding the function y(t)y(t) whose derivative matches the equation. The general solution contains an arbitrary constant; an initial condition pins it down to a particular solution.

Separation of variables

If a first-order equation can be written as dydx=f(x)g(y)\dfrac{dy}{dx} = f(x)g(y), separate the variables so each side involves only one variable, then integrate both sides:

∫1g(y) dy=∫f(x) dx.\int \frac{1}{g(y)}\,dy = \int f(x)\,dx.

Treating dydx\dfrac{dy}{dx} as a ratio of differentials is a valid shortcut for this technique. Include a single constant of integration, then solve algebraically for yy where possible.

Exponential growth and decay

The equation dydt=ky\dfrac{dy}{dt} = ky separates to ∫1y dy=∫k dt\int \frac{1}{y}\,dy = \int k\,dt, giving ln⁡∣y∣=kt+C\ln|y| = kt + C and hence

y=Aekt,A=y(0).y = A e^{kt}, \qquad A = y(0).

Positive kk gives growth; negative kk gives decay (for example radioactive decay, where the half-life is t1/2=ln⁡2∣k∣t_{1/2} = \dfrac{\ln 2}{|k|}).

Newton's law of cooling

The rate of cooling is proportional to the temperature difference from the surroundings TsT_s:

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

Substituting u=T−Tsu = T - T_s reduces this to exponential decay, giving T=Ts+(T0−Ts)e−ktT = T_s + (T_0 - T_s)e^{-kt}. The temperature approaches the ambient TsT_s as t→∞t \to \infty.

The logistic model

When growth is limited by a carrying capacity MM, the model is

dPdt=rP(1−PM).\frac{dP}{dt} = rP\left(1 - \frac{P}{M}\right).

Growth is near-exponential when PP is small and slows to zero as PP approaches MM. Solving by separation (using partial fractions on the left) produces an S-shaped curve that levels off at P=MP = M. The logistic model is a favourite IA1 modelling context because it captures realistic saturation.

General versus particular solutions

Integrating a first-order equation introduces exactly one arbitrary constant, and the family of curves this produces is the general solution. An initial or boundary condition (a known value of yy at a particular xx or tt) fixes the constant and selects the single particular solution that fits the situation. The order matters: keep the constant through the algebra, apply the condition once you have isolated it, and only then simplify. Applying the condition too early, or absorbing it into an exponential before evaluating, is the source of most lost marks.

Velocity, acceleration and the inverse-tangent integral

Specialist Mathematics often phrases differential equations in a kinematics context, where acceleration is dvdt\dfrac{dv}{dt} or vdvdxv\dfrac{dv}{dx}. Separating such an equation can produce standard integrals beyond logarithms, in particular ∫1a2+v2 dv=1aarctan⁡va+C\displaystyle\int \dfrac{1}{a^2 + v^2}\,dv = \dfrac{1}{a}\arctan\dfrac{v}{a} + C and ∫1a2−v2 dv=arcsin⁡va+C\displaystyle\int \dfrac{1}{\sqrt{a^2 - v^2}}\,dv = \arcsin\dfrac{v}{a} + C. Recognising which standard form the separated integrand matches is the planning step, and choosing aa to match the constant under the square or the constant added to v2v^2 completes it.

Checking a solution

Always verify a solution by substituting it back into the original equation. Differentiate your y(t)y(t) and confirm it reproduces the stated rate, and check that the initial condition is satisfied. A second sanity check is the long-term behaviour: a decay model should tend to zero, a cooling model should tend to the ambient temperature, and a logistic model should tend to the carrying capacity. If the long-run behaviour contradicts the physical setup, a sign error in kk or a dropped constant is the usual cause.

Standard solution forms

dydt=ky ⇒ y=Aekt,dTdt=−k(T−Ts) ⇒ T=Ts+(T0−Ts)e−kt.\frac{dy}{dt} = ky \ \Rightarrow\ y = Ae^{kt}, \qquad \frac{dT}{dt} = -k(T - T_s) \ \Rightarrow\ T = T_s + (T_0 - T_s)e^{-kt}.

Worked example

A population grows so that dPdt=0.2P\dfrac{dP}{dt} = 0.2P, with P=500P = 500 at t=0t = 0 (years). Find P(t)P(t) and the time for the population to reach 20002000.

The rate is proportional to the current population, which is the classic exponential growth form. We solve by separation of variables.

Step 1: Separate the variables and integrate

Move PP to the left side and tt to the right, then integrate both sides. One constant of integration on the right side is sufficient:

∫1P dP=∫0.2 dt  ⇒  ln⁡∣P∣=0.2t+C.\int \frac{1}{P}\,dP = \int 0.2\,dt \;\Rightarrow\; \ln|P| = 0.2t + C.

Step 2: Solve for PP

Exponentiate both sides to undo the logarithm. The constant eCe^C is just another positive constant, which we call AA:

P=Ae0.2t.P = A e^{0.2t}.

Step 3: Apply the initial condition

At t=0t = 0, P=500P = 500, which fixes AA uniquely:

500=Ae0=A,soP(t)=500 e0.2t.500 = A e^{0} = A, \quad \text{so} \quad P(t) = 500\, e^{0.2t}.

Step 4: Find when P=2000P = 2000

Set the formula equal to 20002000 and solve for tt by taking logarithms:

2000=500 e0.2t  ⇒  4=e0.2t  ⇒  0.2t=ln⁡4.2000 = 500\,e^{0.2t} \;\Rightarrow\; 4 = e^{0.2t} \;\Rightarrow\; 0.2t = \ln 4.

So t=ln⁡40.2=1.38630.2≈6.93t = \dfrac{\ln 4}{0.2} = \dfrac{1.3863}{0.2} \approx 6.93 years.

Step 5: Check the answer

Substitute t≈6.93t \approx 6.93 back into the formula:

P(6.93)=500 e0.2×6.93=500 e1.386=500×4.00=2000.P(6.93) = 500\,e^{0.2 \times 6.93} = 500\,e^{1.386} = 500 \times 4.00 = 2000.

The check confirms the solution is correct.

Final answer: P(t)=500 e0.2tP(t) = 500\,e^{0.2t}, and the population reaches 2000 after approximately 6.936.93 years.

Common trap

The most penalised error is omitting or mishandling the constant of integration. You need exactly one constant; introduce it on integration and then use the initial condition before, not after, simplifying away the absolute value. A second error is dropping the negative sign in a decay or cooling model, which turns decay into growth. Always check the long-term behaviour of your solution against the physical situation: cooling should approach ambient, decay should approach zero.

Exam-style questions

Questions in the style of QCAA 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.

Practice5 marks
Paper 2 (complex familiar). An object moves so that dvdt=−(4+v2)\dfrac{dv}{dt} = -(4 + v^2), where vv is velocity (m s−1^{-1}) and t≥0t \geq 0 (s). (a) Determine the general solution. The initial velocity is 1.51.5 m s−1^{-1}. (b) Determine the time when the object comes to rest.
Show worked answer →

(a) Separate: dv4+v2=−dt\dfrac{dv}{4 + v^2} = -dt. Using ∫1a2+v2 dv=1aarctan⁡va\displaystyle\int \dfrac{1}{a^2 + v^2}\,dv = \dfrac{1}{a}\arctan\dfrac{v}{a} with a=2a = 2: 12arctan⁡v2=−t+c\tfrac{1}{2}\arctan\dfrac{v}{2} = -t + c. This is the general solution.

(b) Apply v=1.5v = 1.5 at t=0t = 0: c=12arctan⁡(0.75)≈0.32c = \tfrac{1}{2}\arctan(0.75) \approx 0.32. At rest v=0v = 0: 12arctan⁡0=−t+0.32\tfrac{1}{2}\arctan 0 = -t + 0.32, so t≈0.32t \approx 0.32 s.

Markers reward the separation, the inverse-tangent standard integral, applying the initial condition, and solving for tt.

Practice4 marks
Paper 1 (technique). A radioactive sample decays so that dNdt=−0.05N\dfrac{dN}{dt} = -0.05N, with N=800N = 800 at t=0t = 0 (years). (a) Determine N(t)N(t). (b) Determine the half-life of the sample, to the nearest year.
Show worked answer →

(a) Separate and integrate: ln⁡∣N∣=−0.05t+C\ln|N| = -0.05t + C, so N=Ae−0.05tN = Ae^{-0.05t}. With N(0)=800N(0) = 800, A=800A = 800, giving N(t)=800e−0.05t.N(t) = 800e^{-0.05t}.

(b) Half-life solves 400=800e−0.05t400 = 800e^{-0.05t}, so e−0.05t=0.5e^{-0.05t} = 0.5, −0.05t=ln⁡0.5-0.05t = \ln 0.5, t=ln⁡20.05≈13.86t = \dfrac{\ln 2}{0.05} \approx 13.86, so about 1414 years.

Markers reward the general solution form, applying the initial condition, and the half-life calculation.

Practise this

Sources & how we know this

ExamExplained