Engineering Mathematics 2: cram pack

Formula sheets, the key method steps for every topic and 16 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

Matrices & determinants

2×2 determinant
det = ad − bc
2×2 inverse
A⁻¹ = (1/det)·[d −b; −c a]
Product
(AB)ij = Σk AikBkj
Eigenvalues
det(A − λI) = 0
Cramer’s rule
xi = det(Ai)/det(A)

Vectors: dot & cross product

Magnitude
|a| = √(ax² + ay² + az²)
Unit vector
û = a/|a|
Components from an angle
Fx = F·cos θ, Fy = F·sin θ
Dot product
a·b = axbx + ayby + azbz = |a||b|cos θ
Scalar projection of a onto b
a·b/|b|
Work
W = F·d
Cross product
a × b = (aybz − azby, azbx − axbz, axby − aybx)
Areas
parallelogram |a × b|, triangle ½|a × b|
Moment about O
MO = r × F, r = A − O
Volume of a parallelepiped
V = |a·(b × c)|

Method

  1. Write every vector in components; for a force given by a magnitude and an angle, use F·cos θ and F·sin θ.
  2. For lengths and directions use the magnitude √(x² + y² + z²) and the unit vector a/|a|.
  3. Use the dot product for angles, projections, perpendicularity and work: a·b = |a||b|cos θ.
  4. Use the cross product for perpendicular vectors, areas and moments: write the three components carefully, keeping the order a × b.

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).

Multivariable calculus

Partial derivatives
fx = ∂f/∂x (y constant), fxy = fyx
Gradient
∇f = (fx, fy), largest slope |∇f|
Directional derivative
Duf = ∇f·u, u = v/|v|
Tangent plane
z = f(a, b) + fx(x − a) + fy(y − b)
Chain rule
dz/dt = fx·dx/dt + fy·dy/dt
Second-derivative test
D = fxxfyy − fxy²: D > 0, fxx > 0 minimum; D > 0, fxx < 0 maximum; D < 0 saddle
Double integral
∬R f dA = ∫ab∫g(x)h(x) f dy dx
Polar coordinates
x = r cos θ, y = r sin θ, dA = r dr dθ

Method

  1. For a partial derivative, treat the other variable as a constant and use the one-variable rules.
  2. For gradient questions compute ∇f = (fx, fy) at the point; use |∇f| for the steepest slope and ∇f·u with a unit vector u for a directional derivative.
  3. For extrema solve fx = 0 and fy = 0 together, then classify every point with D = fxx·fyy − fxy².
  4. For a double integral sketch the region, choose the order (or polar coordinates for discs) so the inner limits are the bounding curves, and integrate inside first.

Sequences & series

Arithmetic sequence
aₙ = a₁ + (n − 1)d, Sₙ = n(a₁ + aₙ)/2
Geometric sequence
aₙ = a₁rⁿ⁻¹, Sₙ = a₁(1 − rⁿ)/(1 − r)
Infinite geometric series
S = a₁/(1 − r), |r| < 1
p-series
Σ1/kᵖ converges only for p > 1
Ratio test
L = lim |ak+1/ak|: L < 1 converges, L > 1 diverges
Radius of convergence
R = lim |ck/ck+1|
Taylor polynomial
Pₙ(x) = Σ f⁽ᵏ⁾(a)(x − a)ᵏ/k!
Maclaurin series
eˣ = Σxᵏ/k!, sin x = x − x³/3! + …, cos x = 1 − x²/2! + …, ln(1 + x) = x − x²/2 + …
Lagrange error bound
|Rₙ| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!

Method

  1. Identify the type: arithmetic (constant difference), geometric (constant ratio), a standard series, a power series or a Taylor polynomial.
  2. For sums use Sₙ = n(a₁ + aₙ)/2, Sₙ = a₁(1 − rⁿ)/(1 − r), or S = a₁/(1 − r) when |r| < 1, with a₁ the first term of the sum.
  3. For convergence: n-th term test first, then compare with Σ1/kᵖ or a geometric series, and use the ratio or root test for factorials and k-th powers.
  4. For Taylor polynomials list f(a), f′(a), f″(a), …, divide the k-th derivative by k!, and bound the error with M·|x − a|ⁿ⁺¹/(n + 1)!.

