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)
odrkdvodpget/odpset の管理する内部変数 'MAXSTEP' で指定された 回数だけ積分しても tend に到達しないときは, エラーメッセージを出力して, 実行を中止する. 'MAXSTEP'odiset で変更可能である.
(c)
使用するステッパルーチンは, 特に支障がない限り, odrkgr (Runge-Kutta-Gill 再計算) を推奨する.
(d)
fcn は メソッド(元サブルーチン)の形でユーザーが用意する. その形式は以下のとおり.
SUBROUTINE FCN(N,T,X,DX)