ape (version 5.1)

matexpo: Matrix Exponential

Description

This function computes the exponential of a square matrix using a spectral decomposition.

Usage

matexpo(x)

Arguments

x

a square matrix of mode numeric.

Value

a numeric matrix of the same dimensions than `x'.

Examples

# NOT RUN {
### 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))
# }