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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 刚开始学Fortran 帮忙看一下这个程序 为什么算不出结果
41/1返回列表
查看: 825  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +3 哈哈157349 2026-03-21 3/150 2026-03-24 14:11 by zhyzzh
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 291 求调剂 +4 化工2026届毕业� 2026-03-21 5/250 2026-03-23 16:46 by 化工2026届毕业�
[考研] 350求调剂 +6 weudhdk 2026-03-19 6/300 2026-03-23 15:47 by tangyuan0840221
[考研] 求调剂材料学硕080500，总分289分 5+3 @taotao 2026-03-19 21/1050 2026-03-23 10:17 by 冠c哥
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 070300，一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 一志愿华中农业071010，总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
信息提示
关闭
请填处理意见
关闭