I have a 4n system of differential equations. I got all the eigen values of the equations in complex numbers. For real eigen values I can solve the problem in Matlab using Matrix exponentiation method . Can anyone suggest some idea for solving the same in Matlab when eigen values are complex.