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)38Ta​T = ?

A forged steel part at 480 °C is left to cool in a workshop at 20 °C. It follows Newton's law of cooling, dT/dt = −k·(T − Tₐ). With k = 0.05 min⁻¹, find its temperature T after 38 min.

T0 = 480 °CTa = 20 °Ck = 0.05 min⁻¹t = 38 min
T = °C
2.Differential equationsHard
tyK = 2y(0.75) = ?

A system obeys y″ + 3y′ + 27.25y = 54.5 for t ≥ 0, starting from rest: y(0) = 0, y′(0) = 0. Find y(0.75); angles in radians.

y″ + 3y′ + 27.25y = 54.5y(0) = 0y′(0) = 0
y(0.75) =
3.Matrices & determinantsMedium
A =62−3−5B =44−4−6

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

A = [6 2; −3 −5]B = [4 4; −4 −6]
c21 =
4.Matrices & determinantsHard
A =−54−2−4b =−1112

Solve the system −5x + 4y = −11 and −2x − 4y = 12. Find x.

−5x + 4y = −11−2x − 4y = 12
x =
5.Vectors: dot & cross productMedium
xyzdFAB

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

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

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

a = (1, 2, 0)b = (−3, 3, 2)c = (2, 2, −3)
V =
7.IntegralsMedium
xy02

Evaluate I = ∫02 e−0.5x dx.

f(x) = e−0.5x[0, 2]
I =
8.IntegralsHard
xyy = 2x2​y = x

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

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

For f(x, y) = −4x³y² + 5x³ − 4y² − 3xy, 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) = −4x³y² + 5x³ − 4y² − 3xyP = (2, −2)
θ = °
10.Multivariable calculusHard
xyx = 4y = 0.5x42

Evaluate I = ∬T (2x − 2y) dA over the triangle T bounded by y = 0, x = 4 and y = 0.5x.

f = 2x − 2y0 ≤ x ≤ 40 ≤ y ≤ 0.5x
I =
11.Sequences & seriesMedium
AΣ 10ᵏ/k!k = 1 … ∞BΣ ln k / kk = 1 … ∞CΣ (1 + 1/k)ᵏk = 1 … ∞DΣ 1/√kk = 1 … ∞

Exactly one of the four series converges. Click it.

A: Σ 10ᵏ/k!B: Σ ln k / kC: Σ (1 + 1/k)ᵏD: Σ 1/√k
  1. 1A: Σ 10ᵏ/k!
  2. 2B: Σ ln k / k
  3. 3C: Σ (1 + 1/k)ᵏ
  4. 4D: Σ 1/√k
12.Sequences & seriesHard
xyfP₄1.5

Use the Maclaurin polynomial of degree 4 of f(x) = e2x to approximate f(1.5). Give P₄(1.5).

f(x) = e2xn = 4x = 1.5
P₄(1.5) =
13.Laplace transformsMedium
t (s)f(t)t = 0.2f(0.2) = ?

Find f(t) for F(s) = (4s + 2)/(s² + 6s + 9) and evaluate f(0.2).

F(s) = (4s + 2)/(s² + 6s + 9)t = 0.2 s
f(0.2) =
14.Laplace transformsHard
t (s)y(t)t = 0.25y(0.25) = ?

A system obeys y″ + 10y′ + 24y = 72·u(t). It starts at rest (y(0) = 0, y′(0) = 0) and the input is a step u(t) = 10 for t ≥ 0. Use the Laplace transform to find the output y(0.25).

y″ + 10y′ + 24y = 72·u(t)u(t) = 10y(0) = 0, y′(0) = 0t = 0.25 s
y(0.25) =
15.Convolution & LTI systemsMedium
x(t)t10h(t)t0.51

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

h(t) = e−2tu(t)x(t) = 10·u(t)t = 0.5 s
y(0.5) =
16.Convolution & LTI systemsHard
x(t)t (ms)0.25122

A pulse train has height 2, pulse width 0.25 ms and period 1 ms. Find the amplitude c1 of harmonic number 1, at 1000 Hz.

A = 2τ = 0.25 msT = 1 ms
c1 =

3Solutions

1.Differential equations

Answer: T = 88.8 °C

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.

Solution of the cooling law

T(t) = 20 + 460·e−0.05t

Substitute t

T(38) = 20 + 460·0.1496 = 88.8 °C

Answer

T = 88.8 °C
2.Differential equations

Answer: y(0.75) = 2.644

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.

Steady value and roots

yss = 2, r = −1.5 ± 5j

Constants from the initial values

y(0) = 2 + A = 0 → A = −2
y′(0) = −1.5A + 5B = 0 → B = −0.6

