SUBROUTINE VIBSOL(NQ1max,NQ1,RQ1,TMS,eigenvaluor,Vpot,eigenvalue)
    use imsl
    IMPLICIT REAL*8 (A-H,O-Z)
    DIMENSION  RQ1(NQ1),HAR(NQ1,NQ1),Vpot(NQ1),eigenvalue(NQ1),
    $WK1(1),WK2(1),eigenvaluor(NQ1,NQ1)

    RL=RQ1(NQ1)-RQ1(1)+2.D0*(RQ1(2)-RQ1(1))
    NDVR=NQ1+1
    DO  I=1,NQ1
       HAR(I,I)=Vpot(I)+DVRKE1(I,NDVR,RL,TMS)
    DO  J=1,I-1
       HAR(J,I)=DVRKE2(I,J,NDVR,RL,TMS)
    HAR(I,J)=HAR(J,I)
    enddo
    enddo
    eigenvalue=eig(HAR,v=eigenvaluor)
    RETURN
    END
    
运行时提示堆栈出错，程序崩溃。调试时发现崩在了 EIG函数上，提示如截图。第一次用IMSL库，懂的不多。请教高手，代码哪里有问题？ 谢谢！