Re: eigen values

[Nir Krakauer]

[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);