Development of a new numerical scheme for the solution of exponential growth and decay models
This paper presents the development of a new numerical scheme for the solution of exponential growth and decay models emanated from biological sciences. The scheme has been derived via the combination of two interpolants namely, polynomial and exponential functions. The analysis of the local truncation error of the derived scheme is investigated by means of the Taylor’s series expansion. In order to test the performance of the scheme in terms of accuracy in the context of the exact solution, four biological models were solved numerically. The absolute error has been computed successfully at each mesh point of the integration interval under consideration. The numerical results generated via the scheme agree with the exact solution and with the fifth order convergence based upon the analysis carried out. Hence, the scheme is found to be of order five, accurate and is a good approach to be included in the class of linear explicit numerical methods for the solution of initial value problems in ordinary differential equations.