Contd5 Contd5_r Contd8 Contd8_r Contx1 Contx1_r Contx2 Contx2_r Deabm Deabm_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

Contd5_r() Function Contd5_r ( I As  Long, T As  Double, RCont() As  Double, ICont() As  Long  ) Initial value problem of ordinary differential equations (5(4)-th order Dorman-Prince method) (Reverse communication version) (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 Dopri5_r. Dopri5_r must be called with Iout = 1 to enable dense output. Additional function evaluation for dense output is not required. When returned from Dopri5_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 Dopri5_r with IRev = 5, and T2 is the value on this time. [in] RCont() Array RCont(LRCont - 1) (LRCont >= 5*N) Control information for dense output. [in] ICont() Array ICont(LICont - 1) (LICont >= N) Integer control information for dense output. Reference E. Hairer, S.P. Norsett and G. Wanner, "Solving Ordinary Differential Equations. Nonstiff Problems. 2nd edition", Springer Series in Computational Mathematics, Springer-Verlag (1993) Example Program See example program (2) of Dopri5_r.