Compact finite difference scheme for the solution of a time fractional partial integro‐differential equation with a weakly singular kernel

In this paper, we develop a high‐order finite difference scheme for the solution of a time fractional partial integro‐differential equation with a weakly singular kernel. The fractional derivative is used in the Riemann‐Liouville sense. We prove the unconditional stability and convergence of scheme using energy method and show that the convergence order is . We provide some numerical experiments to confirm the efficiency of suggested scheme. The results of numerical experiments are compared with analytical solutions to show the efficiency of proposed scheme. It is illustrated that the numerical results are in good agreement with theoretical ones.