[数学库] 利用MKL求解矩阵方程出错

program main
  use mkl95_precision
  use mkl95_lapack
  implicit none
  integer i,info
  integer ipiv(3)
  real*8 ::a(3,3),b(3,1),X(3,1)
  a=0
  b=1
  
  forall( i=1:3)
  a(i,i)=1
  end forall
  print *,'Matrix A is:'
  do i=1,3
      print *,a(i,:)
  end do
  print *,'Matrix B is:'
  print *,b
  
  call getrf(a,ipiv,info)
  call getrs(a,ipiv,b,'N',info)
  print *,'Solution X is:'
  print *,b
end program