“…Several numerical approaches have been developed and implemented to solve BVP (1.1)- (1.3). In this regard, these approaches include the Adomian decomposition method [4,5], the Laplace transform-homotopy perturbation [6], the homotopy analysis method [7], the finite difference method [8,9], Diagonal block method [10,11], B-spline collocation method [12,13], Cubic Hermite collocation method [14], Gegenbaur integration matrices [15] and the spectral method [16,17].…”