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)42Ta​232 °C

A forged steel part at 500 °C is left to cool in a workshop at 29 °C. After 42 min its temperature is 232 °C. It follows Newton's law of cooling, dT/dt = −k·(T − Tₐ). Find the constant k in min⁻¹.

T0 = 500 °CTa = 29 °CT(42) = 232 °C
k = min⁻¹
2.Differential equationsHard
t (min)m (kg)120 minm = ?

A tank holds 400 L of pure water. Brine with 0.05 kg of salt per litre flows in at 2 L/min, and the well-mixed solution flows out at the same rate. How much salt m is in the tank after 120 min?

V = 400 Lm0 = 0 kgc = 0.05 kg/Lr = 2 L/min
m = kg
3.Matrices & determinantsMedium
A =−2−5−43−23−25−1

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

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

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

3x − 6y = 12−5x − 6y = −7
x =
5.Vectors: dot & cross productMedium
xyzdFAB

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

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

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

F = (−10, −20, 10) NA(0, 3, 3) mO(0, 0, 0) m
Mx = N·m
7.IntegralsMedium
xyπ/2y = x·cos x

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

∫₀π/2 x·cos x dx
I =
8.IntegralsHard
xyy = 2x2​y = 4x

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

y = 2x2y = 4x
A =
9.Multivariable calculusMedium
xy−1−1P

Use the tangent plane of f(x, y) = −2x³e2y at (−1, −1) to approximate f(−1.1, −0.9). Give the linear approximation L.

f(x, y) = −2x³e2y(a, b) = (−1, −1)(x, y) = (−1.1, −0.9)
L =
10.Multivariable calculusHard
xy−20−2−101(0, 0)(−1, −1)

Find the critical points of f(x, y) = x³ + y³ + 3xy, classify them, and give the value fmax of f at its local maximum.

f(x, y) = x³ + y³ + 3xy
fmax =
11.Sequences & seriesMedium
xx₀ = 0??R

Find the radius of convergence R of the power series Σk=1∞ 16ᵏx²ᵏ/k.

Σ 16ᵏx²ᵏ/kx₀ = 0
R =
12.Sequences & seriesHard

f(x) = sin x is approximated at x = 0.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 = 0.5M = 1ε = 0.0005
n =
13.Laplace transformsMedium
t (s)f(t)t = 1.5f(1.5) = ?

Find f(t) for F(s) = 35/(s(s + 1)(s + 5)) by partial fractions, and evaluate f(1.5).

F(s) = 35/(s(s + 1)(s + 5))t = 1.5 s
f(1.5) =
14.Laplace transformsHard
t (s)y(t)t = 0.3y(0.3) = ?

A system obeys y″ + 7y′ + 12y = 6·u(t). It starts at rest (y(0) = 0, y′(0) = 0) and the input is a step u(t) = 5 for t ≥ 0. Use the Laplace transform to find the output y(0.3).

y″ + 7y′ + 12y = 6·u(t)u(t) = 5y(0) = 0, y′(0) = 0t = 0.3 s
y(0.3) =
15.Convolution & LTI systemsMedium
x(t)t5h(t)t0.51

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

h(t) = e−2tu(t)x(t) = 5·u(t)t = 0.2 s
y(0.2) =
16.Convolution & LTI systemsHard
5cos(20t)h(t) = 10e−10t​u(t)y(t)

The input x(t) = 5·cos(20t) is applied to a system with h(t) = 10·e−10tu(t). Find the amplitude Ym of the steady-state output.

h(t) = 10e−10tu(t)x(t) = 5cos(20t)
Ym =

3Solutions

1.Differential equations

Answer: k = 0.02004 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.

Ratio of temperature differences

e−42k = 203/471

Solve for k

k = ln(471/203)/42 = 0.02004 min⁻¹

Answer

k = 0.02004 min⁻¹
2.Differential equations

Answer: m = 9.024 kg

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.

Set up the equation

dm/dt = 2·0.05 − (2/400)·m = 0.1 − 0.005·m

Solve

The integrating factor is ert/V; the solution approaches the equilibrium c·V.

m(t) = 20 + (−20)·e−t/200

Substitute t

m(120) = 20 + (−20)·0.5488 = 9.024 kg

Answer

m = 9.024 kg
3.Matrices & determinants

Answer: det A = −3

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₁₁ = −13, M₁₂ = 3, M₁₃ = 11

