Abstract
An estimate for the fractional reaction subdiffusion equation is presented by a discrete Crank-Nicolson finite element method (FEM), which we can obtain by using the finite difference method (FDM) (in time) and the finite element method (in space). The proposed scheme is obtained at time t(n+1/2) because at this time there are some different coefficients compared to those at time t(n+1), that is, ((k+1/2)(1-beta) - (k-1/2)(1-beta)), instead of ((k+1)(1-beta)-(k)(1-beta)). We studied the stability analysis, truncation error, and convergence analysis of the derived scheme. Numerical examples are presented to illustrate the efficiency of the given algorithm.