工業数学:直前対策パック

公式集、全トピックの解法の要点、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 の 2 次方程式
z = −p/2 ± j√(q − p²/4)

解き方

  1. 演算に合った形を選びましょう:足し算と引き算には直交形式 a + jb、掛け算・割り算・冪・累乗根には極形式 r∠θ です。
  2. 展開するときは j² = −1 を使います。割り算は、分子と分母に分母の共役複素数を掛けて行います。
  3. 変換には r = √(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. 垂直なベクトル、面積、モーメントには外積を使います:3 つの成分を慎重に書き、順序 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 − μ)/σ
正規分布
±1σ 以内に 68%、±2σ 以内に 95%

正規分布

z スコア
z = (x − μ)/σ
x 未満
P(X < x) = Φ(z)
x より大きい
P(X > x) = 1 − Φ(z)
2 つの値の間
Φ(z2) − Φ(z1)
パーセンタイル
x = μ + z·σ

解き方

  1. 値を z スコアに変換します:z = (x − μ)/σ。
  2. z より左側の面積 Φ(z) を表で調べます。
  3. 右側の面積は 1 − Φ(z)、区間の面積は 2 つの左側の面積の差です。
  4. パーセンタイルでは逆にたどります:面積から z を求め、x = μ + z·σ とします。

2演習問題

1.極限ふつう
xyy = L

極限 L = limx→∞ (√(x2 + 2x) − x) を求めなさい。

f(x) = √(x2 + 2x) − xx → ∞
L =
2.極限むずかしい
xyx = 0

極限 L = limx→0 (1 − cos 3x)/(4x²) を求めなさい。x の単位はラジアンとします。

f(x) = (1 − cos 3x)/(4x²)x → 0
L =
3.微分ふつう
xyx = 1

f(x) = x3·e−0.5x とします。f′(1) を求めなさい。

f(x) = x3·e−0.5xx = 1
f′(1) =
4.微分むずかしい
xy

f(x) = x3 + 3x2 − 9x − 2 とします。極小点の x 座標 x を求めなさい。

f(x) = x3 + 3x2 − 9x − 2
x =
5.積分ふつう
xy12

I = ∫12 1/x dx を計算しなさい。

f(x) = 1/x[1, 2]
I =
6.積分むずかしい
xyy = 2x2​y = 3x

y = 2x2 と y = 3x で囲まれた部分の面積 A を求めなさい。

y = 2x2y = 3x
A =
7.複素数ふつう
ReIm22−2−244−4−4z₁z₂w

z₁ = 2.5∠150°、z₂ = 0.5∠45° のとき、w = z₁/z₂ の虚部を求めてください。

z₁ = 2.5∠150°z₂ = 0.5∠45°
Im(w) =
8.複素数むずかしい
ReIm22−2−244−4−4ABCD

点 A–D のうち、ちょうど 1 つが z⁶ = 729∠120° の解です。それをクリックしてください。

w = 729∠120°n = 6
  1. 1A (3∠20°)
  2. 2B (3∠50°)
  3. 3C (5.4∠50°)
  4. 4D (3∠120°)
9.微分方程式ふつう
t (min)T (°C)47Ta​70 °C49 °C

300 °C の鍛造鋼部品を 28 °C の作業場に置いて冷まします。 47 分後の温度は 70 °C です。 温度はニュートンの冷却法則 dT/dt = −k·(T − Tₐ) に従います。 開始から 49 °C に達するまでの時間 t を求めなさい。

T0 = 300 °CTa = 28 °CT(47) = 70 °CT = 49 °C
t = min
10.微分方程式むずかしい
tyK = 10y(0.75) = ?

系が t ≥ 0 で y″ + y′ + 25.25y = 252.5 に従い、静止状態 y(0) = 0、y′(0) = 0 から始まります。y(0.75) を求めなさい。角度はラジアンとします。

y″ + y′ + 25.25y = 252.5y(0) = 0y′(0) = 0
y(0.75) =
11.行列と行列式ふつう
A =50−224410−3

示された 3×3 行列 A の行列式 det A を求めなさい。

A = [5 0 −2; 2 4 4; 1 0 −3]
det A =
12.行列と行列式むずかしい
A =−1421

示された行列 A の大きい方の固有値 λ を求めなさい。

A = [−1 4; 2 1]
λ =
13.ベクトル:内積と外積ふつう
xyzab

