Added the eig method for eigen-system calculation