Runge kutta 4th order method using mathematica