应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 刚开始学Fortran 帮忙看一下这个程序 为什么算不出结果
41/1返回列表
查看: 836  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

caoweiyue

新虫 (小有名气)

[求助] 刚开始学Fortran 帮忙看一下这个程序 为什么算不出结果

下面的程序算的是全选主元的高斯消去法， 我是照徐士良的算法集编的，但是结果都是0
program main
implicit none
integer::i,L
integer,parameter::n=4
real::A(n,n),b(n),x(n)
data A / 0.2368,0.1968,0.1582,1.1161,0.2471,0.2071,1.1675,0.1254,&
         &0.2568,1.2168,0.1768,0.1397,1.2671,0.2271,0.1871,0.1490 /
data b / 1.8471,1.7471,1.6471,1.5471 /
call agaus(A,b,n,L)
if(L==0) then
  write(*,*) "fail"
else
  do i=1,n
    write(*,"('x',i1,'=',f6.4)" i,x(i)
  end do
end if
stop
end program

subroutine agaus(A,b,n,L)
implicit none
integer:: i,j,k,n,is,js,L
real:: D, T
real:: A(n,n),b(n),x(n)

L=1
do k=1,n-1
  D=0.0
!选出了最大的，并保存其位置
   do i=k,n
     do j=k,n
           if( abs(A(i,j))>D ) then
           D=A(i,j)
                   is=i
                    js=j
           end if  
     end do
   end do
  if(D+1.0==1.0) then
    L=0
    return
  end if
!列交换
  if(js/=k) then
    do i=1,n
          T=A(i,k)
          A(i,k)=A(i,js)
      A(i,js)=T
    end do
  end if
!行交换
  if(is/=k) then
    do j=k,n
          T=A(k,j)
          A(k,j)=A(is,j)
          A(is,j)=A(k,j)
        end do
        T=b(k)
        b(k)=b(is)
        b(is)=b(k)
  end if
!化成三角矩阵
  do i=k+1,n
    do j=k,n
          A(i,j)=A(i,j)-A(i,k)/A(k,k)*A(k,j)
        end do
        b(i)=b(i)-A(i,k)/A(k,k)*b(k)
  end do
end do
if(A(n,n)+1.0==1.0) then
  L=0
  return
end if
!回代
x(n)=b(n)/A(n,n)
  do i=n-1,1,-1
    T=0.0
          do j=i+1,n
            T=T+A(i,j)*x(j)
          end do
        x(i)=(b(i)-T)/A(i,i)
  end do
  return
end subroutine agaus
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-16 18:39:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasonluo

至尊木虫 (知名作家)

中央书籍处书记

【答案】应助回帖

感谢参与，应助指数 +1
楼主中间那个害羞的表情是啥意思啊？呵呵
赞一下 回复此楼
朝为田舍郎，暮登天子堂。将相本无种，学子当自强！
2楼2012-05-16 20:59:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jwwang_lq

铁虫 (初入文坛)

【答案】应助回帖

感谢参与，应助指数 +1
write语句里应该输出b(i)，而不是x(i)
赞一下 回复此楼
老实人
3楼2012-05-17 08:51:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daizhijuan

新虫 (初入文坛)

【答案】应助回帖

write(*,"('x',i1,'=',f6.4)" i,x(i)
少了右边的括号
real:: A(n,n),b(n),x(n)
子程序中的这句之前最好也给n一个值：比如 parameter（n=4）
赞一下 回复此楼
4楼2012-12-13 16:52:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 caoweiyue 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 芦lty 2026-03-25 4/200 2026-03-25 23:25 by 芦lty
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 一志愿哈工大，085400，320，求调剂 +4 gdlf9999 2026-03-24 4/200 2026-03-25 23:01 by boxking200
[考研] 0854电子信息求调剂 +7 α____ 2026-03-22 9/450 2026-03-25 13:37 by α____
[考研] 0703化学调剂，求导师收 +7 天天好运来上岸� 2026-03-24 7/350 2026-03-24 20:26 by peike
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 292求调剂 +4 鹅鹅鹅额额额额� 2026-03-24 4/200 2026-03-24 16:41 by peike
[考研] 291求调剂 +3 HanBeiNingZC 2026-03-24 3/150 2026-03-24 16:34 by barlinike
[考研] 279分求调剂 一志愿211 +18 chaojifeixia 2026-03-19 20/1000 2026-03-24 10:34 by dolphin_ycj
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 一志愿陕师大生物学071000，298分，求调剂 +3 SYA！ 2026-03-23 3/150 2026-03-23 19:09 by macy2011
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪！？！
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
信息提示
关闭
请填处理意见
关闭