工程数学:速成包

公式表、每个专题的关键方法步骤和 20 道练习题,完整解答附在最后。

新题目返回课程

在打印窗口中选择“另存为 PDF”,即可在设备上保存一份副本。

1公式与方法

极限

直接代入
先试 x = a;得到 0/0 说明需要化简
因式分解并约去
xⁿ − aⁿ = (x − a)(xⁿ⁻¹ + a·xⁿ⁻² + … + aⁿ⁻¹)
共轭
(√x − √a)(√x + √a) = x − a
无穷远处的极限
除以 x 的最高次幂;次数相同时极限为首项系数之比
重要极限
limx→∞ (1 + a/x)bx = eab
洛必达法则
lim f/g = lim f′/g′ (0/0 or ∞/∞)

导数

幂法则
(xⁿ)′ = n·xⁿ⁻¹
乘积法则
(uv)′ = u′v + uv′
商法则
(u/v)′ = (u′v − uv′)/v²
链式法则
f(g(x))′ = f′(g(x))·g′(x)
常用导数
(eˣ)′ = eˣ, (ln x)′ = 1/x, (sin x)′ = cos x, (cos x)′ = −sin x

方法

  1. 识别结构:和、积、商还是复合。
  2. 对 x 的每个幂应用幂法则。
  3. 需要时使用乘积法则、商法则或链式法则。
  4. 化简,如有要求再代入该点求值。

积分

幂法则
∫xⁿ dx = xⁿ⁺¹/(n + 1) + C
定积分
∫ab f dx = F(b) − F(a)
常用积分
∫eˣ dx = eˣ, ∫1/x dx = ln|x|, ∫sin x dx = −cos x
曲线间的面积
∫ab (f − g) dx

方法

  1. 尽可能把被积函数改写为简单项之和。
  2. 应用幂法则:∫xⁿ dx = xⁿ⁺¹/(n + 1),以及 ex、1/x、sin 和 cos 的基本积分。
  3. 当看到一个函数和它的导数同时出现时,用换元法。
  4. 对于定积分,计算 F(b) − F(a)。

复数

直角坐标形式
z = a + jb, j² = −1, z̄ = a − jb
模与辐角
|z| = √(a² + b²), tan θ = b/a (check the quadrant)
极坐标形式与指数形式
z = r∠θ = r(cos θ + j sin θ) = r·ejθ
除法
z₁/z₂ = z₁z̄₂/|z₂|²
乘积与商
r₁∠θ₁ · r₂∠θ₂ = r₁r₂∠(θ₁ + θ₂), r₁∠θ₁ / r₂∠θ₂ = (r₁/r₂)∠(θ₁ − θ₂)
棣莫弗定理
(r∠θ)ⁿ = rⁿ∠nθ
n 次方根
zk = R1/n∠((φ + 360°·k)/n), k = 0 … n − 1
p² < 4q 的二次方程
z = −p/2 ± j√(q − p²/4)

方法

  1. 根据运算选择合适的形式:加减用直角坐标形式 a + jb,乘除、乘方和开方用极坐标形式 r∠θ。
  2. 展开相乘时使用 j² = −1;做除法时把分子和分母同乘以分母的共轭。
  3. 用 r = √(a² + b²) 转换,θ 通过画出点 (a, b) 的草图得到(注意象限),再用 a = r cos θ,b = r sin θ 换回去。
  4. 乘方用棣莫弗定理,(r∠θ)ⁿ = rⁿ∠nθ;对于 zⁿ = R∠φ,列出全部 n 个根 R1/n∠((φ + 360°k)/n)。

微分方程

指数增长/衰减
dy/dt = ky → y = y0ekt
半衰期
k = ln 2 / t½
牛顿冷却定律
T(t) = Ta + (T0 − Ta)e−kt
积分因子
y′ + p(x)y = q(x): μ = e∫p dx, (μy)′ = μq
混合槽
dm/dt = r·c − (r/V)·m

矩阵与行列式

2×2 行列式
det = ad − bc
2×2 逆矩阵
A⁻¹ = (1/det)·[d −b; −c a]
产品
(AB)ij = Σk AikBkj
特征值
det(A − λI) = 0
克拉默法则
xi = det(Ai)/det(A)

