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
t (min)T (°C)31Ta​118 °C64 °C

A forged steel part at 340 °C is left to cool in a workshop at 28 °C. After 31 min its temperature is 118 °C. It follows Newton's law of cooling, dT/dt = −k·(T − Tₐ). How long after the start, t, does it take to reach 64 °C?

T0 = 340 °CTa = 28 °CT(31) = 118 °CT = 64 °C
t = min
2.Differential equationsHard
t (min)m (kg)t = ?16 kg

A tank holds 400 L of pure water. Brine with 0.1 kg of salt per litre flows in at 4 L/min, and the well-mixed solution flows out at the same rate. After how many minutes t does the tank contain 16 kg of salt?

V = 400 Lm0 = 0 kgc = 0.1 kg/Lr = 4 L/min
t = min
3.Matrices & determinantsMedium
A =−4−64−5B =32−44

Let C = AB for the matrices shown. Find the entry c12 (row 1, column 2).

A = [−4 −6; 4 −5]B = [3 2; −4 4]
c12 =
4.Matrices & determinantsHard
A =−213−6b =−75

Solve the system −2x + y = −7 and 3x − 6y = 5. Find x.

−2x + y = −73x − 6y = 5
x =
5.Vectors: dot & cross productMedium
xyzdFAB

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

F = (20, −15, −20) NA(−3, −3, 5) mB(−1, 2, −1) m
W = J
6.Vectors: dot & cross productHard
xyzabc

Find the volume V of the parallelepiped whose edges from one corner are a = (0, 4, −1), b = (5, −5, −5) and c = (3, 3, −5).

a = (0, 4, −1)b = (5, −5, −5)c = (3, 3, −5)
V =
7.IntegralsMedium
xyπ/3y = x·cos x

Evaluate I = ∫₀π/3 x·cos x dx (x in radians).

∫₀π/3 x·cos x dx
I =
8.IntegralsHard
xy4y = 4xy = x²

The region between the line y = 4x and the parabola y = x² is rotated about the x-axis. Find the volume V of the solid.

y = 4xy = x²
V =
9.Multivariable calculusMedium
xy−202−2−10PQ

For f(x, y) = 3 sin(3x − y), find the directional derivative Duf at P(0.5, −1) in the direction towards Q(−2.5, −2).

f(x, y) = 3 sin(3x − y)P = (0.5, −1)Q = (−2.5, −2)
Duf =
10.Multivariable calculusHard
xyy = 2y = x22

The integral I = ∫02 ∫x2 sin(y²) dy dx cannot be done in this order. Reverse the order of integration and evaluate it.

f = sin(y²)0 ≤ x ≤ y ≤ 2
I =
11.Sequences & seriesMedium
AΣ 3ᵏ/k³k = 1 … ∞BΣ (4/3)ᵏk = 1 … ∞CΣ ln k / kk = 1 … ∞DΣ (−1)ᵏ/√kk = 1 … ∞

Exactly one of the four series converges. Click it.

A: Σ 3ᵏ/k³B: Σ (4/3)ᵏC: Σ ln k / kD: Σ (−1)ᵏ/√k
  1. 1A: Σ 3ᵏ/k³
  2. 2B: Σ (4/3)ᵏ
  3. 3C: Σ ln k / k
  4. 4D: Σ (−1)ᵏ/√k
12.Sequences & seriesHard

Recognise the series and find its sum: S = Σk=1∞ (−1)k+1·0.2ᵏ/k.

ak = (−1)k+1·0.2ᵏ/k
S =
13.Laplace transformsMedium
ReIm−5−105

Find f(t) for F(s) = (3s + 33)/(s² + 4s + 29) by completing the square, and evaluate f(0.2).

F(s) = (3s + 33)/(s² + 4s + 29)t = 0.2 s
f(0.2) =
14.Laplace transformsHard
t (s)y(t)t = 1y(1) = ?

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

y″ + 5y′ + 4y = 12·u(t)y(0) = 6y′(0) = 0u(t) = 2t = 1 s
y(1) =
15.Convolution & LTI systemsMedium
x(t)t1h(t)t0.21

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

h(t) = e−5tu(t)x(t) = 1·u(t)t = 0.2 s
y(0.2) =
16.Convolution & LTI systemsHard
δ(t)y′ + 5y = 5xh(t)

A system obeys y′(t) + 5·y(t) = 5·x(t). Find its impulse response at t = 1 s, h(1).

y′ + 5y = 5xt = 1 s
h(1) =

3Solutions

1.Differential equations

Answer: t = 53.85 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.

Cooling constant

k = ln(312/90)/31 = 0.0401 min⁻¹

Time to reach the target

t = ln(312/36)/0.0401 = 53.85 min

Answer

t = 53.85 min
2.Differential equations

Answer: t = 51.08 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.

Result

m(t) = 40 + (−40)·e−t/100

Solve for t

e−t/100 = (40 − 16)/(40 − 0) = 0.6
t = 100·ln(1.667) = 51.08 min

Answer

t = 51.08 min
3.Matrices & determinants

