Classical Runge-Kutta of order 4.
rk4sys(f, a, b, y0, n, Pars, ...)
List with components x for grid points between a and b and y an n-by-m matrix with solutions for variables in columns, i.e. each row contains one time stamp.
function in the differential equation \(y' = f(x, y)\); defined as a function \(R \times R^m \rightarrow R^m\), where \(m\) is the number of equations.
starting time for the interval to integrate
ending time for the interval to integrate.
starting values at time a
the number of steps from a to b.
named list of parameters passed to f
additional inputs to function f
Classical Runge-Kutta of order 4 for (systems of) ordinary differential equations with fixed step size. Copied from pracma under GPL-3, with small modifications to work with RTMB