De Vogelaere's method for the automatic solution of systems of coupled homogeneous second-order differential equations (Q796271)

This paper is concerned with the numerical solution of systems of coupled second order boundary value problems arising in quantum mechanical scattering problems using an automatic implementation of de Vogelaere's fourth order method. Accurate starting values are provided and a strategy for choosing the stepsize based on the magnitude of the estimated local error is developed. The automatic version is tested on a 2 channel test problem and favourable comparisons are made with a fixed stepsize code.