solving expm

expm(X) = V*diag(exp(diag(D)))/V

scaling expm using eigenvalue https://math.stackexchange.com/questions/3480626/eigenvalues-eigenvectors-of-scaled-matrix