a = (3, 0, −2) と b = (0, −5, −6) の内積 a·b を求めよ。

a = (3, 0, −2)b = (0, −5, −6)
a·b =
14.ベクトル:内積と外積むずかしい
xyzab

a = (5, 5, 5)、b = (−5, −6, 3) のとき、外積 a × b の x 成分 (a × b)x を求めよ。

a = (5, 5, 5)b = (−5, −6, 3)
(a × b)x =
15.確率ふつう
赤: 6青: 6

袋に赤玉 6 個と青玉 6 個が入っています。1 個目を戻さずに 2 個取り出します。2 個とも赤である確率 P はいくつですか?

6 / 6
P =
16.確率むずかしい
€p(x)−10.8510.1190.05

ゲームの参加費は €1 です。確率 0.1 で €2、確率 0.05 で €20 を獲得し、それ以外では何も獲得しません。賭け金は戻りません。1 ゲームあたりの純利得の期待値 E(G) を求めなさい。

賭け金 = €1€2: p = 0.1€20: p = 0.05
E(G) = €
17.記述統計ふつう
15202530354045505560

データの標本標準偏差 s を求めなさい:23, 59, 14, 60, 20, 40, 36。

n = 7
s =
18.記述統計むずかしい
z −2−1012

試験の点数は平均 100、標準偏差 15 の正規分布に従います。70 点の z スコア z を求めなさい。

μ = 100σ = 15x = 70
z =
19.正規分布ふつう
−3−2−10123zz の右側の色つき面積 = 0.7054

標準正規曲線の下で、z の右側の色つき面積が 0.7054 になるように z をドラッグします。z を小数第 2 位まで求めなさい。

A = 0.7054
z =
20.正規分布むずかしい
157164171178185192199cm下端 20%

成人男性の身長 は、平均 178 cm、標準偏差 7 cm の正規分布に従います。色つきの下位 20% が左側になるように、しきい値をドラッグします。しきい値 x を求めなさい。

μ = 178 cmσ = 7 cm20%
x = cm

3溶液

1.極限

答え: L = 1

極限:まず直接代入してみます。0/0 になる場合は、先に式を簡単にします:因数分解して共通因数を約分するか、共役な式を掛けます。x → ∞ では、分子と分母を x の最高次の累乗で割ります。ロピタルの定理:f/g が 0/0 または ∞/∞ になるなら、lim f/g = lim f′/g′ です。分子と分母を別々に微分します(商の微分法は使いません)。

共役な式を掛ける

f(x) = 2x/(√(x2 + 2x) + x)

x で割る

f(x) = 2/(√(1 + 2/x) + 1)

x → ∞ とする

L = 2/(1 + 1) = 1

答え

L = 1
2.極限

答え: L = 1.125

極限:まず直接代入してみます。0/0 になる場合は、先に式を簡単にします:因数分解して共通因数を約分するか、共役な式を掛けます。x → ∞ では、分子と分母を x の最高次の累乗で割ります。ロピタルの定理:f/g が 0/0 または ∞/∞ になるなら、lim f/g = lim f′/g′ です。分子と分母を別々に微分します(商の微分法は使いません)。

1 回目の適用

3·sin 3x/(8x) → 0/0

2 回目の適用

9·cos 3x/8

x = 0 を代入する

L = 9/8 = 1.125

答え

L = 1.125
3.微分

答え: f′(1) = 1.516

微分:導関数 f′(x) は x におけるグラフの傾きです。べき乗の法則:(xⁿ)′ = n·xⁿ⁻¹。積の法則:(uv)′ = u′v + uv′。連鎖律:(f(g(x)))′ = f′(g(x))·g′(x)。f′(x) = 0 のところではグラフの接線が水平になり、極大・極小・鞍点のいずれかです。

積の法則

f′(x) = u′v + uv′ = e−0.5x·(3x2 + (−0.5)·x3)

代入する

f′(1) = 0.6065·(3 + (−0.5)) = 1.516

答え

f′(1) = 1.516
4.微分

答え: x = 1

微分:導関数 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 = −3 または x = 1

どちらがどちらか

f″(x) = 6x + 6:極大では負、極小では正。

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

答え

x = 1
5.積分

答え: I = 0.6931

