R 行列

• A（行列）のn乗

n<-3
A<-matrix(c(1:(n^2)),n,n)
x<-eigen(A)                #固有値と固有ベクトル
D<-diag(x$value)         #対角行列
P<-x$vectors             #固有ベクトルが行列になっている

A^3                       #これは要素を３乗した行列
A%*%A%*%A
P%*%D^3%*%solve(P)　　　　#対角化して計算

• • eigen(A):　$valueと$vectoresがある
  • diag(v):　対角行列
  • solve(A):　逆行列

• 上記記事の修正についてこちらを参照