Abstract
Fractional differential equations describe nature adequately because of the symmetry properties which describe physical and biological processes. In this article, a fourth-order new implicit difference scheme is formulated and applied to solve the two-dimensional time-fractional modified sub-diffusion equation involving two times Riemann-Liouville fractional derivatives. The stability of the fourth-order implicit difference scheme is investigated using the von Neumann technique. The proposed scheme is shown to be unconditionally stable. Numerical examples are given to illustrate the feasibility of the proposed scheme.