Convergence Analysis of Legendre Pseudospectral Scheme for Solving Nonlinear Systems of Volterra Integral Equations

We are concerned with the extension of a Legendre spectral method to the numerical solution of nonlinear systems of Volterra integral equations of the second kind. It is proved theoretically that the proposed method converges exponentially provided that the solution is sufficiently smooth. Also, three biological systems which are known as the systems of Lotka-Volterra equations are approximately solved by the presented method. Numerical results confirm the theoretical prediction of the exponential rate of convergence.