定積分:∫ₐᵇ f(x) dx は、a から b までのグラフと x 軸の間の符号付き面積です。原始関数 F(F′ = f)を求め、F(b) − F(a) を計算します。べき乗の法則:∫xⁿ dx = xⁿ⁺¹/(n + 1)。2 曲線の間の面積は ∫(上 − 下) dx です。

原始関数

F(x) = ln x

上端と下端を代入する

I = (ln 2 − ln 1) = ln 2 = 0.6931

答え

I = 0.6931
6.積分

答え: A = 1.125

定積分:∫ₐᵇ f(x) dx は、a から b までのグラフと x 軸の間の符号付き面積です。原始関数 F(F′ = f)を求め、F(b) − F(a) を計算します。べき乗の法則:∫xⁿ dx = xⁿ⁺¹/(n + 1)。2 曲線の間の面積は ∫(上 − 下) dx です。

交点

2x2 = 3x → x = 0, x = 1.5

上から下を引いて積分する

A = [1.5x2 − 0.6667x3]01.5 = 1.125

答え

A = 1.125
7.複素数

答え: Im(w) = 4.83

極形式:z = r∠θ = r(cos θ + j sin θ) = r·ejθ、r = |z|、θ = arg z です。直交形式に戻すには:a = r cos θ、b = r sin θ。偏角には象限の考慮が必要です:θ = atan(b/a) が使えるのは a > 0 のときだけで、a < 0 のときは 180° を足すか引きます。極形式では掛け算が簡単です:絶対値を掛けて偏角を足します。割り算では絶対値を割って偏角を引きます。

w の絶対値と偏角

|w| = 2.5 / 0.5 = 5
θ = 150° − 45° = 105°

直交形式

Re(w) = 5·cos(105°) = −1.294
Im(w) = 5·sin(105°) = 4.83
Im(w) = 4.83

答え

Im(w) = 4.83
8.複素数

答え: A (3∠20°)

ド・モアブルの定理と累乗根:(r∠θ)ⁿ = rⁿ∠nθ。方程式 zⁿ = R∠φ にはちょうど n 個の解があり、すべて絶対値は R1/n で、偏角は (φ + 360°·k)/n(k = 0, 1, …, n − 1)です:円周上に 360°/n ずつ等間隔に並びます。p² < 4q のとき、2 次方程式 z² + pz + q = 0 は複素数の解 z = −p/2 ± j√(q − p²/4) をもちます。実数係数の多項式の複素数の解は、必ず共役な対 a ± jb として現れます。

解の絶対値と偏角

|z| = 3
θ = 20°, 80°, 140°, −160°, −100°, −40°

その点

A = 3∠20°

答え

A (3∠20°)
9.微分方程式

答え: t = 64.44 min

1 階微分方程式: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(272/42)/47 = 0.03975 min⁻¹

目標温度に達するまでの時間

t = ln(272/21)/0.03975 = 64.44 min

答え

t = 64.44 min
10.微分方程式

答え: y(0.75) = 16.03

2 階線形方程式:y″ + b·y′ + c·y = 0 では y = ert とおきます。これから特性方程式 r² + b·r + c = 0 が得られます。異なる 2 つの実根 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 = 10, r = −0.5 ± 5j

初期値から求める定数

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

計算する

y(0.75) = 10·(1 − e−0.375(cos 3.75 + 0.1·sin 3.75)) = 16.03

答え

y(0.75) = 16.03
11.行列と行列式

答え: det A = −52

行列:2×2 行列の行列式は ad − bc です。3×3 では 1 つの行に沿って展開します:各成分に、その行と列を除いた 2×2 行列の行列式を掛け、符号は + − + と交互にします。固有値 λ は det(A − λI) = 0 を満たし、2×2 では λ² −(トレース)·λ + det = 0 です。

第 1 行の 2×2 小行列式

M₁₁ = −12, M₁₂ = −10, M₁₃ = −4

余因子展開

det A = 5·(−12) − 0·(−10) + (−2)·(−4) = −52

答え

det A = −52
12.行列と行列式

答え: λ = 3

行列:2×2 行列の行列式は ad − bc です。3×3 では 1 つの行に沿って展開します:各成分に、その行と列を除いた 2×2 行列の行列式を掛け、符号は + − + と交互にします。固有値 λ は det(A − λI) = 0 を満たし、2×2 では λ² −(トレース)·λ + det = 0 です。

固有方程式

λ² − (0)λ + (−9) = 0

解く

