“…Because of the non-local properties of fractional operators, obtaining the analytical solutions of the fractional differential equations (FDEs) is more challenging or sometimes even impossible. Hence the proposal, development, and analysis of numerical methods to solve FDEs are at present a quite active field of research, and many methods have been considered, for instance, finite difference method [2][3][4], finite element method [5,6], spectral method [7][8][9], meshless method [10][11][12][13][14], and so on.…”