MatrixExp.eig <-
function(Q)
{
tmp <- eigen(Q, symmetric = FALSE)
P1 <- tmp$vectors %*% diag(exp(tmp$values)) %*% solve(tmp$vectors)
return(P1)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.