Numerical solution of ODE