Learn R Programming

rodeo (version 0.8.2)

step: Numerical Integration Over a Single Time Step

Description

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.

Value

A named numeric vector holding the values of state variables and process rates in all boxes.

Arguments

t0

Numeric. Initial time.

h

Numeric. Length of time step of interest.

hmin

Minimum tolerated internal step size. The default of NULL sets this to 10 times the value of .Machine$double.eps.

maxsteps

Maximum tolerated number of sub-steps.

tol

Numeric. Relative accuracy requested. This is currently a global value, i.e. one cannot set the accuracy per state variable.

method

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.

check

Logical. Can be used to avoid repeated checks of arguments. This may increase performance in repeated calls.

See Also

Use deSolve for advanced solvers with more options and capabilities to handle stiff problems.