“…The numerical solutions on the IVPs and BVPs of high-order ODEs and SODEs have attracted the attention of many scholars. Most of these investigations are focused on two main aspects: one is to simplify the higher order ODE into an equivalent system of first-order ODEs, which involves a lot of computer time and more human effort and the other is to solve directly high-order ODEs, such as block methods (Majid et al, 2012), Runge-Kutta (R-K) methods (Hussain et al, 2016(Hussain et al, , 2017You and Chen, 2013), finite difference methods (Chawla and Katti, 1979), classical orthogonal polynomials (Legendre polynomials, Chebyshev polynomials and Jacobi polynomials), wavelet methods (Haar wavelets (Amin et al, 2021;Islam et al, 2010), Shannon wavelets (Shi and Li, 2014)), Adomian decomposition method with Green's function (Al-Hayani, 2011), spectral algorithms (Doha and Abd-Elhameed, 2009;Doha and Bhrawy, 2002;Doha et al, 2012) etc.…”