Abstract
The solution of a nonlinear system of first order differential equations with nonlinear boundary conditions by implicit Runge-Kutta methods based on interpolatory quadrature formulae is examined. An equivalence between implicit Runge-Kutta and collocation schemes is established. It is shown that the difference equations obtained have a unique solution in a neighbourhood of an isolated solution of the continuous problem, that this solution can be computed by Newton iteration and that it converges to the isolated solution. The order of convergence is equal to the degree of precision of the related quadrature formula plus one. The efficient implementation of the methods is discussed and numerical examples are given.