rk4sys function

Classical Runge-Kutta for system of equations