| 查看: 597 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助关于program的用法已有1人参与
|
||
| 我在电脑上运行一个高斯赛德尔迭代的计算代码，迭代是用子程序给出的。在主程序中，我在首位加入了program main运行后就会出现This name does not match the unit name.不加program就能运行。哪位大神知道只是为什么啊？ |
回复此楼» 猜你喜欢
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有7人回复
-
申请2026年博士
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
|
代码如下 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
4楼2015-11-11 11:50:29
6楼2015-11-11 11:52:48