-- A shift is needed to make sure the QR algorithm
-- converges. Since the roots (and thus eigenvalues) of orthogonal
-- polynomials are real and distinct, there does exist a shift
-- that will make it work. We guess lambda=1 makes it work.
-- A shift is needed to make sure the QR algorithm
-- converges. Since the roots (and thus eigenvalues) of orthogonal
-- polynomials are real and distinct, there does exist a shift
-- that will make it work. We guess lambda=1 makes it work.