A Fourier method for the fractional diffusion equation describing sub-diffusion

In this paper, a fractional partial differential equation (FPDE) describing sub-diffusion is considered. An implicit difference approximation scheme (IDAS) for solving a FPDE is presented. We propose a Fourier method for analyzing the stability and convergence of the IDAS, derive the global accuracy of the IDAS, and discuss the solvability. Finally, numerical examples are given to compare with the exact solution for the order of convergence, and simulate the fractional dynamical systems.