Hi,

I am using dgeevx routine in LAPACK to solve a system of ordinary differential equations (chemical kinetics reaction model). This routine returns the eigenvectors that are normalized to have Euclidean norm equal to 1 and largest component real. How do I fix the sign of each eigenvector so that I can obtain the time course of each chemical species? I know that the concentrations of all species should be positive at all times. Any suggestions or pointers are welcome.

Thanks,

mmgk