Solves numerically an initial conditions problem for a set of Ordinary Differential Equations (ODE) by Runge-Kutta 4 method. Integrates numerically the equations from tmin to tmax by steps of dt.
multirunge(y0, tmin, tmax, dt, param, sys)Numerical vector : initial conditions (as many elements as equations in sys)
Minimum value of the independent variable (generally the time)
Maximum value of the independent variable
Time increment (default = 0.01)
Numerical vector providing the parameters for sys
The set of functions giving the right sides of the ODEs
A numerical vector of nfunct elements. Nfunct is the number of unknown functions of sys determined by the function as length(sys())
sys must be provided by the user. Please edit the object sys (see ?sys and examples)