Inverse Iteration for Finding Eigenvalues

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