Andreas: I actually didn't know about eigfact. I avoided eig since I knew most of the eigenvalues were zero (and positive) so I thought it was a waste to use eig to find all the eigenvalues, then discard the near zeros ones.
I'll look into eigfact. Thanks for all your help!
