一階常微分方程式ソルバー

一階常微分方程式ソルバー

一階常微分方程式ソルバーは、dy/dx = f(x, y) 形式の常微分方程式を入力すると、その構造（変数分離型、線形、自律型、完全微分、または一般形式）を自動的に分類し、可能な場合は記号的な解析解を、不可能な場合でも高精度な数値解を生成します。解曲線が重ね書きされたライブ方向場（スロープフィールド）の可視化により、方程式の幾何学的な意味が即座に明らかになります。解曲線とは、すべての点において方向場の矢印に接する曲線そのものです。

一階常微分方程式とは？

一階常微分方程式は、未知関数 y(x) とその一階導関数 y'(x) のみを含む方程式です。標準的な明示的形式は次の通りです：

初期条件 y(x₀) = y₀ と組み合わせることで、これは初期値問題 (IVP) と定義されます。ピカール・リンデレフの定理により、f が (x₀, y₀) の近傍で y に関してリプシッツ連続であれば、一意の解が保証されます。幾何学的には、IVP は (x₀, y₀) を通り、すべての点での傾きがその点での f と一致する唯一の曲線（方向場に接する曲線）を求めることと同義です。

ソルバーが認識する6つのクラス

解析的手法：定数係数線形微分方程式

右辺が定数 a, b を用いて dy/dx = a·y + b と簡略化される場合、積分因子 μ(x) = e^(-a·x) を用いることで厳密解が得られます。一般解は以下の通りです：

初期条件 y(x₀) = y₀ を適用することで定数 C が決定され、唯一の特殊解が得られます。このクラスは、教科書に登場する非常に多くの問題をカバーしています：

• 指数関数的成長 — dy/dx = k·y、特殊解 y(t) = y₀·e^(k·t)。
• 指数関数的減衰 — dy/dx = -k·y、半減期 ln 2 / k。
• ニュートンの冷却法則 — dy/dx = -k·(y - T_amb)、物体温度は周囲温度に向かって指数関数的に収束します。
• RC回路の充電 — dV/dt = (1/RC)·(V_in - V)、コンデンサ電圧は電源電圧に近づきます。
• 薬物クリアランス — 消失速度定数 k を持つ一階薬物速度論。

方向場の読み方

各グリッド点 (x, y) において、ツールは傾きが f(x, y) に等しい短い線分を描きます。注目すべき3つのポイント：

• 平衡点は f(x, y) = 0 となる点で、方向場は水平になります。自律型方程式の場合、これらは f(y*) = 0 を満たす不動点 y* です。付近の軌道はここへ近づく（安定）か、遠ざかる（不安定）かの挙動を示します。
• 等傾斜線 (Isoclines) は f(x, y) が定数 c に等しくなる曲線で、この曲線上のすべての矢印は同じ傾き c を持ちます。
• 解曲線は決して交差しない（f がリプシッツ連続の場合）。視覚的にも、交差する2つの曲線が交点で異なる傾きを持つことは不可能なため、明らかです。

数値解法：古典的ルンゲ＝クッタ法 (RK4)

(x_n, y_n) が与えられたとき、次の値は4つの傾き推定値の加重平均によって計算されます：

RK4は局所打ち切り誤差が O(h⁵)、全域誤差が O(h⁴) であり、非スティフな方程式においてデフォルトのステップ数で概ね6桁の精度を提供します。ソルバーは初期点から x の両方向に積分を行い、y の絶対値が 10¹⁵ を超えた場合は（dy/dx = y² のように有限時間で爆発する解を想定して）適切に停止します。

この電卓の使い方

1. 右辺を入力: dy/dx = ... のフィールドに式を入力します。変数には x と y を使い、乗算には *、累乗には ^ または **、標準関数（sin, cos, exp, log, sqrt など）を使用します。定数 pi と e も認識されます。
2. 初期条件を指定: (x₀, y₀) を入力します。唯一の解曲線はこの点を必ず通過します。
3. x の範囲を選択: 方向場と解曲線を表示する範囲を決めます。y の範囲は計算された解から自動調整されます。
4. 「計算と可視化」をクリック: まず分類が実行されます。定数係数線形などの解析解パターンに一致すれば、記号解が表示されます。方向場と数値解曲線は常に描画されます。
5. 方向場の切り替え: 解曲線に集中するために方向場を非表示にしたり、アニメーションを再再生して積分がどのように進むかを確認したりできます。

