0
$\begingroup$

I have a n x n matrix. I would like to

(a) take successively higher powers of the matrix and then multiply by projection vectors until the resulting vectors differ by only a scalar factor.

(b) calculate the dominant eigenvalue of the matrix to compare to (a)

and (c) calculate, using the same tactic as in (a), the dominant right and left eigenvectors.

This is too much work to do by hand, so my question is: can anyone recommend a program/language or package that would be ideal for the above calculations?

Thanks.

  • 0
    It seems to me that OP wants nothing more than a simple *power method* implementation, and the answers given thus far do quite a bit more than that... [NAPACK](http://www.netlib.org/napack/) has [a power method implementation](http://www.netlib.org/napack/power.f) (in FORTRAN of course); translation to other computing environments shouldn't be too hard...2012-04-25

3 Answers 3