Cofactor expansion

det A = (−2)·(−13) − (−5)·3 + (−4)·11 = −3

Answer

det A = −3
4.Matrices & determinants

Answer: x = 2.375

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 = 3·(−6) − (−6)·(−5) = −48

Replace the x column

det Aₓ = 12·(−6) − (−6)·(−7) = −114
x = −114/(−48) = 2.375

Answer

x = 2.375
y = −0.8125
5.Vectors: dot & cross product

Answer: W = 160 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 = (−4, −4, −5) − (0, 4, 5) = (−4, −8, −10) m

Work

W = F·d = (−20)·(−4) + (−10)·(−8) + 0·(−10) = 160 J

Answer

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

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.

Position vector from O to A

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

Moment M = r × F

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

Requested component

Mx = 90 N·m

Answer

Mx = 90 N·m
7.Integrals

Answer: I = 0.5708

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

(π/2·sin π/2 + cos π/2) − (0 + 1) = 0.5708

Answer

I = 0.5708
8.Integrals

Answer: A = 2.667

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 = 4x → x = 0, x = 2

Integrate upper minus lower

A = [2x2 − 0.6667x3]02 = 2.667

Answer

A = 2.667
9.Multivariable calculus

Answer: L = 0.406

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.

Value and slopes at the point

f(−1, −1) = 0.2707
fx = −6x²e2y → fx(−1, −1) = −0.812
fy = −4x³e2y → fy(−1, −1) = 0.5413

Tangent plane

z = 0.2707 − 0.812(x + 1) + 0.5413(y + 1)
z = −0.812x + 0.5413y + 0

Linear approximation

L = 0.2707 + (−0.812)·(−0.1) + 0.5413·0.1 = 0.406

Answer

L = 0.406

The exact value is f = 0.440026; the plane is close because the step is small.

10.Multivariable calculus

Answer: fmax = 1

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

Critical points

Solve fx = 0 and fy = 0 together.

fx = 3x² + 3y = 0
fy = 3y² + 3x = 0
(0, 0), (−1, −1)

Second-derivative test

(0, 0): fxx = 0, fyy = 0, fxy = 3, D = −9 → saddle point
(−1, −1): fxx = −6, fyy = −6, fxy = 3, D = 27 → local maximum

Value at the extremum

f(−1, −1) = 1

Answer

fmax = 1
11.Sequences & series

Answer: R = 0.25

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.

Coefficients

The series is in x²: it converges for 16·x² < 1, so |x| < 1/4.

Radius

R = 0.25

Answer

R = 0.25

The ratio test says nothing at the two end points; they have to be checked separately.

12.Sequences & series

Answer: n = 4

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)!bound010.5120.125260.020833240.00260441200.0002604

Bounds for increasing n

n = 3: 0.002604 ≥ 0.0005
n = 4: 0.0002604 < 0.0005

Smallest degree

n = 4

Answer

n = 4
13.Laplace transforms

Answer: f(1.5) = 5.049

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

Partial fractions

F(s) = A/s + B/(s + 1) + C/(s + 5)
A = 35/(1·5) = 7
B = 35/((−1)·4) = −8.75
C = 35/((−5)·(−4)) = 1.75

Inverse transform

f(t) = 7 − 8.75·e−t + 1.75·e−5t

Evaluate

f(1.5) = 7 − 8.75·e−1.5 + 1.75·e−5·1.5 = 5.049

Answer

f(1.5) = 5.049
14.Laplace transforms

Answer: y(0.3) = 0.6933

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

Y(s) = 30/(s(s + 3)(s + 4))

Partial fractions

A = 30/(3·4) = 2.5
B = 30/((−3)·1) = −10
C = 30/((−4)·(−1)) = 7.5

Output

y(t) = 2.5 − 10·e−3t + 7.5·e−4t
y(0.3) = 0.6933

Answer

y(0.3) = 0.6933
15.Convolution & LTI systems

Answer: y(0.2) = 0.8242

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

Evaluate

y(0.2) = 2.5·(1 − e−0.4) = 0.8242

Answer

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

Answer: Ym = 2.236

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

Frequency response

H(jω) = 10/(10 + jω)

Gain at the input frequency

|H(j20)| = 10/√(10² + 20²) = 0.4472
Ym = 5·0.4472 = 2.236

Answer

Ym = 2.236

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