Numerical simulation with high order accuracy for the time fractional reaction-subdiffusion equation

In this paper, we propose the numerical simulation method with second order temporal accuracy and fourth order spatial accuracy for the time fractional reaction-subdiffusion equation; the stability, convergence and solvability of the numerical simulation method respectively are discussed by Fourier analysis and algebraic theory; the theoretical analysis results very consistent with the numerical experiment.