[数学库] 调用MKL的dgesv函数时出现错误

program main
use lapack95
    implicit none

integer,parameter::size=5
integer::i,j
integer::n,lda,ldb
integer::nrhs
integer::info
integer::ipiv(size)
real::a(size,size),b(size,size)
n=size;lda=size;ldb=size
nrhs=size


forall(i=1:size,j=1:size,i>j)     a(i,j)=2
forall(i=1:size,j=1:size,i==j)    a(i,j)=3
forall(i=1:size,j=1:size,i<j)     a(i,j)=4

forall(i=1:size,j=1:size,i>j)     b(i,j)=0
forall(i=1:size,j=1:size,i==j)    b(i,j)=1
forall(i=1:size,j=1:size,i<j)     b(i,j)=0

call dgesv(n,nrhs,a,lda,ipiv,b,ldb,info)
write(*,*)'solution:'
write(*,*)b
write(*,*)'info=',info

end program