High order Bessel fitting methods for the numerical integration of the Schrödinger equation

In this paper we show how to construct explicit multistep algorithms for an accurate and efficient numerical integration of the radial Schrödinger equation. The proposed methods are Bessel fitting, that is to say, they integrate exactly any linear combination of Bessel and Newman functions and ordinary polynomials. They are the first of the like methods that can achieve any order.