The main purpose of this paper is to propose the Legendre spectral-collocation
method to solve the Volterra integral differential equations with nonvanishing delay
which arise in many problems, such as modeling in biosciences and population.
In our method we divide the definition domain of the solution into several
subintervals where the solution is sufficiently smooth. Then we can use the
spectral-collocation method for these equations in each subinterval. We provide
convergence analysis for this method, which shows that the numerical errors decay
exponentially. Numerical examples are presented to confirm these theoretical
results.