Skip to main content

Differential equations: VCE Specialist Mathematics Unit 4

Syllabus dot point

“Formulation and solution of first-order differential equations including those solvable by direct integration and by separation of variables, the use of initial conditions to find particular solutions, and the interpretation of solutions in modelling contexts”

7 min read

Quick answer

A first-order differential equation relates a function to its first derivative. If dydx=f(x)\frac{\mathrm{d}y}{\mathrm{d}x} = f(x), solve by direct integration: y=∫f(x) dx+cy = \int f(x)\,\mathrm{d}x + c. If dydx=g(x)h(y)\frac{\mathrm{d}y}{\mathrm{d}x} = g(x)h(y), solve by separation of variables: rearrange to 1h(y) dy=g(x) dx\frac{1}{h(y)}\,\mathrm{d}y = g(x)\,\mathrm{d}x and integrate both sides. The constant of integration cc gives a family of solution curves; an initial condition such as y(x0)=y0y(x_0) = y_0 selects the one particular solution through that point. Always state the domain and interpret the answer for the modelling context.

Jump to a section
  1. What this dot point is asking
  2. What a differential equation is
  3. Solution by direct integration
  4. Solution by separation of variables
  5. Applying initial conditions
  6. Interpreting solutions in modelling
  7. Examples in context
  8. Try this
  9. Exam-style questions

What this dot point is asking

VCAA wants you to recognise and solve first-order differential equations, principally those solvable by direct integration (where dydx\frac{\mathrm{d}y}{\mathrm{d}x} is a function of xx alone) and by separation of variables (where the equation can be written with all the yy terms on one side and all the xx terms on the other). You also need to apply an initial condition to pin down the arbitrary constant and produce a particular solution, then interpret that solution in a modelling context.

What a differential equation is

A differential equation is an equation involving a function and its derivatives. A first-order equation involves only the first derivative. A solution is a function that satisfies the equation; the general solution contains an arbitrary constant and so describes a whole family of curves, while a particular solution fixes the constant using extra information (an initial or boundary condition).

The order of business is: identify the type, integrate to obtain the general solution, then apply the condition.

Solution by direct integration

When the derivative is a function of the independent variable only,

dydx=f(x),\frac{\mathrm{d}y}{\mathrm{d}x} = f(x),

the solution is found by antidifferentiating once:

y=∫f(x) dx+c.y = \int f(x)\,\mathrm{d}x + c.

Second derivatives are handled by integrating twice, introducing two constants. This is the simplest case and uses the integration techniques from elsewhere in Unit 4.

Solution by separation of variables

Many equations have the derivative depending on both variables in a product or quotient form:

dydx=g(x) h(y).\frac{\mathrm{d}y}{\mathrm{d}x} = g(x)\,h(y).

Treat dydx\frac{\mathrm{d}y}{\mathrm{d}x} as a ratio of differentials, gather all yy terms with dy\mathrm{d}y and all xx terms with dx\mathrm{d}x:

1h(y) dy=g(x) dx,\frac{1}{h(y)}\,\mathrm{d}y = g(x)\,\mathrm{d}x,

then integrate both sides. A single constant of integration suffices (combine the two). Where convenient you can leave the solution in implicit form, or solve for yy explicitly if the question asks.

Applying initial conditions

The general solution carries an arbitrary constant cc. An initial condition of the form y=y0y = y_0 when x=x0x = x_0 is substituted into the general solution to solve for cc, giving the particular solution. In modelling problems the condition is usually a known starting value, such as a population or temperature at time t=0t = 0.

Worked example

A separable equation with an initial condition

Solve dydx=xy\dfrac{\mathrm{d}y}{\mathrm{d}x} = \dfrac{x}{y} given that y=3y = 3 when x=0x = 0, taking y>0y > 0.

Step 1: Separate the variables

Multiply both sides by yy and by dx\mathrm{d}x:

y dy=x dx.y\,\mathrm{d}y = x\,\mathrm{d}x.

The goal of separation is to get all yy terms on one side and all xx terms on the other, so that each side can be integrated independently. Here, multiplying by y dxy\,\mathrm{d}x achieves this cleanly.

Step 2: Integrate both sides

∫y dy=∫x dx  ⇒  y22=x22+c.\int y\,\mathrm{d}y = \int x\,\mathrm{d}x \;\Rightarrow\; \frac{y^2}{2} = \frac{x^2}{2} + c.

Only one constant of integration is needed; constants from both sides can be combined into a single cc on the right. This gives the general solution as a family of curves.

Step 3: Apply the initial condition to find cc

Substitute x=0x = 0, y=3y = 3: 92=0+c\dfrac{9}{2} = 0 + c, so c=92c = \dfrac{9}{2}. Thus

y22=x22+92  ⇒  y2=x2+9.\frac{y^2}{2} = \frac{x^2}{2} + \frac{9}{2} \;\Rightarrow\; y^2 = x^2 + 9.

The initial condition selects the one particular solution from the family. Substituting the known values of xx and yy directly into the general solution is the most reliable method for finding cc.

Step 4: Solve for yy and verify

Since y>0y > 0 we take the positive root:

y=x2+9.y = \sqrt{x^2 + 9}.

Differentiate to check: dydx=12(x2+9)−1/2⋅2x=xx2+9=xy\dfrac{\mathrm{d}y}{\mathrm{d}x} = \dfrac{1}{2}(x^2 + 9)^{-1/2}\cdot 2x = \dfrac{x}{\sqrt{x^2 + 9}} = \dfrac{x}{y}, which matches the original equation, and y(0)=9=3y(0) = \sqrt{9} = 3, confirming the initial condition.

Interpreting solutions in modelling

