| 查看: 1031 | 回复: 3 | |||
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 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 配置一台1万左右的计算机跑Fortran程序以及高斯的计算
已经有8人回复
- 求高人帮忙看一下我这个Fortran程序的错误出在哪里？
已经有12人回复
- 物理模拟是用c语言好还是fortran好？
已经有30人回复
- 我编的Simpson积分法fortran程序给不出结果，大侠们看看哪里出了问题？
已经有4人回复
- Fortran 95／2003科学计算与工程 对应的各章程序源代码
已经有16人回复
- 有关fortran的一次而问题，希望大家能帮帮忙，谢谢
已经有4人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】用fortran怎么描述下面这样的情况，我是新手，请大家帮忙
已经有8人回复
- 【求助】初学fortran，
已经有15人回复
- 【求助】请问一下fortran的可视化编程
已经有7人回复
- 【原创】elk中的状态方程拟合工具eos
已经有128人回复
- 【求助】fortran计算时出现NaN
已经有11人回复
Jasonluo
至尊木虫 (知名作家)
中央书籍处书记
- 应助: 61 (初中生)
- 贵宾: 0.01
- 金币: 22673.1
- 散金: 166
- 红花: 14
- 帖子: 9267
- 在线: 681小时
- 虫号: 518969
- 注册: 2008-03-05
- 性别: GG
- 专业: 自然地理学
2楼2012-05-16 20:59:37
jwwang_lq
铁虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 17
- 帖子: 7
- 在线: 1.6小时
- 虫号: 1719996
- 注册: 2012-03-27
- 性别: GG
- 专业: 生物大分子结构与功能
3楼2012-05-17 08:51:20
daizhijuan
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 12.5
- 帖子: 1
- 在线: 27分钟
- 虫号: 2184351
- 注册: 2012-12-13
- 专业: 凝聚态物性 II ：电子结构
4楼2012-12-13 16:52:02
20