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

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

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

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

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

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

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

f(x) = x3 − 15x2 + 63x − 6 とします。極大点の x 座標 x を求めなさい。

f(x) = x3 − 15x2 + 63x − 6
x =
5.積分ふつう
xy1y = √(x)

y = √(x) の下の x = 0 から x = 1 までの領域を x 軸のまわりに回転させます。立体(放物面体)の体積 V を求めなさい。

y = √(x)0 ≤ x ≤ 1
V =
6.積分むずかしい
xyy = x2​y = 2x

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

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

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

r = 10θ = 140°
a =
8.複素数むずかしい
ReIm11−1−122−2−2z₁z₂

z³ − 7z² + 20z − 24 = 0 の解の 1 つは z₁ = 2 + j2 です。係数は実数です。実数の解 r を求めてください。

z₁ = 2 + j2
r =
9.微分方程式ふつう
t (min)T (°C)51Ta​T = ?

340 °C の鍛造鋼部品を 23 °C の作業場に置いて冷まします。 温度はニュートンの冷却法則 dT/dt = −k·(T − Tₐ) に従います。 k = 0.02 min⁻¹ として、51 分後の温度 T を求めなさい。

T0 = 340 °CTa = 23 °Ck = 0.02 min⁻¹t = 51 min
T = °C
10.微分方程式むずかしい
t (min)m (kg)t = ?41 kg

タンクに、塩 20 kg を溶かした水が 500 L 入っています。 1 L あたり 0.1 kg の塩を含む塩水が 20 L/min で流入し、よく混ざった溶液が同じ流量で流出します。 タンク内の塩が 41 kg になるのは何分後 t ですか?

V = 500 Lm0 = 20 kgc = 0.1 kg/Lr = 20 L/min
t = min
11.行列と行列式ふつう
A =43−45B =621−6

示された行列について C = AB とします。成分 c11(第 1 行、第 1 列)を求めなさい。

A = [4 3; −4 5]B = [6 2; 1 −6]
c11 =
12.行列と行列式むずかしい
A =33−8−7

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

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

ベクトル a = (−2, −5, −6) と b = (−2, 2, −6) のなす角 θ を度で求めよ。

a = (−2, −5, −6)b = (−2, 2, −6)
θ = °
14.ベクトル:内積と外積むずかしい
xyzabc

1 つの頂点から出る辺が a = (−4, 3, 1)、b = (−1, 1, 4)、c = (2, −3, 0) である平行六面体の体積 V を求めよ。

a = (−4, 3, 1)b = (−1, 1, 4)c = (2, −3, 0)
V =
15.確率ふつう
n = 6k = 2N = ?

6 人のグループから 2 人の委員会を選ぶ方法 N は何通りですか?順序は問いません。

n = 6k = 2
N =
16.確率むずかしい
n = 9, p = 0.1: P(X ≤ 1) = ?k0123456789

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

n = 9p = 0.1
P =
17.記述統計ふつう
30354045505560

データの標本標準偏差 s を求めなさい:44, 39, 26, 60, 51。

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

ある量は平均 70、標準偏差 5 の正規分布に従います。75 より大きい値の割合 P はいくつですか?

μ = 70σ = 5x = 75
P =
19.正規分布ふつう
49.567P = ?

試験の得点は、平均 65 点、標準偏差 10 点 の正規分布に従います。49.5 から 67 点 の間にある割合はいくらですか?P を小数で答えなさい。

μ = 65σ = 10x1 = 49.5x2 = 67
P =
20.正規分布むずかしい
0.03 mm0D − d?

穴の直径 D は平均 25.03 mm、標準偏差 0.02 mm の正規分布に従います。軸の直径 d は D と独立で、平均 25 mm、標準偏差 0.015 mm の正規分布に従います。無作為に選んだ軸は D − d > 0 のとき穴にはまります。はまる確率 P を求めなさい。z は小数第 2 位に丸めなさい。

穴: μ = 25.03 mm, σ = 0.02 mm軸: μ = 25 mm, σ = 0.015 mm
P =

3溶液

1.極限

答え: L = −1.25

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

共役な式を掛ける

f(x) = −5x/(√(4x2 − 5x) + 2x)

x で割る

f(x) = −5/(√(4 + (−5)/x) + 2)

x → ∞ とする

L = −5/(2 + 2) = −1.25

答え

L = −1.25
2.極限

答え: L = 3.125

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

1 回目の適用

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

2 回目の適用

25·cos 5x/8

x = 0 を代入する

L = 25/8 = 3.125

答え

L = 3.125
3.微分

答え: f′(2) = 0.5413

微分:導関数 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−x·(3x2 + (−1)·x3)

代入する

f′(2) = 0.1353·(12 + (−8)) = 0.5413