Convolution & LTI systems

Discrete convolution
y[n] = Σ x[k]·h[n − k]
Continuous convolution
y(t) = ∫ x(τ)·h(t − τ) dτ
Length and sum
Ly = Lx + Lh − 1, Σy = Σx·Σh
Start, end and area
t1 + t3 … t2 + t4, ∫y = ∫x·∫h
Step into e^(−at)u(t)
y(t) = (A/a)(1 − e−at)
First-order system
y′ + a·y = b·x → h(t) = b·e−atu(t)
First-order recursion
y[n] = a·y[n − 1] + b·x[n] → h[n] = b·aⁿ
Sinusoid through h(t)
|H(jω)|, ∠H(jω); H(jω) = b/(a + jω)
Fourier series amplitudes
square ±A: 4A/(kπ), odd k; triangle ±A: 8A/(π²k²), odd k; pulse train: A·D and (2A/(kπ))|sin kπD|

Method

  1. Sketch both signals and note where each one is nonzero.
  2. Flip h and slide it: for each output time, find the range where both signals overlap.
  3. Multiply and add (discrete) or multiply and integrate (continuous) over that overlap only.
  4. Check the result: the output starts at the sum of the start times, ends at the sum of the end times, and its area or sum is the product of the input areas or sums.

2Practice problems

1.Differential equationsMedium
tyy(0) = 3y(2) = ?

Solve y″ + 4y′ + 4y = 0 with y(0) = 3 and y′(0) = 0. Find y(2).

y″ + 4y′ + 4y = 0y(0) = 3y′(0) = 0
y(2) =
2.Differential equationsHard
ty2y(2) = ?

Solve y′ + 2y = 2e−3t with y(0) = 2. Find y(2).

y′ + 2y = 2e−3ty(0) = 2
y(2) =
3.Matrices & determinantsMedium
A =−1−2−3325−104

Find the determinant det A of the 3×3 matrix A shown.

A = [−1 −2 −3; 3 2 5; −1 0 4]
det A =
4.Matrices & determinantsHard
A =1−435b =−1−7

Solve the system x − 4y = −1 and 3x + 5y = −7. Find x.

x − 4y = −13x + 5y = −7
x =
5.Vectors: dot & cross productMedium
xyzdFAB

A constant force F = (−10, 5, −5) N acts on an object while it moves in a straight line from point A(−2, 2, −2) to point B(−5, 2, −1). Coordinates are in metres. Find the work W done by the force.

F = (−10, 5, −5) NA(−2, 2, −2) mB(−5, 2, −1) m
W = J
6.Vectors: dot & cross productHard
xyzrFOA

A force F = (20, −20, 0) N acts at point A(3, 0, 1). Coordinates are in metres. Find the z-component Mz of the moment of the force about the origin O(0, 0, 0).

F = (20, −20, 0) NA(3, 0, 1) mO(0, 0, 0) m
Mz = N·m
7.IntegralsMedium
xy01

Evaluate I = ∫01 e0.5x dx.

f(x) = e0.5x[0, 1]
I =
8.IntegralsHard
xyy = 2x2​y = 3x

Find the area A of the region enclosed by y = 2x2 and y = 3x.

y = 2x2y = 3x
A =
9.Multivariable calculusMedium
xy−2−2∇fP

For f(x, y) = xy² − 4y² − 2xy, find the direction of steepest ascent at (−2, −2) as the angle θ from the positive x-axis, in degrees between −180° and 180°.

f(x, y) = xy² − 4y² − 2xyP = (−2, −2)
θ = °
10.Multivariable calculusHard
xyy = 3xy = x²39

Find the area I = ∬R 1 dA of the region R between the parabola y = x² and the line y = 3x.

y = x²y = 3xf = 1
I =
11.Sequences & seriesMedium
AΣ (−1)ᵏ/√kk = 1 … ∞BΣ 3/kk = 1 … ∞CΣ (7/4)ᵏk = 1 … ∞DΣ (1 + 1/k)ᵏk = 1 … ∞