向量:点积和叉积

模
|a| = √(ax² + ay² + az²)
单位向量
û = a/|a|
由角度求分量
Fx = F·cos θ, Fy = F·sin θ
点积
a·b = axbx + ayby + azbz = |a||b|cos θ
a 在 b 上的标量投影
a·b/|b|
功
W = F·d
叉积
a × b = (aybz − azby, azbx − axbz, axby − aybx)
面积
平行四边形 |a × b|,三角形 ½|a × b|
对 O 的力矩
MO = r × F, r = A − O
平行六面体的体积
V = |a·(b × c)|

方法

  1. 把每个向量都写成分量形式;对于由大小和角度给出的力,用 F·cos θ 和 F·sin θ。
  2. 求长度和方向时,用模 √(x² + y² + z²) 和单位向量 a/|a|。
  3. 求夹角、投影、垂直关系和功时,用点积:a·b = |a||b|cos θ。
  4. 求垂直向量、面积和力矩时,用叉积:仔细写出三个分量,并保持 a × b 的顺序。

概率

补事件
P(not A) = 1 − P(A)
加法公式
P(A ∪ B) = P(A) + P(B) − P(A ∩ B)
独立事件
P(A ∩ B) = P(A)·P(B)
条件概率
P(A|B) = P(A ∩ B)/P(B)
二项分布
P(X = k) = C(n,k)·pk(1 − p)n−k

描述统计

平均数
x̄ = Σx/n
样本标准差
s = √(Σ(x − x̄)²/(n − 1))
中位数
排序后数据的中间值
z 分数
z = (x − μ)/σ
正态分布
68% 落在 ±1σ 内,95% 落在 ±2σ 内

正态分布

z 分数
z = (x − μ)/σ
小于 x
P(X < x) = Φ(z)
大于 x
P(X > x) = 1 − Φ(z)
介于两者之间
Φ(z2) − Φ(z1)
百分位数
x = μ + z·σ

方法

  1. 把数值换算为 z 分数:z = (x − μ)/σ。
  2. 查 Φ(z),即 z 左侧的面积。
  3. 右侧面积用 1 − Φ(z);区间面积用两个左侧面积相减。
  4. 求百分位数时反向计算:先由面积求 z,再求 x = μ + z·σ。

2练习题

1.极限中等
xyy = L

求极限 L = limx→∞ (5x2 + 3x − 4)/(2x2 + 5x + 3)。

f(x) = (5x2 + 3x − 4)/(2x2 + 5x + 3)x → ∞
L =
2.极限困难
xyx = 0

求极限 L = limx→0 (ex − 1 − x)/x²。

f(x) = (ex − 1 − x)/x²x → 0
L =
3.导数中等
xyx = 2

设 f(x) = (3x + 3)3。求 f′(2)。

f(x) = (3x + 3)3x = 2
f′(2) =
4.导数困难
xy

设 f(x) = x3 − 3x2 − 9x − 1。求局部极小值点的 x 坐标 x。

f(x) = x3 − 3x2 − 9x − 1
x =
5.积分中等
xy15y = x ln x

计算 I = ∫₁5 x·ln x dx。

∫₁5 x ln x dx
I =
6.积分困难
xy2y = 4xy = 2x²

将直线 y = 4x 与抛物线 y = 2x² 之间的区域绕 x 轴旋转。求该旋转体的体积 V。

y = 4xy = 2x²
V =
7.复数中等
ReIm11−1−122−2−233−3−3θz

把 z = 3∠−75° 写成直角坐标形式 a + jb,并给出实部 a。

r = 3θ = −75°
a =
8.复数困难
ReIm22−2−244−4−4z₁z₂

z³ − 6z² + 4z + 40 = 0 的一个根为 z₁ = 4 + j2。系数为实数。求实根 r。

z₁ = 4 + j2
r =
9.微分方程中等
t (min)T (°C)24Ta​20 °C23 °C

一罐 6 °C 的汽水从冰箱中取出,放到 25 °C 的房间里。 24 min 后其温度为 20 °C。 它遵循牛顿冷却定律 dT/dt = −k·(T − Tₐ)。 从开始算起需要多长时间 t 才能达到 23 °C?

