Contex Contex_r Contr5 Contr5_r Contra Contra_r Contro Contro_r Contrp Contrp_r Dassl Dassl_r Radau Radau5 Radau5_r Radau_r Radaup Radaup_r Rodas Rodas_r Seulex Seulex_r

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