Numerical solution of initial value problems