T0 = 6 °CTa = 25 °CT(24) = 20 °CT = 23 °C
t = min
10.微分方程困难
tyK = 1y(0.5) = ?

一个系统在 t ≥ 0 时满足 y″ + 2y′ + 5y = 5,从静止开始:y(0) = 0,y′(0) = 0。求 y(0.5);角度以弧度为单位。

y″ + 2y′ + 5y = 5y(0) = 0y′(0) = 0
y(0.5) =
11.矩阵与行列式中等
A =02−24255−53

求所示 3×3 矩阵 A 的行列式 det A。

A = [0 2 −2; 4 2 5; 5 −5 3]
det A =
12.矩阵与行列式困难
A =5−1−42

求所示矩阵 A 的较大特征值 λ。

A = [5 −1; −4 2]
λ =
13.向量:点积和叉积中等
xyzdFAB

一个恒力 F = (−10, −30, −30) N 作用在物体上,物体沿直线从点 A(0, −3, −5) 运动到点 B(−5, 1, 5)。坐标单位为米。求该力所做的功 W。

F = (−10, −30, −30) NA(0, −3, −5) mB(−5, 1, 5) m
W = J
14.向量:点积和叉积困难
xyzABC

求顶点为 A(4, 4, 1)、B(−1, −1, 1)、C(−4, −2, 4) 的三角形的面积 A。坐标单位为米。

A(4, 4, 1)B(−1, −1, 1)C(−4, −2, 4)
A = m²
15.概率中等
n = 6k = 3N = ?

从 6 人中选出 3 人组成委员会,有多少种选法 N?顺序无关。

n = 6k = 3
N =
16.概率困难
xp(x)10.3520.3530.140.2μ

X 的分布如图所示,Y = 10X + 5。求 Y 的标准差 σY。

x: 1, 2, 3, 4p(x): 0.35, 0.35, 0.1, 0.2Y = 10X + 5
σY =
17.描述统计中等
20253035404550

求数据的样本标准差 s:37, 34, 39, 50, 52, 20, 23。

n = 7
s =
18.描述统计困难
z −2−1012

某量服从正态分布,平均数为 50,标准差为 15。数值低于 35 的比例 P 是多少?

μ = 50σ = 15x = 35
P =
19.正态分布中等
105 s117.8 sX + Y?

X 服从均值 60 s、标准差 5 s 的正态分布;Y 独立,服从均值 45 s、标准差 4 s 的正态分布。求 X + Y 大于 117.8 s 的概率 P。z 取两位小数。

μX = 60, σX = 5μY = 45, σY = 4X + Y > 117.8 s
P =
20.正态分布困难
0.05 mm0D − d?

孔径 D 服从均值 25.05 mm、标准差 0.02 mm 的正态分布。轴径 d 独立,服从均值 25 mm、标准差 0.015 mm 的正态分布。随机选取的一根轴在 D − d > 0 时能装入。求它能装入的概率 P。z 取两位小数。

孔: μ = 25.05 mm, σ = 0.02 mm轴: μ = 25 mm, σ = 0.015 mm
P =

3溶液

1.极限

答案: L = 2.5

极限:先尝试直接代入。若得到 0/0,先化简:因式分解并约去公因式,或乘以共轭式。当 x → ∞ 时,分子分母同除以 x 的最高次幂。洛必达法则:若 f/g 为 0/0 或 ∞/∞ 型,则 lim f/g = lim f′/g′,分子和分母分别求导(不是用商法则)。

除以最高次幂

(5 + 3/x + (−4)/x2)/(2 + 5/x + 3/x2)

令 x → ∞

L = 5/2 = 2.5

答案

L = 2.5

次数相同时,极限为首项系数之比。

2.极限

答案: L = 0.5

极限:先尝试直接代入。若得到 0/0,先化简:因式分解并约去公因式,或乘以共轭式。当 x → ∞ 时,分子分母同除以 x 的最高次幂。洛必达法则:若 f/g 为 0/0 或 ∞/∞ 型,则 lim f/g = lim f′/g′,分子和分母分别求导(不是用商法则)。

