“…Thus a variety of numerical and approximate methods have been developed to solve these, such as operational matrices, [19] homotopy perturbation method (HPM), [20] homotopy analysis method (HAM), [21] fractional differential transform method, [22] extrapolation method, [23] Legendre wavelets, [24] and Sinc approximation with the single exponential (SE). [25] Li [26] proposed a novel application of Taylor expansion method for approximate solution of linear ordinary differential equations with variable coefficients. The method expanded by Li and his co-authors to solve Abel's integral equation, [18,27] Riccati equation, [28] an integral equation with fixed singularity for a cruciform crack, [29] a class of linear integro-differential equations, [30] and fractional integro-differential equations.…”