行列の指数関数というのがある 簡単に言うと、固有値を指数の肩に乗せて、それを固有ベクトルの行列でサンドイッチする exp.m <- function(A,n=1,eigen=FALSE){ # 固有値分解 eigen.out<-eigen(A) # P=V,P^{-1}=U V<-eigen.out[[2]] U<-solve(V) B<-diag(exp…

複素行列のdeterminant