摘要

We introduce a multistep Legendre-Gauss spectral collocation method for the nonlinear Volterra integral equations of the second kind. This method is easy to implement and possesses high order accuracy. In addition, it is very suitable for long time calculations. We also derive the optimal convergence of the hp-version of the multistep collocation method under the L-2-norm. Numerical experiments confirm the theoretical expectations.