dif_equ3

機能：

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

書式：

dif_equ3(F(x,y,y1,y2),x0,y0,y01,y02,h,[ox],[oy],x始,x終,x増)

解説：

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

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

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

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

　　　　　y01 ･･･････････････ １次微分dy/dx 初期値

　　　　　y02 ･･･････････････ ２次微分dy2/dx2 初期値

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

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

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

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

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