derkf derkf_int derkf_r dop853 dop853_r dopri5 dopri5_r doprin doprin_r dverk dverk_int dverk_r odex odex2 odex2_r odex_r retard retard_r

dverk_int() void dverk_int ( int  n, double  t, double  y[], double  c[], double  work[]  ) Initial value problem of ordinary differential equations (6(5)-th order Runge-Kutta-Verner method) (Interpolation for dense output) NOTE - THIS ROUTINE 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 dverk. For the details of interpolants used here, refer to the reference below. dverk must be called with setting parameter c[8] = 2 to enable dense output. Note that larger work array work[] is required. Further, one additional function evaluation per accepted step is required. When returned from dverk with ind = 5 (ind = 3 for last step) after accepted step, the interporated solution within last step interval can be computed by using this routine. Parameters [in] n Number of differential equations. (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 dverk with ind = 5, and t2 is the value on the last return with info = 5. [out] y[] Array y[ly] (ly >= n) Interpoolated solution at t. [in] c[] Array c[lc] Parameter array c[] used by dverk. [in] work[] Array work[lwork] Work array work[] used by dverk. Reference W H Enright et al. "Interpolants for Runge-Kutta Formulas" ACM Transactions on Mathematical Software Vol.12, No.3, 1986, pp.193-218