Differential Equations: cram pack

Formula sheets, the key method steps for every topic and 8 practice problems, with full solutions at the back.

New problemsBack to the course

In the print window, choose “Save as PDF” to keep a copy on your device.

1Formulas and methods

Differential equations

Exponential growth/decay
dy/dt = ky → y = y0ekt
Half-life
k = ln 2 / t½
Newton's law of cooling
T(t) = Ta + (T0 − Ta)e−kt
Integrating factor
y′ + p(x)y = q(x): μ = e∫p dx, (μy)′ = μq
Mixing tank
dm/dt = r·c − (r/V)·m

Integrals

Power rule
∫xⁿ dx = xⁿ⁺¹/(n + 1) + C
Definite integral
∫ab f dx = F(b) − F(a)
Common integrals
∫eˣ dx = eˣ, ∫1/x dx = ln|x|, ∫sin x dx = −cos x
Area between curves
∫ab (f − g) dx

Method

  1. Rewrite the integrand as a sum of simple terms where possible.
  2. Apply the power rule: ∫xⁿ dx = xⁿ⁺¹/(n + 1), and standard integrals for ex, 1/x, sin and cos.
  3. Use substitution when you see a function and its derivative together.
  4. For a definite integral, evaluate F(b) − F(a).

Derivatives

Power rule
(xⁿ)′ = n·xⁿ⁻¹
Product rule
(uv)′ = u′v + uv′
Quotient rule
(u/v)′ = (u′v − uv′)/v²
Chain rule
f(g(x))′ = f′(g(x))·g′(x)
Common derivatives
(eˣ)′ = eˣ, (ln x)′ = 1/x, (sin x)′ = cos x, (cos x)′ = −sin x

Method

  1. Identify the structure: sum, product, quotient or composition.
  2. Apply the power rule to each power of x.
  3. Use the product, quotient or chain rule where needed.
  4. Simplify, and evaluate at the point if asked.

Numerical methods

Bisection
c = (a + b)/2, width after n steps (b − a)/2ⁿ
Newton–Raphson
xn+1 = xₙ − f(xₙ)/f′(xₙ)
Linear interpolation
y = y₁ + (y₂ − y₁)(x − x₁)/(x₂ − x₁)
Trapezoid rule
T = (h/2)[f₀ + 2f₁ + … + 2fn−1 + fₙ], h = (b − a)/n
Simpson’s rule (n even)
S = (h/3)[f₀ + 4f₁ + 2f₂ + … + 4fn−1 + fₙ]
Trapezoid error bound
|E| ≤ (b − a)h²M/12, M = max|f″|
Euler
yn+1 = yₙ + h·f(tₙ, yₙ)
Heun
yn+1 = yₙ + (h/2)(k₁ + k₂), k₂ = f(tₙ + h, yₙ + hk₁)
RK4
yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄)

Method

  1. For a root: check the sign change for bisection, or compute f and f′ for Newton, and keep six or more digits in a table of iterates.
  2. For a table value: interpolate linearly between the two neighbouring points.
  3. For an integral: find h = (b − a)/n, list f at every node, apply the weights (1, 2, …, 2, 1 for trapezoid; 1, 4, 2, …, 4, 1 for Simpson with even n) and multiply by h/2 or h/3.
  4. For an ODE: step from t₀ with yₙ₊₁ = yₙ + h·f(tₙ, yₙ) (Euler) or the averaged slope (Heun), one row per step, until you reach the requested t.

2Practice problems

1.Differential equationsMedium
tyy(0) = 4y(1.5) = ?

Solve y″ + 8y′ + 15y = 0 with y(0) = 4 and y′(0) = 6. Find y(1.5).

y″ + 8y′ + 15y = 0y(0) = 4y′(0) = 6
y(1.5) =
2.Differential equationsHard
t (min)m (kg)80 minm = ?

A tank holds 1000 L of water in which 20 kg of salt is dissolved. Brine with 0.25 kg of salt per litre flows in at 20 L/min, and the well-mixed solution flows out at the same rate. How much salt m is in the tank after 80 min?

