A linearized difference scheme for a class of fractional partial differential equations with delay

A class of non linear fractional partial differential equations with initial and Dirichlet boundary conditions is under consideration. We seek to obtain numerical solutions for this considered class of equations based on finite difference method. The convergence order will be $2-\alpha$ in time and four in space. A numerical example is given to support the theoretical results.