“…Currently, methods for solving FDEs with initial conditions can be classified into two classes, namely, approximative method and analytical methods. Typical approximative methods include the operational matrix method based on orthogonal functions, the predictor-corrector method, fractional Euler method, and so on ( [11][12][13][14][15][16]). The most practical analytical methods are the Adomian decomposition method, the homotopy analysis method, the homotopy perturbation method, the Laplace transform method, and the variational iteration method ( [17][18][19][20][21]).…”