V = 1000 Lm0 = 20 kgc = 0.25 kg/Lr = 20 L/min
m = kg
3.IntegralsMedium
xy2y = 2x

The region under y = 2x from x = 0 to x = 2 is rotated about the x-axis. Find the volume V of the solid (a cone).

y = 2x0 ≤ x ≤ 2
V =
4.IntegralsHard
xyy = 2x2​y = 3x

Find the area A of the region enclosed by y = 2x2 and y = 3x.

y = 2x2y = 3x
A =
5.DerivativesMedium
xyx = −1

Let f(x) = (4x + 1)4. Find f′(−1).

f(x) = (4x + 1)4x = −1
f′(−1) =
6.DerivativesHard
xy

Let f(x) = x3 − 3x2 − 9x + 2. Find the x-coordinate x of the local maximum.

f(x) = x3 − 3x2 − 9x + 2
x =
7.Numerical methodsMedium
t (s)051015202530v (m/s)491413192221

A vehicle's speed is measured every 5 s (table). Estimate the distance s travelled from t = 0 to t = 30 s with the trapezoid rule.

h = 5 sn = 6
s = m
8.Numerical methodsHard
00.5112ty(0, 0.5)

Solve y′ = t + y, y(0) = 0.5 with Euler's method and step h = 0.25 up to t = 1. The exact solution is y = −t − 1 + (y₀ + 1)eᵗ. Find the global error |E| = |y(1) − y4|.

y′ = t + yy(0) = 0.5h = 0.25
|E| =

3Solutions

1.Differential equations

Answer: y(1.5) = 0.1394

Second-order linear equations: for y″ + b·y′ + c·y = 0 try y = ert: this gives the characteristic equation r² + b·r + c = 0. Two real roots r₁ ≠ r₂ give y = A·er₁t + B·er₂t (overdamped); a double root r gives y = (A + Bt)·ert (critically damped); complex roots −α ± jω give y = e−αt(A cos ωt + B sin ωt) (underdamped, oscillating). The initial values y(0) and y′(0) fix A and B. Written as y″ + 2ζωny′ + ωn²y = 0, ωn is the natural frequency and ζ the damping ratio: ζ < 1 oscillates. With a constant input, y″ + by′ + cy = c·K, the solution is the steady value K plus the free response.

Characteristic roots

(r + 3)(r + 5) = 0 → r = −3, −5

Constants from the initial values

A + B = 4, −3A − 5B = 6
A = 13, B = −9

Evaluate

y(1.5) = 13·e−4.5 + (−9)·e−7.5 = 0.1394

Answer

y(1.5) = 0.1394
2.Differential equations

Answer: m = 203.6 kg

First-order differential equations: dy/dt = k·y has the solution y = y₀·ekt (growth for k > 0, decay for k < 0); a half-life T gives k = ln 2/T. Newton's law of cooling, dT/dt = −k·(T − Tₐ), gives T(t) = Tₐ + (T₀ − Tₐ)·e−kt. A linear equation y′ + p(x)·y = q(x) is solved with the integrating factor μ = e∫p dx: (μy)′ = μq, so μy = ∫μq dx + C.

Set up the equation

dm/dt = 20·0.25 − (20/1000)·m = 5 − 0.02·m

Solve

The integrating factor is ert/V; the solution approaches the equilibrium c·V.

m(t) = 250 + (−230)·e−t/50

Substitute t

m(80) = 250 + (−230)·0.2019 = 203.6 kg

Answer

m = 203.6 kg
3.Integrals

Answer: V = 33.51

Volumes of revolution: rotating the region under y = f(x) from x = a to x = b about the x-axis makes a solid whose cross-section at x is a disc of radius f(x). Adding the discs gives V = π∫ₐᵇ f(x)² dx. When the region lies between an outer curve R(x) and an inner curve r(x), each slice is a washer: V = π∫ₐᵇ (R(x)² − r(x)²) dx. Square the radii before you subtract them.

Disc method

V = π∫₀2 4x² dx

