briefly we need to find all Eigenvectors:

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