“…It is well known that the analytical solutions to the fractional differential equations are usually difficult to derive and (if luckily obtained) always contain some infinite series which make evaluation very expensive. Therefore, we resort to some numerical methods for fractional differential equations such as finite difference method [8][9][10][11], finite element method [12,13], spectral method [14], radial basis function collocation method [15], spectral Tau algorithm [16], a hybrid of lagrange operational matrix and Tau-collocation method [17] and other analytical methods (e.g., variational iteration method [18], homotopy perturbation method [19], a domain decomposition method [20]).…”