計算例：ニュートンの冷却法則

80 °C のコーヒーが 20 °C の部屋で冷める様子を考えます。熱伝達率は温度差に比例します：

これは定数係数線形微分方程式 (a = -0.1, b = 2) です。解析解は：

30分後: T(30) = 20 + 60·e⁻³ ≈ 22.99 °C。方向場を見ると、開始温度に関わらず、すべての解曲線が水平線 T = 20 に漸近していく様子が明確に分かります。

一般的な応用例

• 人口動態 — 指数関数、ロジスティック、アリー効果モデル。
• 薬物速度論 — 薬物の吸収と消失、半減期の計算。
• 熱伝達 — ニュートンの冷却法則、集中熱容量モデル。
• RCおよびRL回路 — 一階線形電気回路の過渡現象。
• 放射性崩壊 — 単一核種の崩壊鎖。
• 混合槽問題 — 流入・流出がある槽内の溶質濃度。
• 空気抵抗がある落下体 — 終端速度分析 dv/dt = g - kv。

よくある質問

一階常微分方程式とは何ですか？

一階常微分方程式（ODE）とは、未知関数 y(x) とその一階導関数を含む dy/dx = f(x, y) という形式の方程式です。常微分方程式を解くとは、その微分が右辺と一致する関数 y(x) を見つけることを意味します。初期条件 y(x₀) = y₀ が与えられれば、穏やかな正則性の仮定（ピカール・リンデレフの定理）の下で解は一意に定まります。

方向場（スロープフィールド）とは何ですか？

方向場（または勾配場）は、各グリッド点 (x, y) において、その傾きが f(x, y) に等しい小さな線分をプロットしたものです。常微分方程式の解曲線は、すべての点においてこれらの線分に接する曲線となります。方向場は、方程式を記号的に解くことなく、解の全体的な振る舞いを瞬時に視覚的に理解させてくれます。

このツールはどの種類の一階常微分方程式を解くことができますか？

ツールは方程式を自動的に次のいずれかに分類します：積分可能（x のみに依存）、定数係数線形 y' = a·y + b（解析解を提供）、自律型（y のみに依存）、変数分離型、変数係数線形、または一般形式。すべてのクラスにおいて、高精度なルンゲ＝クッタ数値解と方向場の可視化が生成されます。

どのような数値解法が使用されていますか？

古典的な四次ルンゲ＝クッタ法（RK4）が、初期点から各方向に300のサブステップで適用されます。RK4は局所打ち切り誤差が O(h⁵) であり、この規模の非スティフな常微分方程式における標準的な手法です。ソルバーは発散（オーバーフローまたは NaN）を検出し、適切に停止します。

線形常微分方程式の積分因子法とは何ですか？

線形一階常微分方程式 y' + P(x)·y = Q(x) の場合、両辺に積分因子 μ(x) = e^∫P(x) dx を掛けます。左辺は完全微分 d/dx[μ·y] となるため、積分することで解が得られます。P と Q が定数の場合、これは解析解 y = -b/a + C·e^(a·x) に帰着します。

このツールはスティフな方程式や連立常微分方程式を扱えますか？

このソルバーは非スティフな一階スカラー常微分方程式を対象としています。非常にスティフな問題には後退オイラー法などの陽的な手法が、結合系にはベクトル値ソルバーが必要です。そのような場合は、SciPy の solve_ivp などの専門的なツールを使用してください。

参考文献

• 常微分方程式 — Wikipedia
• 方向場 — Wikipedia (英語)
• ルンゲ＝クッタ法 — Wikipedia
• 積分因子 — Wikipedia
• ピカール・リンデレフの定理 — Wikipedia