Evaluate

y(0.75) = 2·(1 − e−1.125(cos 3.75 + 0.3·sin 3.75)) = 2.644

Answer

y(0.75) = 2.644
3.Matrices & determinants

Answer: c21 = 8

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

c21 = (−3)·4 + (−5)·(−4) = 8

Answer

c21 = 8
4.Matrices & determinants

Answer: x = −0.1429

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 = (−5)·(−4) − 4·(−2) = 28

Replace the x column

det Aₓ = (−11)·(−4) − 4·12 = −4
x = −4/28 = −0.1429

Answer

x = −0.1429
y = −2.929
5.Vectors: dot & cross product

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

Work

W = F·d = (−15)·(−6) + 30·3 + 0·(−2) = 180 J

Answer

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

Answer: V = 23

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 = 3·(−3) − 2·2 = −13
(b × c)y = bzcx − bxcz = 2·2 − (−3)·(−3) = −5
(b × c)z = bxcy − bycx = (−3)·2 − 3·2 = −12

Dot product with a

a·(b × c) = 1·(−13) + 2·(−5) + 0·(−12) = −23
V = |−23| = 23

Answer

V = 23
7.Integrals

Answer: I = 1.264

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) = e−0.5x/(−0.5)

Evaluate at the limits

I = (e−1 − 1)/(−0.5) = (0.3679 − 1)/(−0.5) = 1.264

Answer

I = 1.264
8.Integrals

Answer: A = 0.04167

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

Integrate upper minus lower

A = [0.5x2 − 0.6667x3]00.5 = 0.04167

Answer

A = 0.04167
9.Multivariable calculus

Answer: θ = 132.4 °

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 = −12x²y² + 15x² − 3y → fx(2, −2) = −126
fy = −8x³y − 8y − 3x → fy(2, −2) = 138

Direction

θ = atan2(138, −126) = 132.4°

Answer

θ = 132.4°
10.Multivariable calculus

Answer: I = 16

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

Limits

0 ≤ x ≤ 4, 0 ≤ y ≤ 0.5x

Inner integral

∫00.5x (2x − 2y) dy = 2x·0.5x − 2·(0.5x)²/2 = 0.75x²

Outer integral

I = 0.75·4³/3 = 16

Answer

I = 16
11.Sequences & series

Answer: A: Σ 10ᵏ/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: Σ 10ᵏ/k! → Ratio test: L = lim 10/(k + 1) = 0 < 1, converges.
B: Σ ln k / k → Larger than 1/k for k ≥ 3: diverges by comparison with the harmonic series.
C: Σ (1 + 1/k)ᵏ → The terms tend to e, not 0: diverges by the n-th term test.
D: Σ 1/√k → p-series with p = 1/2 ≤ 1: diverges.

The answer

A: Σ 10ᵏ/k!

Answer

A: Σ 10ᵏ/k!
12.Sequences & series

Answer: P₄(1.5) = 16.38

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.

Maclaurin polynomial

P₄(x) = 1 + 2x + 2x² + 1.333x³ + 0.6667x⁴

Substitute

P₄(1.5) = 1 + 3 + 4.5 + 4.5 + 3.375 = 16.38

Answer

P₄(1.5) = 16.38

The exact value is f(1.5) = 20.085537; the error 3.710537 shrinks when you take more terms or a point closer to 0.

13.Laplace transforms

Answer: f(0.2) = 1.098

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² + 6s + 9 = (s + 3)²

Partial fractions

4s + 2 = 4(s + 3) + −10
F(s) = 4/(s + 3) − 10/(s + 3)²

Inverse transform

f(t) = 4·e−3t − 10·t·e−3t
f(0.2) = (4 − 10·0.2)·e−3·0.2 = 1.098

Answer

f(0.2) = 1.098
14.Laplace transforms

Answer: y(0.25) = 10.28

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) = 720/(s(s + 4)(s + 6))

Partial fractions

A = 720/(4·6) = 30
B = 720/((−4)·2) = −90
C = 720/((−6)·(−2)) = 60

Output

y(t) = 30 − 90·e−4t + 60·e−6t
y(0.25) = 10.28

Answer

y(0.25) = 10.28
15.Convolution & LTI systems

Answer: y(0.5) = 3.161

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

Evaluate

y(0.5) = 5·(1 − e−1) = 3.161

Answer

y(0.5) = 3.161
16.Convolution & LTI systems

Answer: c1 = 0.9003

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

Duty cycle

D = 0.25/1 = 0.25

Harmonic amplitude

c1 = (2·2/(1π))·|sin(1π·0.25)| = 0.9003

Answer

c1 = 0.9003

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