24小时热门版块排行榜    

查看: 389  |  回复: 2

ilovexiaomu

金虫 (小有名气)

[求助] 麻烦看哈这个很简单的求线性方程的问题?已有1人参与

各位,遇到一个程序,在网上直接下载的,编译没错,但是运行时报错,不知什么原因,还望各位解答
!**********************************************
! THIS EXAMPLE USES THE LAPACK ROUTINE DGESV
! TO SOLVE A SYSTEM OF LINEAR EQUATIONS AX=B
! A = [1, 2, 3; 4, 5, 6; 7, 8 10]            
! B = [1, 0; 0, 1; 0, 0]
! *********************************************
program lapack_test
  integer ipiv(3), info, i, j
  double precision A(3,3), B(3,2)
  A(1,1)=1
  A(1,2)=2
  A(1,3)=3
  A(2,1)=4
  A(2,2)=5
  A(2,3)=6
  A(3,1)=7
  A(3,2)=8
  A(3,3)=10
  
  B(1,1)=1
  B(2,1)=0
  B(3,1)=0
  B(1,2)=0
  B(2,2)=1
  B(3,2)=0

  call dgesv (3,2,A,3,ipiv,B,3,info)
  if(info .EQ. 0) then
     do i=1,3
        write(*,'(2F8.3)') (B(i,j), j=1,2)
     enddo
  endif
end program lapack_test
不知道是自己的安装软件有问题还是那里出了什么问题,还望各位解答
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ilovexiaomu

金虫 (小有名气)

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖


感谢参与,应助指数 +1
ilovexiaomu: 金币+1, 有帮助, 感谢 2015-05-21 16:05:55
相关版块跳转 我要订阅楼主 ilovexiaomu 的主题更新
信息提示
请填处理意见