Block boundary value methods for solving linear neutral Volterra integro-differential equations with weakly singular kernels

Abstract

A class of block boundary value methods is constructed for the solution of linear neutral Volterra integro-differential equations with weakly singular kernels. Under suitable conditions on the data, it is shown that these methods yield optimal convergence rates when implemented on special graded meshes. Furthermore, these methods are easily extended to solve linear Volterra integral equations of the 2nd kind with weakly singular kernels. Numerical experiments confirm the theoretical results and the accuracy of the methods, and a comparison with piecewise polynomial collocation methods is provided.