Exactly one of the four series converges. Click it.

A: Σ (−1)ᵏ/√kB: Σ 3/kC: Σ (7/4)ᵏD: Σ (1 + 1/k)ᵏ
  1. 1A: Σ (−1)ᵏ/√k
  2. 2B: Σ 3/k
  3. 3C: Σ (7/4)ᵏ
  4. 4D: Σ (1 + 1/k)ᵏ
12.Sequences & seriesHard

f(x) = sin x is approximated at x = 1.5 by its Maclaurin polynomial of degree n. With the Lagrange bound M·|x|ⁿ⁺¹/(n + 1)! and M = 1, find the smallest degree n that guarantees an error below 0.0005.

f(x) = sin xx = 1.5M = 1ε = 0.0005
n =
13.Laplace transformsMedium
t (s)f(t)t = 0.1f(0.1) = ?

Find f(t) for F(s) = (s + 17)/(s² + 10s + 25) and evaluate f(0.1).

F(s) = (s + 17)/(s² + 10s + 25)t = 0.1 s
f(0.1) =
14.Laplace transformsHard
t (s)y(t)t = 0.1y(0.1) = ?

A system obeys y″ + 9y′ + 18y = 36·u(t), with y(0) = 4 and y′(0) = 0. At t = 0 a step input u(t) = 5 is applied. Use the Laplace transform to find y(0.1).

y″ + 9y′ + 18y = 36·u(t)y(0) = 4y′(0) = 0u(t) = 5t = 0.1 s
y(0.1) =
15.Convolution & LTI systemsMedium
x(t)t4h(t)t0.251

A system has the impulse response h(t) = e−4tu(t). The input is a step x(t) = 4·u(t). Find the output y(0.05) at t = 0.05 s.

h(t) = e−4tu(t)x(t) = 4·u(t)t = 0.05 s
y(0.05) =
16.Convolution & LTI systemsHard
δ[n]0.5+h[n]−0.8·z⁻¹+

A system obeys y[n] = −0.8·y[n − 1] + 0.5·x[n] and starts at rest. Find its impulse response sample h[5].

y[n] = −0.8y[n − 1] + 0.5x[n]x[n] = δ[n]
h[5] =

3Solutions

1.Differential equations

Answer: y(2) = 0.2747

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.

Double root

(r + 2)² = 0 → r = −2

Constants

A = y(0) = 3
y′(0) = B − 2A → B = 0 + 6 = 6

Evaluate

y(2) = (3 + 6·2)·e−4 = 0.2747

Answer

y(2) = 0.2747
2.Differential equations

Answer: y(2) = 0.06831

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

μ = e2t
(e2ty)′ = 2e−t

General solution

y = 2/(−1)·e−3t + C·e−2t = −2·e−3t + C·e−2t

Initial condition

2 = −2 + C → C = 4

Evaluate

y(2) = −2·0.002479 + 4·0.01832 = 0.06831

Answer

y(2) = 0.06831
3.Matrices & determinants

Answer: det A = 20

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.

2×2 minors of the first row

M₁₁ = 8, M₁₂ = 17, M₁₃ = 2

Cofactor expansion

det A = (−1)·8 − (−2)·17 + (−3)·2 = 20

Answer

det A = 20
4.Matrices & determinants

Answer: x = −1.941

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.

Determinant of the coefficient matrix

det A = 1·5 − (−4)·3 = 17

Replace the x column

det Aₓ = (−1)·5 − (−4)·(−7) = −33
x = −33/17 = −1.941

Answer

x = −1.941
y = −0.2353
5.Vectors: dot & cross product

Answer: W = 25 J

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.

Displacement

d = B − A = (−5, 2, −1) − (−2, 2, −2) = (−3, 0, 1) m

Work

W = F·d = (−10)·(−3) + 5·0 + (−5)·1 = 25 J

Answer

W = 25 J
6.Vectors: dot & cross product

Answer: Mz = −60 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.

Position vector from O to A