Integrate

V = π·4·2³/3 = 33.51

Answer

V = 33.51

Check: a cone of radius 4 and height 2 has (1/3)πr²h = 33.51.

4.Integrals

Answer: A = 1.125

Definite integrals: ∫ₐᵇ f(x) dx is the signed area between the graph and the x-axis from a to b. Find an antiderivative F (with F′ = f) and evaluate F(b) − F(a). Power rule: ∫xⁿ dx = xⁿ⁺¹/(n + 1). The area between two curves is ∫(upper − lower) dx.

Intersection points

2x2 = 3x → x = 0, x = 1.5

Integrate upper minus lower

A = [1.5x2 − 0.6667x3]01.5 = 1.125

Answer

A = 1.125
5.Derivatives

Answer: f′(−1) = −432

Differentiation: the derivative f′(x) is the slope of the graph at x. Power rule: (xⁿ)′ = n·xⁿ⁻¹. Product rule: (uv)′ = u′v + uv′. Chain rule: (f(g(x)))′ = f′(g(x))·g′(x). Where f′(x) = 0 the graph has a horizontal tangent: a maximum, minimum or saddle point.

Chain rule

f′(x) = 4·(4x + 1)3·4

Substitute

f′(−1) = 4·(−3)3·4 = −432

Answer

f′(−1) = −432
6.Derivatives

Answer: x = −1

Differentiation: the derivative f′(x) is the slope of the graph at x. Power rule: (xⁿ)′ = n·xⁿ⁻¹. Product rule: (uv)′ = u′v + uv′. Chain rule: (f(g(x)))′ = f′(g(x))·g′(x). Where f′(x) = 0 the graph has a horizontal tangent: a maximum, minimum or saddle point.

Derivative

f′(x) = 3x2 − 6x − 9 = 3·(x2 − 2x − 3)

Solve f′(x) = 0

x = −1 or x = 3

Which one is which

f″(x) = 6x − 6: negative at a maximum, positive at a minimum.

f″(−1) = −12, f″(3) = 12
x = −1

Answer

x = −1
7.Numerical methods

Answer: s = 447.5 m

Numerical integration: split [a, b] into n strips of width h = (b − a)/n with nodes x₀ … xₙ. Trapezoid rule: T = (h/2)[f₀ + 2f₁ + 2f₂ + … + 2fn−1 + fₙ]. Simpson's rule (n even): S = (h/3)[f₀ + 4f₁ + 2f₂ + 4f₃ + … + 4fn−1 + fₙ]. The trapezoid error shrinks like h² (bound (b − a)h²·M/12 with M = max|f″|), Simpson's like h⁴, so Simpson is usually far more accurate for the same nodes.

Weighted sum

4 + 2·9 + 2·14 + 2·13 + 2·19 + 2·22 + 21 = 179

Trapezoid rule

s = (5/2)·179 = 447.5 m

Answer

s = 447.5 m
8.Numerical methods

Answer: |E| = 0.4153

Numerical solution of y′ = f(t, y): Euler's method steps along the slope at the start of each step: yn+1 = yₙ + h·f(tₙ, yₙ). Heun's method (improved Euler) averages the slope at the start and at an Euler prediction: k₁ = f(tₙ, yₙ), k₂ = f(tₙ + h, yₙ + h·k₁), yn+1 = yₙ + (h/2)(k₁ + k₂). The classical Runge–Kutta method (RK4) combines four slopes: yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄). The global error of Euler is proportional to h, Heun's to h², RK4's to h⁴.

ntₙyₙf(tₙ, yₙ)yₙ₊₁000.50.50.62510.250.6250.8750.8437520.50.843751.343751.17968830.751.1796881.9296881.662109

Euler steps

y4 = 1.662109

Exact value

y(1) = 2.077423

Global error

|E| = |2.077423 − 1.662109| = 0.4153
00.5112ty

Graph

Answer

|E| = 0.4153

Halving h roughly halves the error of Euler's method.

ExamProblems · examproblems.com/courses/differential-equations