適当に行列を作ってやってみる my.matrix.eigen <- function(lambdas,vs){ Vs <- t(t(vs) * lambdas) Vs %*% solve(vs) } exp.m <- function(A,n){ # 固有値分解 eigen.out<-eigen(A) # P=V,P^{-1}=U V<-eigen.out[[2]] U<-solve(V) B<-diag(exp(eigen.out[[…

２変量の常微分方程式があって、２変数の時間変化が２つのベクトルを軸としてその２軸のそれぞれに指数関数の係数を与えた和で表されるとき、そのy=1平面への射影に複比保存が表れるのだが、それの「証明」というか、ひたすらな式変形で納得するためのメモ …