In this paper, we used an interpolation function to derive a Numerical Integrator that can be used for solving first order Initial Value Problems in Ordinary Differential Equation. The numerical quality of the Integrator has been analyzed to authenticate the reliability of the new method. The numerical test showed that the finite difference methods developed possess the same monotonic properties with the analytic solution of the sampled Initial Value Problems.