You can use GENEIG Call  



CALL GENEIG( eval, evecs, sym-matrix1, sym-matrix2 ) ;
The GENEIG subroutine computes eigenvalues and eigenvectors of a generalized eigenproblem.

果真如此，谢谢老兄了。