Finite-difference schemes for a diffusion equation with fractional derivatives in a multidimensional domain

In a multidimensional domain, we consider a partial differential equation with fractional space and time derivatives. For the first initial-boundary value problem, we consider a purely implicit scheme based on the approximate factorization method. We prove the stability of the implicit scheme for the considered class of problems.