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
t (min)T (°C)34Ta​350 °C389 °C

A forged steel part at 480 °C is left to cool in a workshop at 30 °C. After 34 min its temperature is 350 °C. It follows Newton's law of cooling, dT/dt = −k·(T − Tₐ). How long after the start, t, does it take to reach 389 °C?

T0 = 480 °CTa = 30 °CT(34) = 350 °CT = 389 °C
t = min
2.Differential equationsHard
xy13(1, 4)y(3) = ?

Solve the initial value problem y′ + 3y/x = x for x > 0, with y(1) = 4. Find y(3).

y′ + 3y/x = xy(1) = 4
y(3) =
3.IntegralsMedium
xy02

Evaluate I = ∫02 2x·(x² + 1)3 dx.

f(x) = 2x·(x² + 1)3[0, 2]
I =
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 = −2

Let f(x) = ln(x2 + 5). Find f′(−2).

f(x) = ln(x2 + 5)x = −2
f′(−2) =
6.DerivativesHard
xyx = 3

The tangent line to f(x) = 3x2 + 1 at x = 3 crosses the y-axis at y = b. Find b.

f(x) = 3x2 + 1x = 3
b =
7.Numerical methodsMedium
xy00.66671.3332

Approximate ∫02 √(1 + x³) dx with the trapezoid rule using n = 3 strips. Give T.

f(x) = √(1 + x³)[a, b] = [0, 2]n = 3
T =
8.Numerical methodsHard
00.511.50.40.60.8ty(0, 0.5)

Use Heun's method (improved Euler) with h = 0.5 for y′ = t − y, y(0) = 0.5 to estimate y3 ≈ y(1.5).

y′ = t − yy(0) = 0.5h = 0.5
y3 =

3Solutions

1.Differential equations

Answer: t = 22.53 min

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.

Cooling constant

k = ln(450/320)/34 = 0.01003 min⁻¹

Time to reach the target

t = ln(450/359)/0.01003 = 22.53 min

Answer

t = 22.53 min
2.Differential equations

Answer: y(3) = 1.941

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.

Integrating factor

μ = e∫3/x dx = x3

Integrate

(x3y)′ = x4 → x3y = 1x5/5 + C

Initial condition

4 = 1/5 + C → C = 3.8

Evaluate

y(3) = (1·35/5 + 3.8)/33 = 1.941

Answer

y(3) = 1.941
3.Integrals

Answer: I = 156

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.

Substitution

u = x² + 1, du = 2x dx
I = ∫15 u3 du

Integrate and evaluate

I = [u4/4]15 = (625 − 1)/4 = 156

Answer

I = 156
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′(−2) = −0.4444

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 for the logarithm

f′(x) = 2x/(x2 + 5)

Substitute

f′(−2) = −4/9 = −0.4444

Answer

f′(−2) = −0.4444
6.Derivatives

Answer: b = −26

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.

Point and slope

f(3) = 28
f′(x) = 6x → f′(3) = 18

Tangent line at x = 0

b = f(x₀) − f′(x₀)·x₀ = 28 − 18·3 = −26

Answer

b = −26
7.Numerical methods

Answer: T = 3.316

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.

ixᵢf(xᵢ)weight001110.6666671.13855221.3333331.83585723231

Nodes and values

h = 0.6667
Σ wᵢfᵢ = 1 + 2·1.13855 + 2·1.835857 + 3 = 9.948814

Trapezoid rule

T = (0.6667/2)·9.948814 = 3.316

Answer

T = 3.316

This integral has no elementary antiderivative, which is exactly when a numerical rule is needed.

8.Numerical methods

Answer: y3 = 0.8662

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ₙk₁k₂yₙ₊₁000.5−0.50.250.437510.50.43750.06250.531250.585938210.5859380.4140630.7070310.866211

Heun steps

y₁ = 0.5 + (0.5/2)(−0.5 + 0.25) = 0.4375
y₂ = 0.4375 + (0.5/2)(0.0625 + 0.53125) = 0.585938
y₃ = 0.585938 + (0.5/2)(0.414063 + 0.707031) = 0.866211
00.511.50.40.60.8ty

Graph

y3 = 0.8662

Answer

y3 = 0.8662

The exact value is y(1.5) = 0.834695.

ExamProblems · examproblems.com/courses/differential-equations