dif_equ1

機能：

任意１階常微分方程式の解を求む（初期値問題）

書式：

dif_equ1(F(x,y),x0,y0,h,[ox],[oy],x始,x終,x増)

解説：

dy/ dx = F(x,y)　初期条件：x＝x0時､y＝y0 ４次のルンゲクッタ法（進み幅h）を使い､上記の１階微分方程式を解き､指定範囲[x始,x終]と刻み値での数値解を数値組[ox],[oy]に出力します。 x初＜x0の場合でも､正しい解が得られます。

　　　　　F(x,y) ･･･････････････ 関数式

　　　　　x0 ･･･････････････ Ｘ初期値

　　　　　y0 ･･･････････････ Ｙ初期値

　　　　　h ･･･････････････ ルンゲクッタ法の進み幅(h＞0:小さいほど精度は良いが､計算時間が長くなる)

　　　　　[ox],[oy] ･･･････････････数値解のＸＹ座標の出力数値組

　　　　　x始,x終 ･･･････････････数値解の出力区間

　　　　　x増 ･･･････････････正:刻み値､負整数:均分点数

例：dy/dx = 1/(1+x^2)-2*y^2 の解を求めて下さい。
厳密解　y = x / (1 + x^2)