A new twelfth-order four-step formula containing fourth derivatives for the numerical integration of the one-dimensional Schrödinger equation has been developed. It was found that by adding multi-derivative terms, the stability of a linear multi-step method can be improved and the interval of periodicity of this new method is larger than that of the Numerov's method. The numerical test shows that the new method is superior to the previous lower orders in both accuracy and efficiency and it is specially applied to the problem when an increasing accuracy is requested.
Explicit two-step high-accuracy hybrid methods with minimal phase-lag for y″ = f(x, y) and their application to the one-dimensional Schrödinger equation
