On the numerical solution of singular integro-differential equations