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 can of soda at 4 °C is taken out of the fridge into a room at 30 °C. After 9 min its temperature is 17 °C. It follows Newton's law of cooling, dT/dt = −k·(T − Tₐ). Find the constant k in min⁻¹.
Solve y′ + 3y = 8e−0.5t with y(0) = 3. Find y(0.5).
Let C = AB for the matrices shown. Find the entry c21 (row 2, column 1).
Find the larger eigenvalue λ of the matrix A shown.
Find the dot product a·b of a = (7, −7, 0) and b = (2, 5, 1).
The unit vector n̂ is perpendicular to both a = (−3, 3, −2) and b = (1, −5, 1) and has a positive z-component. Find its y-component n̂y.
Evaluate I = ∫₀1 x·sin x dx (x in radians).
Find the area A of the region enclosed by y = 0.5x2 and y = 2x.
For f(x, y) = −xy − 4y³, find the directional derivative Duf at (−2, 1) in the direction of v = (4, −1).
Find the critical points of f(x, y) = −x³ + 12x − y² + 4y, classify them, and give the value fmax of f at its local maximum.
Apply the ratio test to Σk=1∞ k!/kᵏ. Find L = lim |ak+1/ak|.
f(x) = ln(1 + x) is approximated by its Maclaurin polynomial of degree 3 at x = 0.5. Use the Lagrange error bound M·|x|ⁿ⁺¹/(n + 1)! with M = 3! = 6 to find the largest possible error |R₃|.
Find f(t) for F(s) = (s + 12)/(s² + 6s + 10) by completing the square, and evaluate f(0.6).
A system is described by the differential equation y″ + 6y′ + 18y = 9u, 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.25 s. Find g(0.25).
The pulse x(t) = 3 for 0 ≤ t < 2 s (zero elsewhere) is applied to a system with the impulse response h(t) = 3 for 0 ≤ t < 3 s (zero elsewhere). Find the output y(2) at t = 2 s.
A system obeys y[n] = 0.5·y[n − 1] + 2·x[n] and starts at rest. The input is the unit step x[n] = u[n]. Find y[2].
Answer: k = 0.07702 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.5) = 2.448
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: c21 = −6
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: λ = 5
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 = −21
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: n̂y = 0.0718
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.
The z-component of a × b is already positive.
Answer: I = 0.3012
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 = 5.333
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: Duf = 1.455
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: fmax = 20
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).
Solve fx = 0 and fy = 0 together.
Answer: L = 0.3679
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.
L < 1, so the series converges.
Answer: |R₃| = 0.01563
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 true error is smaller: |f(0.5) − P₃(0.5)| = 0.0112.
Answer: f(0.6) = 0.9764
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.25) = 0.966
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) = 18
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 18 and falls back to 0 at t = 5 s.
Answer: y[2] = 3.5
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