第一次应用

(1·ex − 1)/(2x) → 0/0

第二次应用

1·ex/2

代入 x = 0

L = 1·1/2 = 0.5

答案

L = 0.5
3.导数

答案: f′(2) = 729

求导:导数 f′(x) 是图像在 x 处的斜率。幂法则:(xⁿ)′ = n·xⁿ⁻¹。乘积法则:(uv)′ = u′v + uv′。链式法则:(f(g(x)))′ = f′(g(x))·g′(x)。在 f′(x) = 0 处,图像有水平切线:极大值、极小值或鞍点。

链式法则

f′(x) = 3·(3x + 3)2·3

代入

f′(2) = 3·92·3 = 729

答案

f′(2) = 729
4.导数

答案: x = 3

求导:导数 f′(x) 是图像在 x 处的斜率。幂法则:(xⁿ)′ = n·xⁿ⁻¹。乘积法则:(uv)′ = u′v + uv′。链式法则:(f(g(x)))′ = f′(g(x))·g′(x)。在 f′(x) = 0 处,图像有水平切线:极大值、极小值或鞍点。

导数

f′(x) = 3x2 − 6x − 9 = 3·(x2 − 2x − 3)

解 f′(x) = 0

x = −1 或 x = 3

判断哪个是哪个

f″(x) = 6x − 6:在极大值处为负,在极小值处为正。

f″(−1) = −12, f″(3) = 12
x = 3

答案

x = 3
5.积分

答案: I = 14.12

分部积分:∫u·dv = u·v − ∫v·du,由乘积法则得出。把 u 选为求导后会变简单的因式(x 的幂或 ln x),把 dv 选为你能积分的因式(eˣ、sin x、cos x,或 u = ln x 时的 x)。对定积分,要在两个上下限处计算 u·v:∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du。含 x² 的因式需要做两轮。

选取 u 和 dv

u = ln x, du = dx/x; dv = x dx, v = x²/2

原函数

∫x ln x dx = (x²/2)ln x − ∫(x/2) dx = (x²/2)ln x − x²/4

极限

(12.5·ln 5 − 6.25) − (0 − 0.25) = 14.12

答案

I = 14.12
6.积分

答案: V = 53.62

旋转体体积:把 y = f(x) 之下、从 x = a 到 x = b 的区域绕 x 轴旋转,得到一个旋转体,它在 x 处的横截面是半径为 f(x) 的圆盘。把圆盘相加得 V = π∫ₐᵇ f(x)² dx。当区域位于外曲线 R(x) 与内曲线 r(x) 之间时,每个薄片是一个垫圈:V = π∫ₐᵇ (R(x)² − r(x)²) dx。相减之前要先把半径平方。

交点

4x = 2x² → x = 0, x = 2

垫圈法

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

积分

V = π(16·2³/3 − 4·2⁵/5) = π·17.07 = 53.62

答案

V = 53.62
7.复数

答案: a = 0.7765

极坐标形式:z = r∠θ = r(cos θ + j sin θ) = r·ejθ,其中 r = |z|,θ = arg z。换回直角坐标形式:a = r cos θ,b = r sin θ。求辐角需要考虑象限:仅当 a > 0 时 θ = atan(b/a);当 a < 0 时要加上或减去 180°。在极坐标形式下乘法很简单:模相乘,辐角相加;除法则模相除,辐角相减。

分量

a = 3·cos(−75°) = 0.7765
b = 3·sin(−75°) = −2.898

直角坐标形式

z = 0.7765 − j2.898
a = 0.7765

答案

a = 0.7765
8.复数

答案: r = −2

棣莫弗定理与方根:(r∠θ)ⁿ = rⁿ∠nθ。方程 zⁿ = R∠φ 恰有 n 个根,模都是 R1/n,辐角为 (φ + 360°·k)/n,k = 0, 1, …, n − 1:它们均匀分布在一个圆上,相隔 360°/n。当 p² < 4q 时,二次方程 z² + pz + q = 0 有复根 z = −p/2 ± j√(q − p²/4);实系数多项式的复根总是成共轭对 a ± jb 出现。

共轭对

