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

Contr5_r() Function Contr5_r ( I As  Long, T As  Double, Cont() As  Double  ) Initial value problem of ordinary differential equations (5-th order implicit Runge-Kutta method (Radau IIA)) (Interpolation for dense output) Purpose This is the support routine which interpolates the solution for dense output when an initial value problem of a system of ordinary differential equations is solved by Radau5_r. Radau5_r must be called with setting Iout = 1 to enable dense output. When returned from Radau5_r with IRev = 5 after each accepted step, the interpolated solution at any T within latest step interval can be computed by using this routine. Returns Double Interpolated solution Y(I). Parameters [in] I Element number of solution to be interpolated. (0 <= I <= N - 1) [in] T T at which the interpolated solution is computed. T must be in the interval of the last step. That is, T1 <= T <= T2 where T1 is the value of T on the previous return from Radau5_r with IRev = 5, and T2 is the value on this time. [in] Cont() Array Cont(LCont - 1) (LCont >= 4*N) Control information for dense output. Reference 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) Example Program See example program (2) of Radau5_r.