r = A − O = (3, 0, 1) − (0, 0, 0) = (3, 0, 1) m

Moment M = r × F

(r × F)x = ryFz − rzFy = 0·0 − 1·(−20) = 20 N·m
(r × F)y = rzFx − rxFz = 1·20 − 3·0 = 20 N·m
(r × F)z = rxFy − ryFx = 3·(−20) − 0·20 = −60 N·m

Requested component

Mz = −60 N·m

Answer

Mz = −60 N·m
7.Integrals

Answer: I = 1.297

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.

Antiderivative

F(x) = e0.5x/0.5

Evaluate at the limits

I = (e0.5 − 1)/0.5 = (1.649 − 1)/0.5 = 1.297

Answer

I = 1.297
8.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
9.Multivariable calculus

Answer: θ = 74.05 °

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.

Partial derivatives

fx = y² − 2y → fx(−2, −2) = 8
fy = 2xy − 8y − 2x → fy(−2, −2) = 28

Direction

θ = atan2(28, 8) = 74.05°

Answer

θ = 74.05°
10.Multivariable calculus

Answer: I = 4.5

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).

Intersections

x² = 3x → x = 0, x = 3

Inner integral

∫x²3x 1 dy = 3x − x²

Outer integral

I = ∫03 (3x − x²) dx = 4.5

Answer

I = 4.5
11.Sequences & series

Answer: A: Σ (−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.

Test each series

A: Σ (−1)ᵏ/√k → Alternating, with 1/√k decreasing to 0: converges (conditionally).
B: Σ 3/k → A multiple of the harmonic series (p = 1): diverges.
C: Σ (7/4)ᵏ → Geometric with r = 7/4 > 1: the terms grow, so it diverges.
D: Σ (1 + 1/k)ᵏ → The terms tend to e, not 0: diverges by the n-th term test.

The answer

A: Σ (−1)ᵏ/√k

Answer

A: Σ (−1)ᵏ/√k
12.Sequences & series

Answer: n = 8

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.

n(n + 1)!bound011.5121.125260.56253240.210941200.0632857200.01582650400.003397403200.000635683628800.0001059

Bounds for increasing n

n = 7: 0.0006356 ≥ 0.0005
n = 8: 0.0001059 < 0.0005

Smallest degree

n = 8

Answer

n = 8
13.Laplace transforms

Answer: f(0.1) = 1.334

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).

Factor the denominator

s² + 10s + 25 = (s + 5)²

Partial fractions

s + 17 = 1(s + 5) + 12
F(s) = 1/(s + 5) + 12/(s + 5)²

Inverse transform

f(t) = e−5t + 12·t·e−5t
f(0.1) = (1 + 12·0.1)·e−5·0.1 = 1.334

Answer

f(0.1) = 1.334
14.Laplace transforms

Answer: y(0.1) = 4.403

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).

Transform the equation

(s² + 9s + 18)·Y(s) = 4s + 36 + 180/s
Y(s) = (4s² + 36s + 180)/(s(s + 3)(s + 6))

Partial fractions

A = 180/(3·6) = 10
B = 108/((−3)·3) = −12
C = 108/((−6)·(−3)) = 6

Result

y(t) = 10 − 12·e−3t + 6·e−6t
y(0.1) = 4.403

Answer

y(0.1) = 4.403

Check: y(0) = A + B + C = 4.

15.Convolution & LTI systems

Answer: y(0.05) = 0.1813

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).

Convolution with a step

y(t) = 4·∫₀ᵗ e−4τ dτ = (4/4)·(1 − e−4t)

Evaluate

y(0.05) = 1·(1 − e−0.2) = 0.1813

Answer

y(0.05) = 0.1813
16.Convolution & LTI systems

Answer: h[5] = −0.1638

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)|.

Iterate

h[0] = 0.5
h[1] = −0.4
h[2] = 0.32
h[3] = −0.256
h[4] = 0.2048
h[5] = −0.1638

Result

h[5] = 0.5·(−0.8)5 = −0.1638

Answer

h[5] = −0.1638

ExamProblems · examproblems.com/courses/engineering-mathematics-2