テイラー多項式:x = a の近くで、滑らかな関数は Pₙ(x) = f(a) + f′(a)(x − a) + f″(a)(x − a)²/2! + … + f⁽ⁿ⁾(a)(x − a)ⁿ/n! で近似されます。a = 0 のときがマクローリン多項式です。標準的な級数: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ᵏ。ラグランジュの誤差評価:|f(x) − Pₙ(x)| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!。ここで M は a と x の間での |f⁽ⁿ⁺¹⁾| の最大値です。
n を増やしたときの評価値
n = 12: 1.052 × 10⁻⁵ ≥ 1 × 10⁻⁵
n = 13: 1.503 × 10⁻⁶ < 1 × 10⁻⁵
最小の次数
n = 13
答え
n = 13
5.多変数の微分積分
答え: Duf = −2.774
勾配:∇f = (fx, fy) は最も急に増加する方向を向き、|∇f| は変化率の最大値です。方向 v に沿った方向微分は、単位ベクトル u = v/|v| を用いて Duf = ∇f·u です。点 (a, b) での接平面は z = f(a, b) + fx(a, b)(x − a) + fy(a, b)(y − b) で、その点の近くでの f の線形近似を与えます。連鎖律:dz/dt = fx·dx/dt + fy·dy/dt。
その点での勾配
fx = −2y + 3x² → fx(2, −1) = 14
fy = −2x − 10y → fy(2, −1) = 6
単位ベクトル
|v| = √((−2)² + 3²) = 3.606
u = (−0.5547, 0.8321)
内積
Duf = 14·(−0.5547) + 6·0.8321 = −2.774
答え
Duf = −2.774
6.多変数の微分積分
答え: D = −24
臨界点と二重積分:臨界点では fx = fy = 0 です。D = fxxfyy − fxy² を用いた 2 階微分判定法:D > 0 かつ fxx > 0 なら極小、D > 0 かつ fxx < 0 なら極大、D < 0 なら鞍点、D = 0 なら何も決まりません。一般の領域上の二重積分では、内側の積分の限界を領域を囲む曲線にします。円板や扇形では極座標 x = r cos θ、y = r sin θ、dA = r dr dθ を使います(係数 r を忘れないでください)。
2 階偏導関数
fxx = 12, fyy = −2, fxy = 0
判別式
D = 12·(−2) − 0² = −24
→ 鞍点
答え
D = −24
7.微分方程式
答え: y(1) = 2.943
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)² = 0 → r = −1
定数
A = y(0) = 2
y′(0) = B − 1A → B = 4 + 2 = 6
計算する
y(1) = (2 + 6·1)·e−1 = 2.943
答え
y(1) = 2.943
8.微分方程式
答え: ymax = 8.376
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 に自由応答を加えたものです。