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 forged steel part at 300 °C is left to cool in a workshop at 30 °C. After 26 min its temperature is 104 °C. It follows Newton's law of cooling, dT/dt = −k·(T − Tₐ). How long after the start, t, does it take to reach 60 °C?
Solve the initial value problem y′ + 2y/x = 6x2 for x > 0, with y(1) = 2. Find y(3).
Find the determinant det A of the 3×3 matrix A shown.
Solve the system 3x − 5y = −4 and −x + 2y = −9. Find x.
Find the value of k for which the vectors a = (1, 3, 0) and b = (−4, k, 4) are perpendicular.
Find the area A of the triangle with corners A(1, 1, 2), B(−1, −2, 1) and C(0, 4, 3). Coordinates are in metres.
Evaluate I = ∫₁1.5 x·ln x dx.
Find the area A of the region enclosed by y = 0.5x2 and y = x.
The tangent plane to z = 2xe0.5y at the point where (x, y) = (−1, 1) can be written as z = px + qy + c. Find c.
The integral I = ∫00.8 ∫y0.8 ex² dx dy cannot be done in this order. Reverse the order of integration and evaluate it.
Find the radius of convergence R of the power series Σk=1∞ 3ᵏ·xᵏ/(k + 1).
Use the Maclaurin polynomial of degree 2 of f(x) = eˣ to approximate f(1.2). Give P₂(1.2).
Find f(t) for F(s) = 45/(s(s + 5)(s + 6)) by partial fractions, and evaluate f(0.3).
A system obeys y″ + 7y′ + 6y = 3·u(t). It starts at rest (y(0) = 0, y′(0) = 0) and the input is a step u(t) = 2 for t ≥ 0. Use the Laplace transform to find the output y(0.3).
A system has the impulse response h(t) = e−1tu(t). The input is a step x(t) = 4·u(t). Find the output y(0.4) at t = 0.4 s.
A square wave switches between +10 and −10 with period T = 5 ms (half a period at each level). Find the amplitude c3 of its Fourier component at 600 Hz.
Answer: t = 44.14 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(3) = 32.49
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: det A = −4
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 = −53
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: k = 1.333
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 = 4.743 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 = 0.1436
Integration by parts: ∫u·dv = u·v − ∫v·du, from the product rule. Choose u as the factor that gets simpler when you differentiate it (a power of x, or ln x), and dv as the factor you can integrate (eˣ, sin x, cos x, or x when u = ln x). For a definite integral, evaluate u·v at both limits: ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. A factor x² needs two rounds.
Answer: A = 0.6667
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: c = 1.649
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.4482
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: R = 0.3333
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.
The ratio test says nothing at the two end points; they have to be checked separately.
Answer: P₂(1.2) = 2.92
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.
The exact value is f(1.2) = 3.320117; the error 0.400117 shrinks when you take more terms or a point closer to 0.
Answer: f(0.3) = 0.7316
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(0.3) = 0.1441
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(0.4) = 1.319
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: c3 = 4.244
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