微分方程式
- 指数関数的な増加・減少
- 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
公式集、全トピックの解法の要点、8 問の演習問題。完全な解答は巻末にあります。
460 °C の鍛造鋼部品を 22 °C の作業場に置いて冷まします。 50 分後の温度は 81 °C です。 温度はニュートンの冷却法則 dT/dt = −k·(T − Tₐ) に従います。 定数 k を min⁻¹ で求めなさい。
y′ + 0.5y = 9e−3t、y(0) = 4 を解きなさい。y(1) を求めなさい。
I = ∫02 2x·(x² + 1)1 dx を計算しなさい。
y = 0.5x2 と y = 2x で囲まれた部分の面積 A を求めなさい。
f(x) = ln(x2 + 5) とします。f′(3) を求めなさい。
f(x) = −x2 − 2x − 3 の x = 3 における接線は、y 軸と y = b で交わります。b を求めなさい。
台形公式を n = 4 で使って ∫02 sin x dx を近似し、厳密な値と比べてください。絶対誤差 |E| = |T − I| を求めてください。
ホイン法(改良オイラー法)を h = 0.5 で y′ = e−ty、y(0) = 0.5 に使い、y3 ≈ y(1.5) を推定してください。
答え: k = 0.04009 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 です。
答え: y(1) = 4.43
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 です。
答え: I = 12
定積分:∫ₐᵇ f(x) dx は、a から b までのグラフと x 軸の間の符号付き面積です。原始関数 F(F′ = f)を求め、F(b) − F(a) を計算します。べき乗の法則:∫xⁿ dx = xⁿ⁺¹/(n + 1)。2 曲線の間の面積は ∫(上 − 下) dx です。
答え: A = 5.333
定積分:∫ₐᵇ f(x) dx は、a から b までのグラフと x 軸の間の符号付き面積です。原始関数 F(F′ = f)を求め、F(b) − F(a) を計算します。べき乗の法則:∫xⁿ dx = xⁿ⁺¹/(n + 1)。2 曲線の間の面積は ∫(上 − 下) dx です。
答え: f′(3) = 0.4286
微分:導関数 f′(x) は x におけるグラフの傾きです。べき乗の法則:(xⁿ)′ = n·xⁿ⁻¹。積の法則:(uv)′ = u′v + uv′。連鎖律:(f(g(x)))′ = f′(g(x))·g′(x)。f′(x) = 0 のところではグラフの接線が水平になり、極大・極小・鞍点のいずれかです。
答え: b = 6
微分:導関数 f′(x) は x におけるグラフの傾きです。べき乗の法則:(xⁿ)′ = n·xⁿ⁻¹。積の法則:(uv)′ = u′v + uv′。連鎖律:(f(g(x)))′ = f′(g(x))·g′(x)。f′(x) = 0 のところではグラフの接線が水平になり、極大・極小・鞍点のいずれかです。
答え: |E| = 0.02963
数値積分:[a, b] を、幅 h = (b − a)/n、節点 x₀ … xₙ の n 個の短冊に分けます。台形公式:T = (h/2)[f₀ + 2f₁ + 2f₂ + … + 2fn−1 + fₙ]。シンプソンの公式(n は偶数):S = (h/3)[f₀ + 4f₁ + 2f₂ + 4f₃ + … + 4fn−1 + fₙ]。台形公式の誤差は h² のように小さくなり(M = max|f″| として評価式 (b − a)h²·M/12)、シンプソンの公式は h⁴ のように小さくなるので、同じ節点ならシンプソンの公式のほうが通常はずっと正確です。
相対誤差:2.092%。
答え: y3 = 1.251
y′ = f(t, y) の数値解法:オイラー法は、各ステップの開始点での傾きに沿って進みます:yn+1 = yₙ + h·f(tₙ, yₙ)。ホイン法(改良オイラー法)は、開始点での傾きと、オイラー法による予測点での傾きを平均します:k₁ = f(tₙ, yₙ)、k₂ = f(tₙ + h, yₙ + h·k₁)、yn+1 = yₙ + (h/2)(k₁ + k₂)。古典的ルンゲ–クッタ法(RK4)は 4 つの傾きを組み合わせます:yn+1 = yₙ + (h/6)(k₁ + 2k₂ + 2k₃ + k₄)。大域誤差は、オイラー法では h に、ホイン法では h² に、RK4 では h⁴ に比例します。
ExamProblems · examproblems.com/courses/differential-equations