戻る|進む

DCL:MATH２:ODELIB : 常微分方程式 (ルンゲクッタ):メソッド(元サブルーチン)の説明

3.3.4 ODRKDV

1.

機能

各ステッパルーチンにより, 要求精度を満たすように積分を行う.

2.

呼び出し方法

CALL ODRKDV(N,STEPER,FCN,T,TINT,DT,X,WORK)

3.

パラメーターの説明

n	(I)	被積分変数(方程式)の数. (i)
steper	手続き名	使用するスッテパルーチン名.
fcn	手続き名	dxを計算するメソッド(元サブルーチン)名 (i)
t	(R)	積分を始める独立変数tの値. RETURN 時には最後のステップt (TEND) が出力される. (i/o)
tend	(R)	積分を終了するtの値. (i)
dt	(I)	ステップ幅の初期値. (i/o)
x	R(N)	被積分変数の t=T における値を入力し, t=TEND における値を出力する. (i/o)
work	R(N,M)	作業変数. M=7 ( odrk4), M=5 (その他).

4.

備考

(a)

精度は odpget/odpset の 管理する内部変数 'EPSILON'によって指定できる. 最初の呼び出し時に, 引数 dt に適当な値を代入しておくと 要求精度を満たすように dt が設定されて返って来るので, これを, そのまま次回の呼び出し時に指定すればよい.

(b)

odrkdv は odpget/odpset の管理する内部変数 'MAXSTEP' で指定された 回数だけ積分しても tend に到達しないときは, エラーメッセージを出力して, 実行を中止する. 'MAXSTEP' は odiset で変更可能である.

(c)

使用するステッパルーチンは, 特に支障がない限り, odrkgr (Runge-Kutta-Gill 再計算) を推奨する.

(d)

fcn は メソッド(元サブルーチン)の形でユーザーが用意する. その形式は以下のとおり.

SUBROUTINE FCN(N,T,X,DX)