Fortran求n*n阶矩阵束 (A,B)的最小特征对 (θ, u) ，给出程序

以下是Fortran语言的程序，用于求解n*n阶矩阵对(A,B)的最小特征对(θ, u)。

PROGRAM MIN_EIGENPAIR
    IMPLICIT NONE
    INTEGER :: N, LWORK, INFO
    REAL :: THETA
    REAL, DIMENSION(:,:), ALLOCATABLE :: A, B, U, WORK
    ! 获取矩阵维数
    WRITE(*,*) '请输入矩阵维数N：'
    READ(*,*) N
    ! 动态分配矩阵空间
    ALLOCATE(A(N,N), B(N,N), U(N,N), WORK(2*N), STAT=INFO)
    IF (INFO /= 0) THEN
        WRITE(*,*) '动态分配矩阵空间失败！'
        STOP
    END IF
    ! 获取矩阵A和B的值
    WRITE(*,*) '请输入矩阵A的值：'
    READ(*,*) A
    WRITE(*,*) '请输入矩阵B的值：'
    READ(*,*) B
    ! 调用LAPACK库中的dsygv函数求解最小特征对
    LWORK = 2*N
    CALL DSYGV(1, 'V', 'L', N, A, N, B, N, THETA, WORK, LWORK, INFO)
    IF (INFO /= 0) THEN
        WRITE(*,*) '求解最小特征对失败！'
        STOP
    END IF
    ! 输出结果
    WRITE(*,*) '最小特征值为：', THETA
    WRITE(*,*) '最小特征向量为：'
    U = A
    WRITE(*,*) U
END PROGRAM MIN_EIGENPAIR

需要注意的是，该程序依赖于LAPACK库中的dsygv函数，因此需要在编译时链接该库。另外，程序中使用了动态分配矩阵空间的方式，需要在程序结束时释放空间，以避免内存泄漏。