| 查看: 671 | 回复: 6 | |||
[求助]
求助关于program的用法 已有1人参与
|
| 我在电脑上运行一个高斯赛德尔迭代的计算代码，迭代是用子程序给出的。在主程序中，我在首位加入了program main运行后就会出现This name does not match the unit name.不加program就能运行。哪位大神知道只是为什么啊？ |
回复此楼» 猜你喜欢
调剂
已经有6人回复
-
280求调剂
已经有6人回复
-
302求调剂一志愿华中师范大学
已经有5人回复
-
346求调剂
已经有4人回复
-
材料专硕322分
已经有9人回复
-
372分材料与化工（085600）一志愿湖南大学求调剂
已经有5人回复
-
318求调剂，计算材料方向
已经有10人回复
-
310求调剂
已经有14人回复
-
315求调剂
已经有3人回复
-
265求调剂
已经有13人回复
|
代码如下 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 (初中生)
- 金币: 26160.7
- 散金: 50
- 红花: 10
- 帖子: 5672
- 在线: 830.9小时
- 虫号: 1331748
- 注册: 2011-06-26
- 专业: 金属结构材料
3楼2015-11-11 11:49:32
4楼2015-11-11 11:50:29
78506174
至尊木虫 (知名作家)
- 应助: 70 (初中生)
- 金币: 26160.7
- 散金: 50
- 红花: 10
- 帖子: 5672
- 在线: 830.9小时
- 虫号: 1331748
- 注册: 2011-06-26
- 专业: 金属结构材料
5楼2015-11-11 11:52:15
6楼2015-11-11 11:52:48
78506174
至尊木虫 (知名作家)
- 应助: 70 (初中生)
- 金币: 26160.7
- 散金: 50
- 红花: 10
- 帖子: 5672
- 在线: 830.9小时
- 虫号: 1331748
- 注册: 2011-06-26
- 专业: 金属结构材料
7楼2015-11-11 11:55:47
