[Fortran] 纯文本查看 复制代码
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