square matrix exponential
expm(X)
:X square matrix with real or complex entries. :
X is a square matrix expm(X) is the matrix
expm(X) = I + X + X^2 /2 + ...
The computation is performed by first block-diagonalizing X and then applying a Pade approximation on each block.
X=[1 2;3 4]
expm(X)
`logm`_(expm(X))