Contd5 Contd5_r Contd8 Contd8_r Contx1 Contx1_r Contx2 Contx2_r Derkf Derkf_r DerkfInt Dop853 Dop853_r Dopri5 Dopri5_r Doprin Doprin_r Dverk Dverk_r DverkInt Odex Odex2 Odex2_r Odex_r Retard Retard_r Ylag Ylag_r

Contx2_r() Function Contx2_r ( I As  Long, T As  Double, RCont() As  Double, ICont() As  Long  ) 常微分方程式の初期値問題 (補外法) (2階微分方程式用) (リバースコミュニケーション版) (密出力のための補間) 注 - 本プログラムは次バージョンで廃止予定です. 目的 本ルーチンは, Odex2_r を使って2階の常微分方程式の初期値問題の解を求める際に, 密出力を行うための補間を行う補助ルーチンである. Odex2_r を用いて密出力を行うためには, Iout = 1 として Odex2_r を呼び出す必要がある. ステップが成功すると IRev = 5 で Odex2_r から戻るので, 本ルーチンを使って直近のステップの区間内の任意の T における解の補間値を求めることができる. 戻り値 Double 解の補間値 Y(I). 引数 [in] I 求める解の要素番号. (0 <= I <= N - 1) [in] T 補間値を求める T の値. T は直近のステップの区間内になければならない. すなわち, 前回 IRev = 5 で Odex2_r から戻ったときの T を T1, 今回を T2 とすると, T1 <= T <= T2 となっている必要がある. [in] RCont() 配列 RCont(LRCont - 1) (LRCont >= (2*Km + 6)*N) 密出力のための制御情報. [in] ICont() 配列 ICont(LICont - 1) (LICont >= N) 密出力のための整数制御情報. 出典 E. Hairer, S.P. Norsett and G. Wanner, "Solving Ordinary Differential Equations. Nonstiff Problems. 2nd edition", Springer Series in Computational Mathematics, Springer-Verlag (1993) 邦訳: 「常微分方程式の数値解法Ⅰ 基礎編」スプリンガージャパン (2007) 使用例 Odex2_r の使用例(2)を参照せよ.