Algorithm 689: Discretized collocation and iterated collocation for nonlinear Volterra integral equations of the second kind

This paper describes a FORTRAN code for calculating approximate solutions to systems of nonlinear Volterra integral equations of the second kind. The algorithm is based on polynomial spline collocation, with the possibility of combination with the corresponding iterated collocation. It exploits certain local superconvergence properties for the error estimation and the stepsize strategy.