オイラー法（1階常微分方程式）
dy/dx = f(x, y) と初期条件・刻み幅を入力すると、オイラー法で数値解を計算。終点での y の値、各ステップの推移表、解曲線のグラフをまとめて表示します。
入力
右辺の例（タップで入力）
使える記号：+ - * / ^（べき乗）／関数：sin cos tan asin acos atan sinh cosh tanh exp log(自然対数) ln log10 sqrt cbrt abs pow(a,b)／定数：pi（π）, e／変数は x と y。
計算結果
dy/dx = y、 y(0) = 1 のとき
y(1) ≒
2.593742
終点の y
2.593742
x = 1
ステップ数
10 回
h = 0.1
初期条件
y(0) = 1
区間 0 → 1
数値解の推移
| n | x | y | 傾き f(x, y) |
|---|---|---|---|
| 0 | 0 | 1 | 1 |
| 1 | 0.1 | 1.1 | 1.1 |
| 2 | 0.2 | 1.21 | 1.21 |
| 3 | 0.3 | 1.331 | 1.331 |
| 4 | 0.4 | 1.4641 | 1.4641 |
| 5 | 0.5 | 1.61051 | 1.61051 |
| 6 | 0.6 | 1.771561 | 1.771561 |
| 7 | 0.7 | 1.948717 | 1.948717 |
| 8 | 0.8 | 2.143589 | 2.143589 |
| 9 | 0.9 | 2.357948 | 2.357948 |
| 10 | 1 | 2.593742 | 2.593742 |
計算方法・使い方
- 右辺 f(x, y)、初期値 x0・y0、刻み幅 h、終了 x（xend）を入力すると、オイラー法 y_{n+1} = y_n + h·f(x_n, y_n) によって x0 から xend まで順番に y を更新し、各ステップの値を求めます。
- 右辺の数式では変数 x と y の両方を使えます。記号は + - * / と ^（べき乗）、関数は sin・cos・tan・exp・log（自然対数）・ln・sqrt・abs・pow(a,b) など、定数は pi（π）と e に対応しています。三角関数の角度はラジアンで扱います。
- 主要な結果として終了 x における y の近似値を大きく表示し、あわせてステップ数・刻み幅・初期条件を確認できます。下の表ではステップ番号 n ごとの x・y・傾き f(x, y) の推移を、グラフでは解曲線と各点をたどれます。
- オイラー法は1次精度の最も基本的な解法で、刻み幅 h に比例した誤差が残ります。h を半分にすると誤差もおよそ半分になりますが、計算ステップ数は増えます。精度が必要な場合は h を十分小さくしてください。
- 刻み幅が大きすぎると解が振動したり発散したりすることがあります（特に -2y のように急減衰する式で顕著です）。値が発散した場合は警告を表示しますので、h を小さくするか区間を見直してください。
- 入力した式や計算結果は外部に送信されず、すべてブラウザ内で処理されます。安心してご利用ください。
関連する計算ツール
オイラー法（2階常微分方程式）
数値解析y''=f(x, y, y') の形の2階常微分方程式を、連立1階に変換してオイラー法で数値的に解くツール。初期条件と刻み幅を入れるだけで、終点での y の近似値・解のグラフ・各ステップの (x, y, y') 表を表示します。
計算する →ルンゲ・クッタ法（2次・1階常微分方程式）
数値解析dy/dx=f(x,y) を2次ルンゲ・クッタ法（中点法）で数値的に解き、終点での近似解・数値解の表・解曲線グラフを表示します。
計算する →ルンゲ・クッタ法（4次・1階常微分方程式）
数値解析dy/dx=f(x,y) の右辺の数式と初期値 x0・y0、刻み幅 h、終点 xEnd を入力すると、古典的4次ルンゲ・クッタ法（RK4）で常微分方程式を数値的に解きます。終点での近似解 y、各ステップの k1〜k4 と y の値の表、解曲線のグラフを表示します。
計算する →ルンゲ・クッタ法（2次・2階常微分方程式）
数値解析y''=f(x,y,y') を1階連立に変換し、2次ルンゲ・クッタ法で数値的に解きます。終点での y と y' を求め、解曲線のグラフと各ステップの数表を表示します。
計算する →
お客様の声
このツールを使った感想をお聞かせください。
レビューを投稿する
- ホーム
オイラー法（1階常微分方程式）