Answer: c12 = −32

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.

Row times column

c12 = (−4)·2 + (−6)·4 = −32

Answer

c12 = −32
4.Matrices & determinants

Answer: x = 4.111

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 = (−2)·(−6) − 1·3 = 9

Replace the x column

det Aₓ = (−7)·(−6) − 1·5 = 37
x = 37/9 = 4.111

Answer

x = 4.111
y = 1.222
5.Vectors: dot & cross product

Answer: W = 85 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 = (−1, 2, −1) − (−3, −3, 5) = (2, 5, −6) m

Work

W = F·d = 20·2 + (−15)·5 + (−20)·(−6) = 85 J

Answer

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

Answer: V = 10

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.

Cross product b × c

(b × c)x = bycz − bzcy = (−5)·(−5) − (−5)·3 = 40
(b × c)y = bzcx − bxcz = (−5)·3 − 5·(−5) = 10
(b × c)z = bxcy − bycx = 5·3 − (−5)·3 = 30

Dot product with a

a·(b × c) = 0·40 + 4·10 + (−1)·30 = 10
V = |10| = 10

Answer

V = 10
7.Integrals

Answer: I = 0.4069

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.

Choose u and dv

u = x, du = dx; dv = cos x dx, v = sin x

Antiderivative

∫x cos x dx = x sin x − ∫sin x dx = x sin x + cos x

Limits

(π/3·sin π/3 + cos π/3) − (0 + 1) = 0.4069

Answer

I = 0.4069
8.Integrals

Answer: V = 428.9

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.

Intersection points

4x = 1x² → x = 0, x = 4

Washer method

V = π∫₀4 ((4x)² − (1x²)²) dx

Integrate

V = π(16·4³/3 − 1·4⁵/5) = π·136.5 = 428.9

Answer

V = 428.9
9.Multivariable calculus

Answer: Duf = 6.08

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.

Gradient at the point

fx = 9 cos(3x − y) → fx(0.5, −1) = −7.21
fy = −3 cos(3x − y) → fy(0.5, −1) = 2.403

Unit vector

|v| = √((−3)² + (−1)²) = 3.162
u = (−0.9487, −0.3162)

Dot product

Duf = (−7.21)·(−0.9487) + 2.403·(−0.3162) = 6.08

Answer

Duf = 6.08
10.Multivariable calculus

Answer: I = 0.8268

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

Region

0 ≤ x ≤ y, 0 ≤ y ≤ 2

Reversed order

I = ∫02 ∫0y sin(y²) dx dy = ∫02 y·sin(y²) dy

Substitution u = y²

I = [−cos(y²)/2]02 = 0.8268

Answer

I = 0.8268
11.Sequences & series

Answer: D: Σ (−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: Σ 3ᵏ/k³ → Ratio test: L = 3 > 1, diverges.
B: Σ (4/3)ᵏ → Geometric with r = 4/3 > 1: the terms grow, so it diverges.
C: Σ ln k / k → Larger than 1/k for k ≥ 3: diverges by comparison with the harmonic series.
D: Σ (−1)ᵏ/√k → Alternating, with 1/√k decreasing to 0: converges (conditionally).

The answer

D: Σ (−1)ᵏ/√k

Answer

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

Answer: S = 0.1823

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.

Standard series

Σ (−1)k+1xᵏ/k = ln(1 + x)

Substitute

x = 0.2
S = ln 1.2 = 0.1823

Answer

S = 0.1823
13.Laplace transforms

Answer: f(0.2) = 4.132

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

Complete the square

s² + 4s + 29 = (s + 2)² + 25 = (s + 2)² + 5²

Match the numerator

F(s) = 3·(s + 2)/((s + 2)² + 5²) + (27/5)·5/((s + 2)² + 5²)

Inverse transform

f(t) = e−2t·(3·cos 5t + 5.4·sin 5t)
f(0.2) = e−2·0.2·(3·cos(5·0.2) + 5.4·sin(5·0.2)) = 4.132

Answer

f(0.2) = 4.132
14.Laplace transforms

Answer: y(1) = 6

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² + 5s + 4)·Y(s) = 6s + 30 + 24/s
Y(s) = (6s² + 30s + 24)/(s(s + 1)(s + 4))

Partial fractions

A = 24/(1·4) = 6
B = 0/((−1)·3) = 0
C = 0/((−4)·(−3)) = 0

Result

y(t) = 6
y(1) = 6

Answer

y(1) = 6

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

15.Convolution & LTI systems

Answer: y(0.2) = 0.1264

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) = 1·∫₀ᵗ e−5τ dτ = (1/5)·(1 − e−5t)

Evaluate

y(0.2) = 0.2·(1 − e−1) = 0.1264

Answer

y(0.2) = 0.1264
16.Convolution & LTI systems

Answer: h(1) = 0.03369

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

Transfer function

H(s) = Y(s)/X(s) = 5/(s + 5)

Impulse response

h(t) = 5·e−5tu(t)
h(1) = 5·e−5 = 0.03369

Answer

h(1) = 0.03369

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