λ = (0 ± √(36))/2 → λ = −3, 3
λ = 3

答え

λ = 3
13.ベクトル:内積と外積

答え: a·b = 12

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.

対応する成分を掛けて足す

a·b = 3·0 + 0·(−5) + (−2)·(−6) = 0 + 0 + 12 = 12

答え

a·b = 12

正:2 つのベクトルのなす角は 90° より小さい。

14.ベクトル:内積と外積

答え: (a × b)x = 45

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.

外積を成分ごとに計算する

(a × b)x = aybz − azby = 5·3 − 5·(−6) = 45
(a × b)y = azbx − axbz = 5·(−5) − 5·3 = −40
(a × b)z = axby − aybx = 5·(−6) − 5·(−5) = −5

求める成分

(a × b)x = 45

答え

(a × b)x = 45

検算:a × b は両方のベクトルに垂直なので、a および b との内積は 0 です。

15.確率

答え: P = 0.2273

確率:どの結果も同様に確からしいとき、P(A) = 該当する場合の数/全体の場合の数。独立な事象は掛け合わせます:P(A かつ B) = P(A)·P(B)。余事象がよく役立ちます:P(少なくとも 1 回) = 1 − P(1 回もない)。二項分布:成功確率 p の独立な試行を n 回行うと、P(X = k) = C(n, k)·pᵏ·(1 − p)ⁿ⁻ᵏ。

1 回目と 2 回目

P = 6/12 · 5/11 = 30/132 = 0.2273

答え

P = 0.2273
16.確率

答え: E(G) = 0.2 €

離散確率変数:確率分布は、変数がとりうるすべての値 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 となり、定数を加えてもばらつきは変わりません。

純利得の結果

G: 1 (p = 0.1), 19 (p = 0.05), −1 (p = 0.85)

期待値

E(G) = 1·0.1 + 19·0.05 − 1·0.85 = 0.2 €

答え

E(G) = 0.2 €

正の値は、平均して得をすることを意味します。

17.記述統計

答え: s = 18.39

記述統計:平均は合計を n で割ったもの、中央値は並べ替えたデータの真ん中の値です(n が偶数なら真ん中 2 つの平均)。標本標準偏差は s = √(Σ(x − x̄)²/(n − 1))。z スコア z = (x − μ)/σ は、値が平均から標準偏差何個分離れているかを表します。

平均

x̄ = 252/7 = 36

偏差の 2 乗

Σ(x − x̄)² = (−13)² + 23² + (−22)² + 24² + (−16)² + 4² + 0² = 2030

標準偏差

s = √(2030/6) = 18.39

答え

s = 18.39
18.記述統計

答え: z = −2

記述統計:平均は合計を n で割ったもの、中央値は並べ替えたデータの真ん中の値です(n が偶数なら真ん中 2 つの平均)。標本標準偏差は s = √(Σ(x − x̄)²/(n − 1))。z スコア z = (x − μ)/σ は、値が平均から標準偏差何個分離れているかを表します。

標準化する

z = (70 − 100)/15 = −2

答え

z = −2
19.正規分布

答え: z = −0.54

面積から z へ戻る:Φ(z) は z の左側の面積です。左側の裾の面積 A なら、表で Φ(z) = A となる z を見つけます。右側の裾の面積 A なら Φ(z) = 1 − A です。−z から z の間の中央の面積 A なら、各裾に (1 − A)/2 が残るので、Φ(z) = (1 + A)/2 です。値 x には x = μ + zσ を使います。

z の左側の面積

Φ(z) = 1 − 0.7054 = 0.2946

z 表を逆に読む

表の中で 0.2946 を見つけ、z を読み取ります。

z = −0.54

答え

z = −0.54
20.正規分布

答え: x = 172.1 cm

面積から z へ戻る:Φ(z) は z の左側の面積です。左側の裾の面積 A なら、表で Φ(z) = A となる z を見つけます。右側の裾の面積 A なら Φ(z) = 1 − A です。−z から z の間の中央の面積 A なら、各裾に (1 − A)/2 が残るので、Φ(z) = (1 + A)/2 です。値 x には x = μ + zσ を使います。

左側の面積

Φ(z) = 0.2

表からの z

z ≈ −0.84

元の単位に戻す

x = 178 − 0.84·7 = 172.1 cm

答え

x = 172.1 cm

ExamProblems · examproblems.com/courses/engineering-mathematics