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.
Solve y″ + 7y′ + 12y = 0 with y(0) = 1 and y′(0) = 6. Find y(0.5).
A system obeys y″ + 3y′ + 18.25y = 18.25 for t ≥ 0, starting from rest: y(0) = 0, y′(0) = 0. Find y(0.5); angles in radians.
Find the determinant det A of the 3×3 matrix A shown.
Solve the system −3x − 3y = −9 and 5x − 4y = 8. Find x.
Find the dot product a·b of a = (−5, 0, 3) and b = (6, −6, −2).
A force F = (0, 30, 30) N acts at point A(−2, 0, 3). Coordinates are in metres. Find the x-component Mx of the moment of the force about the origin O(0, 0, 0).
Evaluate I = ∫01 2x·(x² + 1)3 dx.
The region between the line y = x and the parabola y = x² is rotated about the x-axis. Find the volume V of the solid.
For f(x, y) = 2xe3y, find the maximum rate of change of f at the point (2, −1), that is |∇f|.
Evaluate I = ∬T (3x − y) dA over the triangle T bounded by y = 0, x = 4 and y = 2x.
Exactly one of the four series converges. Click it.
Find the coefficient c₃ of (x − 1)³ in the Taylor series of f(x) = e−x about a = 1.
Find f(t) for F(s) = 30/(s(s + 2)(s + 4)) by partial fractions, and evaluate f(0.6).
A system is described by the differential equation y″ + 9y′ + 18y = 14u, 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.5 s. Find g(0.5).
The pulse x(t) = 2 for 0 ≤ t < 3 s (zero elsewhere) is applied to a system with the impulse response h(t) = 5 for 0 ≤ t < 4 s (zero elsewhere). Find the output y(2.5) at t = 2.5 s.
The input x(t) = 5·cos(4t) is applied to a system with h(t) = 5·e−2tu(t). Find the amplitude Ym of the steady-state output.
Answer: y(0.5) = 1.013
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: y(0.5) = 1.036
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 = −35
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: x = 2.222
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: a·b = −36
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.
Negative: the angle between the vectors is more than 90°.
Answer: Mx = −90 N·m
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.
Answer: I = 3.75
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: V = 0.4189
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.
Answer: |∇f| = 0.6057
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 = 85.33
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: C: Σ 1/k²
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: c₃ = −0.06131
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.6) = 1.831
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.5) = 0.8089
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: y(2.5) = 25
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).
h(t − τ) is the pulse of h flipped and shifted so that it ends at τ = t.
y(t) is a trapezoid: it rises from 0, stays at 30 and falls back to 0 at t = 7 s.
Answer: Ym = 5.59
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