LR algorithm with Laguerre shift for symmetric tridiagonal matrices. - J. Grad; E. Zakrajsek (Comput. J., 1972)