“…Many numerical solution techniques were developed to obtain analytical, semianalytical, and/or numerical solutions of fractional dynamical systems, for example, the finite difference method [8,9], predictor-corrector approach [10,11], operational matrix method [12,13], variational iteration method [14,15], homotopy perturbation method [16,17], Adomian's decomposition method [18,19], to mention a few. Due to the nonlocal character of the fractional derivative, storing the past responses requires a large amount of computer memory.…”