我一直在研究使用R的与特征值和特征向量相关的一些矩阵问题。例如：其中A是可对角化的3x3矩阵求解A ^ 740。
我发现R中的eigen（）对于快速求解特征值非常有用
A = matrix(c(1,-2,4,0,-1,0,0,0,-1), byrow=T, nrow=3, ncol=3)
eigen(A)$values
[1] 1 -1 -1
然而，当特征向量通过eigen（）中的向量长度归一化时，我发现自己手动解算它们（A-λI）= 0 这些形成了一个适当的解决方案：
P = (v1, v2, v3) #eigenvectors
PI = solve(P) #inverse
D=diag(3)*c(1,-1,-1)
A740 = P %*% D**740 %*% PI
我想知道是否有一种方法可以使用eigen（）并返回非标准化的特征向量，用于这样的简单应用。