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
Formula sheets, the key method steps for every topic and 16 practice problems, with full solutions at the back.
A cup of coffee at 90 °C is left in a room at 20 °C. After 15 min its temperature is 58 °C. It follows Newton's law of cooling, dT/dt = −k·(T − Tₐ). Find the constant k in min⁻¹.
A system obeys y″ + 2y′ + 10y = 20 for t ≥ 0, starting from rest: y(0) = 0, y′(0) = 0. Find y(0.25); angles in radians.
Find the determinant det A of the 3×3 matrix A shown.
Find the larger eigenvalue λ of the matrix A shown.
Find the angle θ between the vectors a = (1, −4, −1) and b = (5, 3, 6), in degrees.
For a = (−2, 1, 1) and b = (−6, −5, −5), find the y-component (a × b)y of the cross product a × b.
Evaluate I = ∫02 2x·(x² + 1)1 dx.
Find the area A of the region enclosed by y = 2x2 and y = 2x.
z = 3 sin(−x + 3y) with x = −t + 2 and y = −2t² + 2. Use the chain rule to find dz/dt at t = −1.
Use polar coordinates to evaluate I = ∬R e−(x² + y²) dA, where R is the quarter of the disc x² + y² ≤ 2.25 in the first quadrant.
Find the sum S = Σk=119 4/(k(k + 1)).
Recognise the series and find its sum: S = Σk=0∞ (−1)ᵏ·1.22k/(2k)!.
Find f(t) for F(s) = 15/(s(s + 3)(s + 6)) by partial fractions, and evaluate f(0.25).
A system is described by the differential equation y″ + 4y′ + 29y = 7u, with input u(t) and output y(t). Find its transfer function G(s) = Y(s)/U(s) (zero initial conditions) and evaluate the impulse response g(t) at t = 0.3 s. Find g(0.3).
The signal x(t) is nonzero only between t = 0 s and t = 2 s, and h(t) only between t = 2 s and t = 6 s. Find the duration D of y(t) = x(t) ∗ h(t).
A system obeys y[n] = −0.8·y[n − 1] + 0.5·x[n] and starts at rest. The input is the unit step x[n] = u[n]. Find y[2].
Answer: k = 0.04073 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.
Answer: y(0.25) = 0.5064
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.
Answer: det A = −100
Matrices: the determinant of a 2×2 matrix is ad − bc. For 3×3, expand along a row: each entry times the determinant of the 2×2 matrix left after removing its row and column, with alternating signs + − +. Eigenvalues λ solve det(A − λI) = 0; for 2×2 that is λ² − (trace)·λ + det = 0.
Answer: λ = −1
Matrices: the determinant of a 2×2 matrix is ad − bc. For 3×3, expand along a row: each entry times the determinant of the 2×2 matrix left after removing its row and column, with alternating signs + − +. Eigenvalues λ solve det(A − λI) = 0; for 2×2 that is λ² − (trace)·λ + det = 0.
Answer: θ = 111.5 °
Dot product: a·b = axbx + ayby + azbz = |a|·|b|·cos θ, a scalar. It gives the angle between two vectors, cos θ = a·b/(|a|·|b|), and it is zero when they are perpendicular. The scalar projection of a onto b is a·b/|b|, and a constant force F moving through a displacement d does work W = F·d.
Answer: (a × b)y = −16
Cross product: a × b = (aybz − azby, azbx − axbz, axby − aybx) is perpendicular to both vectors, and |a × b| = |a|·|b|·sin θ is the area of the parallelogram they span (the triangle has half of it). Order matters: b × a = −(a × b). The moment of a force F applied at point A about point O is MO = r × F, with r = A − O.
Check: a × b is perpendicular to both vectors, so its dot product with a and with b is zero.
Answer: I = 12
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.
Answer: A = 0.3333
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.
Answer: dz/dt = −38.61
Gradient: ∇f = (fx, fy) points in the direction of steepest ascent, and |∇f| is the largest rate of change. The directional derivative along a direction v is Duf = ∇f·u with the unit vector u = v/|v|. The tangent plane at (a, b) is z = f(a, b) + fx(a, b)(x − a) + fy(a, b)(y − b); it gives the linear approximation of f near that point. Chain rule: dz/dt = fx·dx/dt + fy·dy/dt.
Answer: I = 0.7026
Critical points and double integrals: a critical point has fx = fy = 0. Second-derivative test with D = fxxfyy − fxy²: D > 0 and fxx > 0 is a local minimum, D > 0 and fxx < 0 a local maximum, D < 0 a saddle point, and D = 0 decides nothing. For a double integral over a general region, let the inner limits be the curves that bound the region; for discs and sectors use polar coordinates, x = r cos θ, y = r sin θ and dA = r dr dθ (do not forget the factor r).
Answer: S = 3.8
Convergence: if the terms do not tend to 0 the series diverges (n-th term test). Σ1/kᵖ converges only for p > 1, and Σrᵏ only for |r| < 1. Compare a series with one of these when its terms behave like them for large k. Ratio test: L = lim |ak+1/ak|; root test: L = lim |ak|1/k. L < 1 converges, L > 1 diverges, L = 1 decides nothing. A power series Σck(x − x₀)ᵏ converges for |x − x₀| < R with R = lim |ck/ck+1|; an alternating series with terms decreasing to 0 converges.
Answer: S = 0.3624
Taylor polynomials: near x = a a smooth function is approximated by Pₙ(x) = f(a) + f′(a)(x − a) + f″(a)(x − a)²/2! + … + f⁽ⁿ⁾(a)(x − a)ⁿ/n!; for a = 0 this is the Maclaurin polynomial. Standard series: eˣ = Σxᵏ/k!, sin x = x − x³/3! + x⁵/5! − …, cos x = 1 − x²/2! + x⁴/4! − …, ln(1 + x) = x − x²/2 + x³/3 − …, 1/(1 − x) = Σxᵏ. Lagrange error bound: |f(x) − Pₙ(x)| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!, where M is the largest value of |f⁽ⁿ⁺¹⁾| between a and x.
Answer: f(0.25) = 0.232
Laplace transform: F(s) = ∫₀∞ f(t)·e−st dt turns derivatives into multiplications by s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Standard pairs: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), and a shift s → s + σ multiplies the time function by e−σt. To invert a rational F(s), split it into partial fractions: the residue of a simple pole p is A = [(s − p)·F(s)] at s = p. Final value theorem: f(∞) = lims→0 s·F(s), valid only when all poles of s·F(s) lie in the left half-plane. Initial value theorem: f(0⁺) = lims→∞ s·F(s).
Answer: g(0.3) = 0.7664
Laplace transform: F(s) = ∫₀∞ f(t)·e−st dt turns derivatives into multiplications by s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Standard pairs: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), and a shift s → s + σ multiplies the time function by e−σt. To invert a rational F(s), split it into partial fractions: the residue of a simple pole p is A = [(s − p)·F(s)] at s = p. Final value theorem: f(∞) = lims→0 s·F(s), valid only when all poles of s·F(s) lie in the left half-plane. Initial value theorem: f(0⁺) = lims→∞ s·F(s).
Answer: D = 6 s
Continuous convolution: y(t) = ∫ x(τ)·h(t − τ) dτ. Flip h, slide it to position t and integrate the product of the overlapping parts. For causal signals the integral runs from 0 to t. Useful results: the area of y is the area of x times the area of h; if x lasts from t₁ to t₂ and h from t₃ to t₄, y lasts from t₁ + t₃ to t₂ + t₄. For h(t) = e−atu(t) and a step of height A, y(t) = (A/a)(1 − e−at).
Answer: y[2] = 0.42
Impulse response and Fourier series: the impulse response h(t) describes an LTI system completely: for y′ + a·y = b·x it is h(t) = b·e−atu(t), and for y[n] = a·y[n − 1] + b·x[n] it is h[n] = b·aⁿu[n]. A sinusoid cos ωt passes through with gain |H(jω)| and phase ∠H(jω), where H is the Fourier transform of h. A periodic signal is a sum of harmonics at k/T: a square wave of ±A has only odd harmonics with amplitudes 4A/(kπ), a triangle wave of peak A has 8A/(π²k²), and a pulse train of height A and duty cycle D has the mean A·D and the harmonics (2A/(kπ))·|sin(kπD)|.
ExamProblems · examproblems.com/courses/engineering-mathematics-2