z₂ = 4 − j2
(z − z₁)(z − z₂) = z² − 2·4·z + (4² + 2²) = z² − 8z + 20

提取二次因式

z³ − 6z² + 4z + 40 = (z² − 8z + 20)(z + 2)

实根

r = −2

答案

r = −2

用根之和检验:z₁ + z₂ + r = 8 + (−2) = 6 = −(−6)。

9.微分方程

答案: t = 40.47 min

一阶微分方程:dy/dt = k·y 的解为 y = y₀·ekt(k > 0 时增长,k < 0 时衰减);半衰期 T 对应 k = ln 2/T。牛顿冷却定律 dT/dt = −k·(T − Tₐ) 给出 T(t) = Tₐ + (T₀ − Tₐ)·e−kt。线性方程 y′ + p(x)·y = q(x) 用积分因子 μ = e∫p dx 求解:(μy)′ = μq,所以 μy = ∫μq dx + C。

冷却常数

k = ln(-19/-5)/24 = 0.05563 min⁻¹

达到目标温度的时间

t = ln(-19/(−2))/0.05563 = 40.47 min

答案

t = 40.47 min
10.微分方程

答案: y(0.5) = 0.4171

二阶线性方程:对 y″ + b·y′ + c·y = 0,试 y = ert:得到特征方程 r² + b·r + c = 0。两个不等实根 r₁ ≠ r₂ 给出 y = A·er₁t + B·er₂t(过阻尼);二重根 r 给出 y = (A + Bt)·ert(临界阻尼);复根 −α ± jω 给出 y = e−αt(A cos ωt + B sin ωt)(欠阻尼,振荡)。初值 y(0) 和 y′(0) 决定 A 和 B。写成 y″ + 2ζωny′ + ωn²y = 0 时,ωn 是固有频率,ζ 是阻尼比:ζ < 1 时振荡。在常值输入下,y″ + by′ + cy = c·K 的解是稳态值 K 加上自由响应。

稳态值与根

yss = 1, r = −1 ± 2j

由初值求常数

y(0) = 1 + A = 0 → A = −1
y′(0) = −1A + 2B = 0 → B = −0.5

求值

y(0.5) = 1·(1 − e−0.5(cos 1 + 0.5·sin 1)) = 0.4171

答案

y(0.5) = 0.4171
11.矩阵与行列式

答案: det A = 86

矩阵:2×2 矩阵的行列式为 ad − bc。对于 3×3,按某一行展开:每个元素乘以划去其所在行和列后剩下的 2×2 矩阵的行列式,符号按 + − + 交替。特征值 λ 满足 det(A − λI) = 0;对于 2×2,即 λ² −(迹)·λ + det = 0。

第一行的 2×2 余子式

M₁₁ = 31, M₁₂ = −13, M₁₃ = −30

代数余子式展开

det A = 0·31 − 2·(−13) + (−2)·(−30) = 86

答案

det A = 86
12.矩阵与行列式

答案: λ = 6

矩阵:2×2 矩阵的行列式为 ad − bc。对于 3×3,按某一行展开:每个元素乘以划去其所在行和列后剩下的 2×2 矩阵的行列式,符号按 + − + 交替。特征值 λ 满足 det(A − λI) = 0;对于 2×2,即 λ² −(迹)·λ + det = 0。

特征方程

λ² − (7)λ + (6) = 0

求解

λ = (7 ± √(25))/2 → λ = 1, 6
λ = 6

答案

λ = 6
13.向量:点积和叉积

答案: W = −370 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.

位移

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

功

W = F·d = (−10)·(−5) + (−30)·4 + (−30)·10 = −370 J

答案

W = −370 J

负功:力主要与运动方向相反。

14.向量:点积和叉积

答案: A = 11.73 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.

边向量

AB = B − A = (−5, −5, 0)
AC = C − A = (−8, −6, 3)

叉积

(AB × AC)x = AByACz − ABzACy = (−5)·3 − 0·(−6) = −15
(AB × AC)y = ABzACx − ABxACz = 0·(−8) − (−5)·3 = 15
(AB × AC)z = ABxACy − AByACx = (−5)·(−6) − (−5)·(−8) = −10
|AB × AC| = √((−15)² + 15² + (−10)²) = 23.45

