Abstract
Fractional analog of the reaction diffusion equation is used to model the subdiffusion process. Diffusion equation with fractional Riemann–Liouville operator is analyzed in this paper. We offer finite-difference methods that can be used to solve the initial-boundary value problems for some time-fractional order differential equations. Stability and convergence theorems are proved.