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

Abstract

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