Spectral Lanczos’ tau method for systems of nonlinear integro-differential equations