If 
is symmetric one can use an old, but efficient
method. Note that eigenvalues and eigenvectors of
are real.
The power method can be also rigged up to find a series of
eigenvalues: the idea is to remove the eigenvalues once these are
found. This is called ``DEFLATION'' and must be used with caution
since round off becomes ever more significant.
