I have been trying to compute the eig decomposition of 13000 x 13000 sparse matrix D using matlab command eigs (single (D)), but I get 'out of error' although memory wasn't fully utilized (with a 4GB RAM) , the help I saw online suggested using numerical methods like ARPACK (Fotran77 software)...