三角形的面积

A = ½·23.45 = 11.73 m²

答案

A = 11.73 m²
15.概率

答案: N = 20

概率:在等可能结果下,P(A) = 有利结果数/总数。独立事件相乘:P(A 且 B) = P(A)·P(B)。补事件常常有用:P(至少一个) = 1 − P(一个都没有)。二项分布:在 n 次独立试验中,成功概率为 p,P(X = k) = C(n, k)·pᵏ·(1 − p)ⁿ⁻ᵏ。

组合

C(6, 3) = (6·5·4)/(3·2·1) = 20

答案

N = 20
16.概率

答案: σY = 11.08

离散随机变量:概率分布列出变量能取的每个值 x 及其概率 p(x);概率介于 0 和 1 之间,且总和为 1。期望值(均值)是 E(X) = Σx·p(x),即长期平均值。方差为 Var(X) = Σ(x − μ)²·p(x) = E(X²) − μ²,标准差为 σ = √Var(X)。对线性变换 Y = aX + b:E(Y) = a·E(X) + b,Var(Y) = a²·Var(X),σY = |a|·σX;加上常数不会改变离散程度。

X 的离散程度

μ = 2.15, Var(X) = 1.228, σX = 1.108

线性变换

σY = |10|·1.108 = 11.08

答案

σY = 11.08

均值确实会改变:E(Y) = 10·2.15 + 5 = 26.5。

17.描述统计

答案: s = 12.18

描述统计:平均数是总和除以 n;中位数是排序后数据的中间值(n 为偶数时取中间两个数的平均)。样本标准差 s = √(Σ(x − x̄)²/(n − 1))。z 分数 z = (x − μ)/σ 表示某个值离平均数多少个标准差。

平均数

x̄ = 255/7 = 36.43

离差平方

Σ(x − x̄)² = 0.5714² + (−2.429)² + 2.571² + 13.57² + 15.57² + (−16.43)² + (−13.43)² = 889.7

标准差

s = √(889.7/6) = 12.18

答案

s = 12.18
18.描述统计

答案: P = 0.1587

描述统计:平均数是总和除以 n;中位数是排序后数据的中间值(n 为偶数时取中间两个数的平均)。样本标准差 s = √(Σ(x − x̄)²/(n − 1))。z 分数 z = (x − μ)/σ 表示某个值离平均数多少个标准差。

标准化

z = −1

标准正态分布表

Φ(−1) = 0.1587
P = 0.1587

答案

P = 0.1587
19.正态分布

答案: P = 0.02275

正态变量的和与差:若 X 和 Y 相互独立且服从正态分布,则 X + Y 与 X − Y 也服从正态分布。均值相加或相减:E(X ± Y) = μX ± μY。方差总是相加,差也是如此:Var(X ± Y) = σX² + σY²,所以 σ = √(σX² + σY²)。标准差从不直接相加。n 个独立同分布项的总和,均值为 n·μ,标准差为 σ·√n。然后照常标准化:z = (x − 均值)/σ,再查 Φ(z)。

和的分布

μ = 105 s, σ = √(5² + 4²) = 6.403 s

标准化

z = (117.8 − 105)/6.403 = 2

上尾

P = 1 − Φ(2) = 1 − 0.9772 = 0.02275

答案

P = 0.02275
20.正态分布

答案: P = 0.9772

正态变量的和与差:若 X 和 Y 相互独立且服从正态分布,则 X + Y 与 X − Y 也服从正态分布。均值相加或相减:E(X ± Y) = μX ± μY。方差总是相加,差也是如此:Var(X ± Y) = σX² + σY²,所以 σ = √(σX² + σY²)。标准差从不直接相加。n 个独立同分布项的总和,均值为 n·μ,标准差为 σ·√n。然后照常标准化:z = (x − 均值)/σ,再查 Φ(z)。

差的分布

μ = 25.05 − 25 = 0.05 mm, σ = 0.025 mm

标准化

z = (0 − 0.05)/0.025 = −2

概率

P = 1 − Φ(−2) = 0.9772

答案

P = 0.9772

ExamProblems · examproblems.com/courses/engineering-mathematics