New P-stable high-order methods with minimal phase-lag for the numerical integration of the radial Schrodinger equation

A new family of P-stable two-step high-order methods with minimal phase-lag are developed for the numerical integration of the special second-order initial value problem. An application to the one-dimensional Schrodinger equation, indicates that these new methods are generally more accurate than other previously developed finite difference methods.