Hybrid Numerical Method with Block Extension for Direct Solution of Third Order Ordinary Differential Equations

This paper focuses on the development of a hybrid method with block extension for direct solution of initial value problems (IVPs) of general third-order ordinary differential equations. Power series was used as the basis function for the solution of the IVP. An approximate solution from the basis function was interpolated at some selected off-grid points while the third derivative of the approximate solution was collocated at all grid and off-grid points to generate a system of linear equations for the determination of the unknown parameters. The derived method was tested for consistency, zero stability, convergence and absolute stability. The method was implemented with five test problems including the Genesio equation to confirm its accuracy and usability. The rate of convergence (ROC) reveals that the method is consistent with the theoretical order of the proposed method. Comparison of the results with some existing methods shows the superiority of the accuracy of the method.