| 查看: 634 | 回复: 4 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 lvhuidan 的 5 个金币 | |||
[交流]
【求助】求助：关于compag visual fortran6 的运行问题
|
|||
| 我不懂编程，只是最近上《化工计算》这门课需要用到书上给出的fortran程序。我把书上的程序编好后，编译和链接都没有错误，可是一运行，就没有结果出现，试了2段程序，都是包括一个主程序和一个子程序的，用的是compag visual fortran6。不知道怎么回事，请高手指点！ |
回复此楼» 猜你喜欢
实验室接单子
已经有3人回复
-
要不要辞职读博？
已经有5人回复
-
不自信的我
已经有10人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555220
- 散金: 27828
- 红花: 2147
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
2楼2009-11-21 16:34:07
jjdg(金币+0,VIP+0):欢迎讨论！请继续支持程序软件版！ 11-23 17:11
|
这是试过的一段程序： DIMENSION C(4,5), X(4) REAL *8 C READ(*,*) M WRITE(*,10) M N=M+1 DO 2 I=1,M 2 READ(*,*) (C(I,J), J=1,N) WRITE(*,11) DO 3 I=1,M 3 WRITE(*,13) (C(I,J),J=1,N) CALL SIMUL(M,N,C,X) WRITE(*,12) DO 1 J=1,M 1 WRITE(*,14) X(J) 10 FORMAT(1X,'1','M=',I5) 11 FORMAT(//,1X,'COEFFICIENT MATRIX=',/) 12 FORMAT(//,1X,'SOLUTION VECTOR=',/) 13 FORMAT(8F10.4) 14 FORMAT(F15.4) STOP END SUBROUTINE SIMUL(M,N,C,X) IMPLICIT REAL*8 (A-H,O-Z) REAL*4X DIMENSION C(4,5),X(4) DO 1 K=1,M PIVOT=C(K,K) IF(DABS(PIVOT).LE.0.001) GOTO 2 DO 3 J=K,N 3 C(K,J)=C(K,J)/PIVOT DO 4 I=1,M IF(I.EQ.K) GOTO 4 PIVI=C(I,K) DO 5 J=K,N 5 C(I,J)=C(I,J)-PIVI*C(K,J) 4 CONTINUE 1 CONTINUE DO 6 J=1,M 6 X(J)=C(J,N) RETURN 2 WRITE(*,100) K 100 FORMAT(1X,'PIVOT ELEMENT',I5,1X,'IS CLOSE TO 0') RETURN END 数据 1M= 4 .47 0 -.012 -1 0 0 .92 0 -1.7857 5 .49 0 -.73 0 5 0 .08 -.01 2.5 990 从“数据..........990” 这一段在程序里不知怎么输入，可能就是这里的问题，因为程序里没有给出M的值，不知我想的对不对？请楼上的指点 |
3楼2009-11-22 10:40:23
4楼2009-11-24 16:31:52
senlia
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 贵宾: 3.521
- 金币: 2178.1
- 红花: 1
- 帖子: 223
- 在线: 6.7小时
- 虫号: 849005
- 注册: 2009-09-16
- 性别: GG
- 专业: 计算数学与科学工程计算
★ ★ ★
nono2009(金币+3,VIP+0):无语。 11-24 19:40
nono2009(金币+0,VIP+0):Perfect. 11-24 20:51
nono2009(金币+3,VIP+0):无语。 11-24 19:40
nono2009(金币+0,VIP+0):Perfect. 11-24 20:51
|
首先有个错误 REAL*4X X前面应该有个空格 READ(*,*) (C(I,J), J=1,N)是需要从键盘输入数据 你可以把 .47到990这一段复制下来 在程序运行界面的最上方蓝色块右击 选择 编辑--粘贴 或者从文件读入 open(41,file='data.txt') READ(41,*) (C(I,J), J=1,N) 另外 在需要的地方输入提示语 不然在程序运行过程中你需要做什么自己都不知道 另外 在你粘贴上来的程序中没有program这一行 不知道你是不知道还是没有粘贴过来 附件里面是程序、数据文件和运行结果 不知道是不是你想要的结果 |
5楼2009-11-24 18:13:10