| 查看: 984 | 回复: 11 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
icecream8508金虫 (著名写手)
|
[交流]
【求助】新人请教fortran编程问题！[求助完成]
|
||
|
刚接触Fortran不久。因为自己求解的东西需要迭代，而自己对迭代的认识就是附初值，然后更新每次的值，直到收敛。大致思路也该如此，但是在程序实现的时候，总是不对的。 因此，今天就想从最简单的迭代开始算起。想先把最基本的一个迭代的例子走通，但是结果仍然不对。 也问了几个人，但是大家对于程序的编写都有自己习惯的语言，觉得问完之后，更茫然了。 在检查程序的时候，不知道错在哪里了？ 不知，大家在审查程序时有什么好方法吗？ [ Last edited by wangmonk on 2009-6-7 at 20:44 ] |
回复此楼» 猜你喜欢
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
40+中年人失业后，还能找到好工作吗？
已经有5人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
icecream8508
金虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 939.6
- 散金: 72
- 红花: 18
- 帖子: 1537
- 在线: 557.8小时
- 虫号: 630888
- 注册: 2008-10-19
- 专业: 污染控制化学
7楼2009-06-04 09:25:01
icecream8508
金虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 939.6
- 散金: 72
- 红花: 18
- 帖子: 1537
- 在线: 557.8小时
- 虫号: 630888
- 注册: 2008-10-19
- 专业: 污染控制化学
|
我今天做的简单例子的编程附着如下： !******* 练习迭代算法******* !******** 所需要求解的方程****** ! x(1)=0.2*x(2)+0.1*x(3)+0.3 ! x(2)=0.2*x(1)+0.1*x(3)+1.5 ! x(3)=0.2*x(1)+0.1*x(2)+2 !******* 1.jACOBI 迭代******* program main real B(3) real A(3,3),x(3),s(30) integer N,k integer m real(kind=8) :: eps real(kind=8) :: xo integer,parameter:: NM=1000 eps=1.0e-6 !***********系数矩阵************* B(1)=0.3 B(2)=1.5 B(3)=2 A(1,2)=0.2 A(1,3)=0.1 A(2,1)=0.2 A(2,3)=0.1 A(3,1)=0.2 A(3,2)=0.1 !***************迭代公式************ ! do k=1,NM,1 !Do i=1,3 ! x(i)=0.0 !end do 20 Do k=1,NM,1 Do i=1,3 s(i)=0.0 Xo=x(i) x(i)=0.0 Do j=1,3 If(i.NE.j) then ! S(i)=S(i)+A(i,j)*Xo S(i)=S(i)+A(i,j)*X(j) End if End do x(i)=B(i)-s(i) If(abs(xo-x(i)).gt.eps) then Go to 20 else exit end if end do end do !write(*,*) (x(i),i=1,3) do j=1,3 write(*,*) x(j) end do end program main |
2楼2009-06-03 11:50:00
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1432.8
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 442.1小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
3楼2009-06-03 21:38:35
maomao1210
金虫 (正式写手)
- 程序强帖: 5
- 应助: 2 (幼儿园)
- 金币: 1432.8
- 散金: 242
- 红花: 16
- 沙发: 1
- 帖子: 991
- 在线: 442.1小时
- 虫号: 253215
- 注册: 2006-05-20
- 性别: MM
- 专业: 考古理论
4楼2009-06-03 21:39:31
200
