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

公式集、全トピックの解法の要点、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→∞ (6x2 − 6x + 2)/(5x2 + 5x + 5) を求めなさい。

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

極限 L = limx→0 (e−x − 1 + x)/x² を求めなさい。

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

f(x) = (4x − 1)3 とします。f′(2) を求めなさい。

f(x) = (4x − 1)3x = 2
f′(2) =
4.微分むずかしい
xyx = 1

f(x) = −x2 + 5x + 4 の x = 1 における接線は、y 軸と y = b で交わります。b を求めなさい。

f(x) = −x2 + 5x + 4x = 1
b =
5.積分ふつう
xy12y = x ln x

I = ∫₁2 x·ln x dx を計算しなさい。

∫₁2 x ln x dx
I =
6.積分むずかしい
xy−22

放物線 y = 4 − x² と x 軸で囲まれた部分の面積 A を求めなさい。

y = 4 − x²
A =
7.複素数ふつう
ReIm55−5−5θz

z = 10∠140° を直交形式 a + jb で表し、実部 a を答えてください。

r = 10θ = 140°
a =
8.複素数むずかしい
ReIm22−2−244−4−4z₁z₂

z² − 2z + 17 = 0 を解いてください。解は a ± jb です。実部 a を求めてください。

p = −2q = 17
a =
9.微分方程式ふつう
tyy(0) = 2y(2) = ?

y(0) = 2、y′(0) = -2 として y″ + 6y′ + 5y = 0 を解きなさい。y(2) を求めなさい。

y″ + 6y′ + 5y = 0y(0) = 2y′(0) = -2
y(2) =
10.微分方程式むずかしい
tyK = 2y(0.75) = ?

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

y″ + 3y′ + 27.25y = 54.5y(0) = 0y′(0) = 0
y(0.75) =
11.行列と行列式ふつう
A =−52−1−415−42−1

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

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

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

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

a = (1, −7, 1) と b = (−1, 6, 0) の内積 a·b を求めよ。

a = (1, −7, 1)b = (−1, 6, 0)
a·b =
14.ベクトル:内積と外積むずかしい
xyzABC

頂点 A(0, 0, −2)、B(2, 0, −3)、C(2, 4, −3) をもつ三角形の面積 A を求めよ。座標の単位はメートルである。

A(0, 0, −2)B(2, 0, −3)C(2, 4, −3)
A = m²
15.確率ふつう
P(少なくとも 1 回 6 が出る) = ??1?2P(6) = 各回 1/6

公正なさいころを 2 回振ります。少なくとも 1 回 6 が出る確率 P はいくつですか?

n = 2
P =
16.確率むずかしい
n = 8, p = 0.7: P(X ≤ 1) = ?k012345678

部品を 8 個試験し、各部品は独立に確率 0.7 で故障します。故障が多くとも 1 個である確率 P はいくつですか?

n = 8p = 0.7
P =
17.記述統計ふつう
1520253035404550

データの標本標準偏差 s を求めなさい:11, 21, 53, 25, 24。

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

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

μ = 500σ = 100x = 380
z =
19.正規分布ふつう
27.242.8P = ?

配達時間は、平均 32 min、標準偏差 6 min の正規分布に従います。27.2 から 42.8 min の間にある割合はいくらですか?P を小数で答えなさい。

μ = 32σ = 6x1 = 27.2x2 = 42.8
P =
20.正規分布むずかしい
380420460500540580620h上端 60%

電池の寿命 は、平均 500 h、標準偏差 40 h の正規分布に従います。色つきの上位 60% が右側になるように、しきい値をドラッグします。しきい値 x を求めなさい。

μ = 500 hσ = 40 h60%
x = h

3溶液

1.極限

答え: L = 1.2

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

最高次の累乗で割る

(6 + (−6)/x + 2/x2)/(5 + 5/x + 5/x2)

x → ∞ とする

L = 6/5 = 1.2

答え

L = 1.2

次数が等しければ、極限は最高次の係数の比になります。

2.極限

答え: L = 0.5

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

1 回目の適用

(−1·e−x + 1)/(2x) → 0/0

2 回目の適用

1·e−x/2

x = 0 を代入する

L = 1·1/2 = 0.5

答え

L = 0.5
3.微分

答え: f′(2) = 588

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

連鎖律

f′(x) = 3·(4x − 1)2·4

代入する

f′(2) = 3·72·4 = 588

答え

f′(2) = 588
4.微分

答え: b = 5

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

接点と傾き

f(1) = 8
f′(x) = −2x + 5 → f′(1) = 3

x = 0 での接線の値

b = f(x₀) − f′(x₀)·x₀ = 8 − 3·1 = 5

答え

b = 5
5.積分

答え: I = 0.6363

部分積分:∫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² があるときは 2 回行います。

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

