大家好：

      fortran中有什么可以求复矩阵特征向量的？我目前用的是eig内部函数，先用实矩阵的特征向量验证，但求出来与matlab的不相符，但经验证matlab的是对的，特征值一样，但特征向量不同，而且不是简单符号的问题，请问问题出在哪里呢？以下是程序：

use imsl

  implicit none

real:: a(3,3)=(/2,0,0,1,2,-1,1,0,1/)
real::eigenvalue(3)
real::eigenvector(3,3)
integer i

eigenvalue=eig(a, v=eigenvector)

do i=1,3
write(*,*)  eigenvalue(i)
write(*,*)  eigenvector(:,i)

write(*,*)

enddo
end

还有什么方法是求复数矩阵的特征值与特征向量的？