Нужно представить Вашу матрицу в виде A = Р^(-1) * D * P, где D - диагональная. Потом использовать тот факт, что exp(A) = Р^(-1) * exp(D) * P.