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

Contra_r() Function Contra_r ( I As  Long, T As  Double, Cont() As  Double  ) Initial value problem of ordinary differential equations (Variable (5, 9, 13-th) order implicit Runge-Kutta method (Radau IIA)) (Interpolation for dense output) NOTE - THIS PROGRAM IS DEPRECATED AND WILL BE REMOVED IN THE NEXT VERSION. 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 Radau_r. Radau_r must be called with setting Iout = 1 to enable dense output. When returned from Radau_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 Radau_r with IRev = 5, and T2 is the value on this time. [in] Cont() Array Cont(LCont - 1) (LCont >= (Nsmax + 1)*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 Radau_r.