“…In many fractional-order models, the analytical solution is more complicated, therefore the numerical solutions for these models are appropriate [1,6,20]. These numerical so-lutions depend on several techniques such as finite difference, finite volume, variational iteration, Legendre polynomials, Chebyshev collocation, homotopy perturbation, operational matrix, variational iteration, Adams-Bashforth, nonstandard finite difference, sinccollocation, compact finite difference, tau method, block pulse, decomposition, radial basis, Taylor collocation, and wavelets spectral (see, for example, [5,11,24,25,35,38]).…”