For a nonzero real vector , the Rayleigh quotient of the matrix is
The Rayleigh quotient iteration starts from a unit vector . Given , compute the shift
solve the linear system
and normalize in the Euclidean norm:
The next eigenvalue estimate is . The process stops when the eigenpair residual is sufficiently small.