This paper presents the development and analysis of a proposed scheme to solve Initial Value Problems (IVPs). The proposed scheme is devised by means of the interpolating function. The properties of the proposed scheme such as the local truncation error, order of accuracy, stability, consistency, and convergence are analyzed. Furthermore, the performance of the proposed scheme is tested on ?ve numerical examples. Moreover, the comparative study of the results generated via the proposed scheme and the exact solution is presented. Hence, the proposed scheme has ?fth order convergence and is a good tool for approximating the solution of IVPs.