a=matrix(c(1,2,3,5),2)
> a
     [,1] [,2]
[1,]    1    3
[2,]    2    5

b=eigen(a)
> b
$values
[1]  6.1622777 -0.1622777
$vectors
           [,1]       [,2]
[1,] -0.5024547 -0.9324647
[2,] -0.8646035  0.3612610
> lam=b$values
> lam
[1]  6.1622777 -0.1622777
> vec=b$vectors
> vec
           [,1]       [,2]
[1,] -0.5024547 -0.9324647
[2,] -0.8646035  0.3612610
>
n次方就是

t(lam)***vec^n***lam ##矩阵相乘，忘了哪个符号了。