極限

(2·ln 2 − 1) − (0 − 0.25) = 0.6363

答え

I = 0.6363
6.積分

答え: A = 10.67

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

零点

4 − x² = 0 → x = ±2

積分する

A = [4x − x³/3]−22 = 5.333 − (−5.333) = 10.67

答え

A = 10.67
7.複素数

答え: a = −7.66

極形式: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° を足すか引きます。極形式では掛け算が簡単です:絶対値を掛けて偏角を足します。割り算では絶対値を割って偏角を引きます。

成分

a = 10·cos(140°) = −7.66
b = 10·sin(140°) = 6.428

直交形式

z = −7.66 + j6.428
a = −7.66

答え

a = −7.66
8.複素数

答え: a = 1

ド・モアブルの定理と累乗根:(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 として現れます。

判別式

p² − 4q = 4 − 68 = −64 < 0

2 次方程式の解の公式

z = (2 ± j√64)/2 = 1 + j4, 1 − j4

答え

a = 1

答え

a = 1
9.微分方程式

答え: y(2) = 0.2707

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 に自由応答を加えたものです。

特性方程式の根

(r + 1)(r + 5) = 0 → r = −1, −5

初期値から求める定数

A + B = 2, −1A − 5B = -2
A = 2, B = 0

計算する

y(2) = 2·e−2 + 0·e−10 = 0.2707

答え

y(2) = 0.2707
10.微分方程式

答え: y(0.75) = 2.644

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 = 2, r = −1.5 ± 5j

初期値から求める定数

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

計算する

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

答え

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

答え: det A = 11

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

第 1 行の 2×2 小行列式

M₁₁ = −11, M₁₂ = 24, M₁₃ = −4

余因子展開

det A = (−5)·(−11) − 2·24 + (−1)·(−4) = 11

答え

det A = 11
12.行列と行列式

答え: λ = −1

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

固有方程式

λ² − (−3)λ + (2) = 0

解く

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

答え

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

答え: a·b = −43

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 = 1·(−1) + (−7)·6 + 1·0 = (−1) + (−42) + 0 = −43

答え

a·b = −43

負:2 つのベクトルのなす角は 90° より大きい。

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

答え: A = 4.472 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 = (2, 0, −1)
AC = C − A = (2, 4, −1)

外積

(AB × AC)x = AByACz − ABzACy = 0·(−1) − (−1)·4 = 4
(AB × AC)y = ABzACx − ABxACz = (−1)·2 − 2·(−1) = 0
(AB × AC)z = ABxACy − AByACx = 2·4 − 0·2 = 8
|AB × AC| = √(4² + 0² + 8²) = 8.944

三角形の面積

A = ½·8.944 = 4.472 m²

答え

A = 4.472 m²
15.確率

答え: P = 0.3056

確率:どの結果も同様に確からしいとき、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)ⁿ⁻ᵏ。

余事象

P(6 が出ない) = (5/6)2 = 0.6944
P = 1 − 0.6944 = 0.3056

答え

P = 0.3056
16.確率

答え: P = 0.00129

確率:どの結果も同様に確からしいとき、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)ⁿ⁻ᵏ。

故障 0 個と 1 個

P(X = 0) = 0.38 = 6.561 × 10⁻⁵
P(X = 1) = 8·0.7·0.37 = 0.001225

足し合わせる

P = 6.561 × 10⁻⁵ + 0.001225 = 0.00129

答え

P = 0.00129
17.記述統計

答え: s = 15.66

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

平均

x̄ = 134/5 = 26.8

偏差の 2 乗

Σ(x − x̄)² = (−15.8)² + (−5.8)² + 26.2² + (−1.8)² + (−2.8)² = 980.8

標準偏差

s = √(980.8/4) = 15.66

答え

s = 15.66
18.記述統計

答え: z = −1.2

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

標準化する

z = (380 − 500)/100 = −1.2

答え

z = −1.2
19.正規分布

答え: P = 0.7522

正規分布:値 x を z スコア z = (x − μ)/σ に変換します。これは平均から標準偏差何個分離れているかを表します。z 表(または Φ(z))は z より左側の面積、つまり P(X < x) を与えます。右側の面積は 1 − Φ(z)、2 つの値の間の面積はそれぞれの左側の面積の差です。

z スコア

z₁ = (27.2 − 32)/6 = −0.8
z₂ = (42.8 − 32)/6 = 1.8

間の面積

P = Φ(1.8) − Φ(−0.8) = 0.9641 − 0.2119 = 0.7522

答え

P = 0.7522
20.正規分布

答え: x = 490 h

面積から 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.4

表からの z

z ≈ −0.25

元の単位に戻す

x = 500 − 0.25·40 = 490 h

答え

x = 490 h

ExamProblems · examproblems.com/courses/engineering-mathematics