答え

f′(2) = 0.5413
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 − 30x + 63 = 3·(x2 − 10x + 21)

f′(x) = 0 を解く

x = 3 または x = 7

どちらがどちらか

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

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

答え

x = 3
5.積分

答え: V = 1.571

回転体の体積:x = a から x = b までの y = f(x) の下の領域を x 軸のまわりに回転させると、x での断面が半径 f(x) の円板である立体ができます。円板を足し合わせると V = π∫ₐᵇ f(x)² dx になります。領域が外側の曲線 R(x) と内側の曲線 r(x) の間にあるときは、各スライスはワッシャー(円環)で、V = π∫ₐᵇ (R(x)² − r(x)²) dx です。引く前に半径を 2 乗しましょう。

円板法

V = π∫₀1 x dx

積分する

V = π·1·1²/2 = 1.571

答え

V = 1.571
6.積分

答え: A = 1.333

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

交点

x2 = 2x → x = 0, x = 2

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

A = [x2 − 0.3333x3]02 = 1.333

答え

A = 1.333
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.複素数

答え: r = 3

ド・モアブルの定理と累乗根:(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₂ = 2 − j2
(z − z₁)(z − z₂) = z² − 2·2·z + (2² + 2²) = z² − 4z + 8

2 次式でくくり出す

z³ − 7z² + 20z − 24 = (z² − 4z + 8)(z − 3)

実数の解

r = 3

答え

r = 3

解の和で確かめましょう:z₁ + z₂ + r = 4 + 3 = 7 = −(−7)。

9.微分方程式

答え: T = 137.3 °C

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 です。

冷却法則の解

T(t) = 23 + 317·e−0.02t

t を代入する

T(51) = 23 + 317·0.3606 = 137.3 °C

答え

T = 137.3 °C
10.微分方程式

答え: t = 30.1 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 です。

結果

m(t) = 50 + (−30)·e−t/25

t について解く

e−t/25 = (50 − 41)/(50 − 20) = 0.3
t = 25·ln(3.333) = 30.1 min

答え

t = 30.1 min
11.行列と行列式

答え: c11 = 27

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

行 × 列

c11 = 4·6 + 3·1 = 27

答え

c11 = 27
12.行列と行列式

答え: λ = −1

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

固有方程式

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

解く

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

答え

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

答え: θ = 55.88 °

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 = (−2)·(−2) + (−5)·2 + (−6)·(−6) = 30

大きさ

|a| = √((−2)² + (−5)² + (−6)²) = 8.062
|b| = √((−2)² + 2² + (−6)²) = 6.633

位相角

cos θ = 30/(8.062·6.633) = 0.561
θ = arccos(0.561) = 55.88°

答え

θ = 55.88°
14.ベクトル:内積と外積

答え: 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.

外積 b × c

(b × c)x = bycz − bzcy = 1·0 − 4·(−3) = 12
(b × c)y = bzcx − bxcz = 4·2 − (−1)·0 = 8
(b × c)z = bxcy − bycx = (−1)·(−3) − 1·2 = 1

a との内積

a·(b × c) = (−4)·12 + 3·8 + 1·1 = −23
V = |−23| = 23

答え

V = 23
15.確率

答え: N = 15

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

組合せ

C(6, 2) = (6·5)/(2·1) = 15

答え

N = 15
16.確率

答え: P = 0.7748

確率:どの結果も同様に確からしいとき、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.99 = 0.3874
P(X = 1) = 9·0.1·0.98 = 0.3874

足し合わせる

P = 0.3874 + 0.3874 = 0.7748

答え

P = 0.7748
17.記述統計

答え: s = 12.79

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

平均

x̄ = 220/5 = 44

偏差の 2 乗

Σ(x − x̄)² = 0² + (−5)² + (−18)² + 16² + 7² = 654

標準偏差

s = √(654/4) = 12.79

答え

s = 12.79
18.記述統計

答え: P = 0.1587

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

標準化する

z = 1

標準正規分布表

Φ(1) = 0.8413
P = 1 − 0.8413 = 0.1587

答え

P = 0.1587
19.正規分布

答え: P = 0.5187

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

z スコア

z₁ = (49.5 − 65)/10 = −1.55
z₂ = (67 − 65)/10 = 0.2

間の面積

P = Φ(0.2) − Φ(−1.55) = 0.5793 − 0.06057 = 0.5187

答え

P = 0.5187
20.正規分布

答え: P = 0.8849

正規変数の和と差: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.03 − 25 = 0.03 mm, σ = 0.025 mm

標準化する

z = (0 − 0.03)/0.025 = −1.2

確率

P = 1 − Φ(−1.2) = 0.8849

答え

P = 0.8849

ExamProblems · examproblems.com/courses/engineering-mathematics