“…They are increasingly used in such fields as fluid flow, control theory of dynamical systems, diffusive transport akin to diffusion, probability, and statistics [3,4]. Moreover, since most fractional differential equations do not have exact analytical solutions, approximate and numerical techniques are used extensively to treat the fractional differential models, including, for example, the homotopy analysis method [5], homotopy perturbation approach [6], variational iteration technique [7], Chebyshev spectral method [8], orthogonal polynomial method [9], Grunwald-Letnikov method [10], fractional Adams method [11], and several other methods; read [12,13] and the given references therewith.…”