briefly we need to find all Eigenvectors:
- find all Eigenvalues of the coefficient matrix by solving
- for each find where is identify matrix and is correspond Eigenvector
- if the is repeated eigenvalue (equivalent to repeated root), need to use to find the generalized eigenvector, this process also called
Jordan Chain
- lastly construct the homogeneous solution