Thanks for your advices. The fact is that even when the input of the eigenvalues is optional in g_anaeig, they must be provided besides the eigenvectors in order to obtain the same result as with the perl script. Williams