expmat function

Matrix exponential