In this paper, we present and analyze a single interval
Legendre-Gauss spectral collocation method for solving the second
order nonlinear delay differential equations with variable delays. We also propose a novel algorithm for the single interval
scheme and apply it to the multiple interval scheme for more
efficient implementation. Numerical examples are provided to
illustrate the high accuracy of the proposed methods.