Given some , we set:

This converges to an eigenvector with eigenvalue closest to We can make it faster by adjusting throughout the iteration Take to be the minimizer of

Note that we will need LU factorisation of This is very efficient if is Upper Heisenberg

We can bring in this form using Givens Rotation or Householder Reflection