R/MatrixExp.eig.R

Defines functions MatrixExp.eig

MatrixExp.eig <-
function(Q)
  {
    tmp <- eigen(Q, symmetric = FALSE)
    P1 <- tmp$vectors %*% diag(exp(tmp$values)) %*% solve(tmp$vectors)
    return(P1)
    
  }
lukejharmon/netphy documentation built on May 21, 2019, 8:58 a.m.