“…It was found that various applications can be modeled with the help of the fractional derivatives [6,7]. For example, the nonlinear oscillation of earthquake [8], the fractional optimal control problems for dynamic systems [9,10,11,12], and the fluid-dynamic models with fractional derivatives can eliminate the deficiency arising from the assumption of continuous traffic flow [13,14,15]. During the last decades, several methods have been used to solve fractional differential equations, fractional partial differential equations, fractional integro-differential equations, the initial and boundary value problems, and dynamic systems containing fractional derivatives, such as Adomian's decomposition method [16,17], fractional-order Legendre functions [18], fractional-order Chebyshev functions of the second kind [19], Homotopy analysis method [20], Bessel functions and spectral methods [21], Legendre and Bernstein polynomials [22], finite element methods [23], Legendre collocation [24], modified spline collocation [25], multiquadratic radial basis functions [26], and other methods [27,28,29,30,31,32,33].…”