Eigen decomp

We need a method to find eigenvalues/eigenvectors that works for non-Hermitian matrices.