Performs integration over a single time step using a built-in ODE solver. At present, a single solver is implement with limited options. The interface of this method may change when support for other solvers is added.
A named numeric vector holding the values of state variables and process rates in all boxes.
Numeric. Initial time.
Numeric. Length of time step of interest.
Minimum tolerated internal step size. The default of NULL
sets this to 10 times the value of .Machine$double.eps.
Maximum tolerated number of sub-steps.
Numeric. Relative accuracy requested. This is currently a global value, i.e. one cannot set the accuracy per state variable.
String. Currently, 'rk5' is the only method implemented. This is a Runge-Kutta Cash-Karp solver adapted from Press et al. (2002), Numerical recipes in Fortran 90. It is designed to handle non-stiff problems only.
Logical. Can be used to avoid repeated checks of arguments. This may increase performance in repeated calls.
Use deSolve for advanced solvers with more
options and capabilities to handle stiff problems.