This article presents a novel continuous numerical method designed for the numerical integration of general third-order initial value problems (IVPs) of ordinary differential equations (ODEs). A combination of power series and exponential function was formulated for the purpose of collocation and interpolation at nodal and off-nodal points to generate system of linear equations necessary for the method. The resulting hybrid linear multistep method was implemented using block mode approach. Consistency, stability, and convergence of the method were verified using established criteria. The developed method was applied directly to solve linear and nonlinear third order ODEs without reducing them to systems of first-order equations. Computational results demonstrated better accuracy when compared with existing numerical methods in the literature.
Read full abstract