Crank-Nicolson and Legendre spectral collocation methods for a partial integro-differential equation with a singular kernel

In this paper we present an efficient numerical method for the solution of a partial integro-differential equation with a singular kernel. In the time direction, a Crank-Nicolson finite difference scheme is used to approximate the differential term and the product trapezoidal method is employed to treat the integral term. Also for space discretization we apply Legendre spectral collocation method. We discuss the stability and convergence of proposed method and show that the method is unconditionally stable and convergent with order O(τ32+N−s) where τ, N and s are time step size, number of collocation points and regularity of exact solution respectively. We compare the numerical results of proposed method with the results of other schemes in the literature in terms of accuracy, computational order and CPU time to show the efficiency and applicability of it.