This example computes all the eigenvalues of the matrix
, where
and those eigenvectors which correspond to eigenvalues
such that
. Here
is general and must first be reduced to upper Hessenberg form
by
zgehrd. The program then calls
zhseqr to compute the eigenvalues, and
zhsein to compute the required eigenvectors of
by inverse iteration. Finally
zunmhr is called to transform the eigenvectors of
back to eigenvectors of the original matrix
.