Re: eigen values

[george brida]

Thank you very much Nir.

On Thu, Mar 24, 2016 at 3:20 PM, Nir Krakauer <address@hidden> wrote:

[eigenvecs, eigenvals] = eig (V);

eigenvals_keep = diag(eigenvals) != 0; #or abs(diag(eigenvals)) > tol to also neglect eigenvalues that are close to zero

R = eigenvals(eigenvals_keep, eigenvals_keep);

V = eigenvecs(:, eigenvals_keep);