在fortran6.6中有imsl数据库，今天看帖子看到：
见http://www.programfan.com/club/showtxt.asp?id=303821

求矩阵特征值和特征向量以及为它们自动排序的程序。想问下：这个EVCRG在fortran6.6中也可以用吗？我用eig(a,w=vector)求的特征值和特征向量，发现特征值不是按照顺序输出的，是按照特征值的绝对值从大到小输出的，如果这个函数可以用那多好啊，省去好多功夫啊


CALL EVCRG (A, EVAL, EVEC [,…])

A — Floating-point array containing the matrix. (Input) 
那个实矩阵

EVAL — Complex array of size N containing the eigenvalues of A in decreasing order of
magnitude. (Output) 按降序排列的特征值数组

EVEC — Complex array containing the matrix of eigenvectors. (Output)
特征向量组成的矩阵

The J-th eigenvector, corresponding to EVAL(J), is stored in the J-th column. 
EVEC的第J列，对应第J个特征值EVAL(J),

Each vector is normalized to have Euclidean length equal to the value one.
每个特征向量都是单位长度