powered by
This function computes the exponential of a square matrix using a spectral decomposition.
matexpo(x)
a numeric matrix of the same dimensions than `x'.
a square matrix of mode numeric.
Emmanuel Paradis
### a simple rate matrix: m <- matrix(0.1, 4, 4) diag(m) <- -0.3 ### towards equilibrium: for (t in c(1, 5, 10, 50)) print(matexpo(m*t))
Run the code above in your browser using DataLab