戻る|進む

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

3.3.3 ODRKDU

1.

機能

各アルゴリズムルーチンにより, 固定刻幅で積分を行う.

2.

呼び出し方法

CALL ODRKDU(N,ALGOR,FCN,T,TEND,ISTEP,X,WORK)

3.

パラメーターの説明

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

4.

備考

(a)

使用するアルゴリズムルーチンは, 特に支障がない限り, odrkg (Runge-Kutta-Gill) を推奨する.

(b)

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

SUBROUTINE FCN(N,T,X,DX)