A symmetric matrix method for Schrodinger eigenstates
Dennis Dunn, Brian Grieves · Journal of Physics A Mathematical and General · 1989
It is shown how the one-dimensional Schrodinger eigenvalue equation can be transformed into a symmetric generalised matrix eigenproblem with a local truncation error of only delta 4 , where delta is the step size, as in the Numerov algorithm. Standard sparse matrix library packages are available for the solution and it is demonstrated that groups of eigenvalues and their associated eigenvectors (wavefunctions) can be determined simultaneously with high precision and speed.