ROCK4

explicit solver for stiff equations

This is a MATLAB-implementation of ROCK4 (Abdulle & Medovikov)for the numerical solution of a (mildly) stiff system of first order differential equations ODEs y'=f(t,y).

ROCK4 is based on a family of second order explicit Runge-Kutta methods with nearly optimal stability domain on the negative real axis.

Usage is the same as ode23s:
[TOUT,YOUT] = rock(ODEFUN,TSPAN,Y0)
or
[TOUT,YOUT] = rock4(DEFUN,TSPAN,Y0,OPTIONS)
or
[TOUT,YOUT,TE,YE,IE] = rock4(ODEFUN,TSPAN,Y0,OPTIONS)