library(Biodem)
`%^%` <- function(mat, n) Biodem::mtx.exp(mat, n)