Numerical method with high order accuracy for solving a anomalous subdiffusion equation

In this paper, a numerical method with second order temporal accuracy and fourth order spatial accuracy is developed to solve a anomalous subdiffusion equation; by Fourier analysis, the convergence, stability and solvability of the numerical method are analyzed; the theoretical results are strongly supported by the numerical experiment.