Inverse iteration (source code)

= Inverse iteration
{wiki}

Inverse iteration repeatedly solves $(A-sI)y_{k+1}=q_k$ and normalizes. It is the power method for $(A-sI)^{-1}$ and converges toward an eigenvector whose eigenvalue is closest to the shift $s$.