Numerical Method for Ordinary Differential Equations