| 查看: 610 | 回复: 6 | ||
[求助]
求助关于program的用法 已有1人参与
|
| 我在电脑上运行一个高斯赛德尔迭代的计算代码，迭代是用子程序给出的。在主程序中，我在首位加入了program main运行后就会出现This name does not match the unit name.不加program就能运行。哪位大神知道只是为什么啊？ |
回复此楼» 猜你喜欢
基金申报
已经有5人回复
-
基金委咋了？2026年的指南还没有出来？
已经有7人回复
-
国自然申请面上模板最新2026版出了吗？
已经有17人回复
-
纳米粒子粒径的测量
已经有8人回复
-
疑惑？
已经有5人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有5人回复
-
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
推荐一本书
已经有12人回复
-
常年博士招收(双一流，工科)
已经有4人回复
|
代码如下 IMPLICIT NONE DIMENSION A(4,4),B(4),X(4) DOUBLE PRECISION A,B,X REAL:: EPS INTEGER::I,L DATA A/7.0,9.0,-2.0,1.0,2.0,15.0,-2.0,3.0,1.0,3.0,11.0,2.0,-2.0,-2.0,5.0,13.0/ DATA B/4.0,7.0,-1.0,0.0/ EPS=1.0E-06 CALL AGSDL(A,B,4,X,EPS,L) IF (L.NE.0) THEN WRITE(*,10) (X(I),I=1,4) END IF 10 FORMAT(1X,D15.6) END SUBROUTINE AGSDL(A,B,N,X,EPS,L) DIMENSION A(N,N),B(N),X(N) DOUBLE PRECISION A,B,X,T,S,P,Q INTEGER::N DO 5 I=1,N IF (ABS(A(I,I))+1.0.EQ.1.0) THEN L=0 WRITE(*,100) RETURN END IF 5 CONTINUE 100 FORMAT(1X,'FAIL') L=100 DO 10 I=1,N 10 X(I)=0.0 20 P=0.0 L=L-1 DO 50 I=1,N T=X(I) S=0.0 DO 30 J=1,N IF(J.NE.I) S=S+A(I,J)*X(J) 30 CONTINUE X(I)=(B(I)-S)/A(I,I) Q=ABS(X(I)-T)/(1+ABS(X(I))) IF (Q.GT.P) P=Q 50 CONTINUE IF((P.GE.EPS).AND.(L.NE.0)) GOTO 20 IF(L.EQ.0) WRITE(*,100) RETURN END SUBROUTINE AGSDL |
2楼2015-11-10 21:19:59
78506174
至尊木虫 (知名作家)
- 应助: 70 (初中生)
- 金币: 25884.2
- 散金: 50
- 红花: 10
- 帖子: 5616
- 在线: 829.7小时
- 虫号: 1331748
- 注册: 2011-06-26
- 专业: 金属结构材料
3楼2015-11-11 11:49:32
4楼2015-11-11 11:50:29
78506174
至尊木虫 (知名作家)
- 应助: 70 (初中生)
- 金币: 25884.2
- 散金: 50
- 红花: 10
- 帖子: 5616
- 在线: 829.7小时
- 虫号: 1331748
- 注册: 2011-06-26
- 专业: 金属结构材料
5楼2015-11-11 11:52:15
6楼2015-11-11 11:52:48
78506174
至尊木虫 (知名作家)
- 应助: 70 (初中生)
- 金币: 25884.2
- 散金: 50
- 红花: 10
- 帖子: 5616
- 在线: 829.7小时
- 虫号: 1331748
- 注册: 2011-06-26
- 专业: 金属结构材料
7楼2015-11-11 11:55:47