Many first-order equations model change. For example, exponential growth or decay arises from dNdt=kN\frac{\mathrm{d}N}{\mathrm{d}t} = kN, whose separable solution is N=N0ektN = N_0 e^{kt}; the constant k>0k > 0 gives growth and k<0k < 0 gives decay. Newton's law of cooling, dTdt=−k(T−Ts)\frac{\mathrm{d}T}{\mathrm{d}t} = -k(T - T_s), models a body cooling towards the surrounding temperature TsT_s. When interpreting, state the units, the long-term behaviour (limiting value as t→∞t \to \infty where relevant), and whether the solution's domain is restricted by the physical situation.

Common exam traps
Forgetting the constant before applying the condition
You must include +c+ c in the general solution, then use the initial condition to find it. Omitting cc makes a particular solution impossible.
Splitting variables incorrectly
Separation only works when the right side factors as g(x)h(y)g(x)h(y). If xx and yy are tangled in a sum (such as dydx=x+y\frac{\mathrm{d}y}{\mathrm{d}x} = x + y), separation fails and a different method is needed.
Choosing the wrong root
After integrating you may get y2=…y^2 = \dots; the sign of the square root must match the given condition or stated domain.
Mishandling the constant in log solutions
When ∫1y dy=ln⁡∣y∣\int \frac{1}{y}\,\mathrm{d}y = \ln|y| appears, exponentiate carefully and absorb constants; ln⁡∣y∣=x+c\ln|y| = x + c becomes y=Aexy = A e^{x} with A=±ecA = \pm e^{c}.

Examples in context

Example 1. Direct integration. dydx=6x2\frac{\mathrm{d}y}{\mathrm{d}x} = 6x^2 with y(0)=1y(0) = 1 gives y=2x3+cy = 2x^3 + c, and c=1c = 1, so y=2x3+1y = 2x^3 + 1.

Example 2. Exponential decay. dNdt=−0.2N\frac{\mathrm{d}N}{\mathrm{d}t} = -0.2N with N(0)=50N(0) = 50 separates to ln⁡N=−0.2t+c\ln N = -0.2t + c, giving N=50e−0.2tN = 50 e^{-0.2t}; as t→∞t \to \infty, N→0N \to 0.

Try this

Q1. Solve dydx=cos⁡x\frac{\mathrm{d}y}{\mathrm{d}x} = \cos x with y(0)=2y(0) = 2. [2 marks]

  • Cue. y=sin⁡x+cy = \sin x + c; c=2c = 2, so y=sin⁡x+2y = \sin x + 2.

Q2. Solve dydx=2xy\frac{\mathrm{d}y}{\mathrm{d}x} = 2xy with y(0)=1y(0) = 1, y>0y > 0. [3 marks]

  • Cue. Separate: 1ydy=2x dx\frac{1}{y}\mathrm{d}y = 2x\,\mathrm{d}x, so ln⁡y=x2+c\ln y = x^2 + c; c=0c = 0, giving y=ex2y = e^{x^2}.

Q3. A tank's water height satisfies dhdt=−3\frac{\mathrm{d}h}{\mathrm{d}t} = -3 with h(0)=12h(0) = 12. When is the tank empty? [2 marks]

  • Cue. h=12−3th = 12 - 3t; empty when h=0h = 0, so t=4t = 4.

Exam-style questions

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

2025 VCAA-style1 mark
A tank initially contains 5 kg of salt dissolved in 3000 litres of water. Salty water containing 0.1 kg of salt per litre enters at 20 litres per minute. The solution is kept thoroughly mixed and drains at the same rate of 20 litres per minute. Let Q be the quantity of salt (kg) in the tank at time t minutes. Show that Q satisfies the differential equation dQ/dt = (300 - Q)/150.
Show worked answer →

Build the rate equation: dQ/dt = (rate salt enters) - (rate salt leaves).

Rate in: the inflow is 20 L/min at 0.1 kg/L, so salt enters at 0.1 x 20 = 2 kg/min.

Rate out: the volume stays at 3000 L (in rate equals out rate), so the concentration is Q/3000 kg/L. The outflow is 20 L/min, so salt leaves at (Q/3000) x 20 = Q/150 kg/min.

Therefore dQ/dt = 2 - Q/150. Writing 2 as 300/150 and combining over the common denominator gives dQ/dt = (300 - Q)/150, as required.

2025 VCAA-style3 marks
Use calculus to solve the differential equation dQ/dt = (300 - Q)/150, expressing Q in terms of t. (The tank initially contains 5 kg of salt.)
Show worked answer →

Separate the variables: 1/(300 - Q) dQ = 1/150 dt.

Integrate both sides: the left gives -log_e|300 - Q| and the right gives t/150 + c, so -log_e(300 - Q) = t/150 + c (the quantity 300 - Q stays positive here).

Apply the initial condition Q = 5 when t = 0: -log_e(295) = c.

Substitute back: -log_e(300 - Q) = t/150 - log_e(295), so log_e(300 - Q) = log_e(295) - t/150.

Exponentiate: 300 - Q = 295 e^(-t/150), giving Q = 300 - 295 e^(-t/150).

2023 VCAA-style1 mark
A fish farmer releases 200 fish into a pond that originally contained no fish. The fish population P grows according to the logistic model dP/dt = P(1 - P/1000), where t is in years. One form of the solution for P is P = 1000/(1 + D e^(-t)), where D is a real constant. Find the value of D.
Show worked answer →

Use the initial condition: the farmer releases 200 fish at t = 0, so P = 200 when t = 0.

Substitute into P = 1000/(1 + D e^(-t)): 200 = 1000/(1 + D e^0) = 1000/(1 + D).

Solve for D: 1 + D = 1000/200 = 5, so D = 4.

Therefore D = 4, and the particular solution is P = 1000/(1 + 4 e^(-t)).

Practise this

Sources & how we know this

ExamExplained