| 查看: 416 | 回复: 0 | |||
[交流]
随机数生成
|
|
请教大神，这个生成随机数的代码里的J是干嘛的？ 我运行时J总是小于零（等于-2405）导致程序暂停，按回车键后报错 C RANDOM NUMBER GENERATORS: C C RAN2(IDUM) RETURNS A UNIFORM RANDOM NUMBER ON (0,1). C THE ROUTINE IS FROM: PRESS, W.H., B.P. FLANNERY, S.A. C TEUKOLSKY, AND W.T. VETTERLING. 1987. NUMERICAL RECIPES. C CAMBRIDGE UNIVERSITY PRESS, CAMBRIDGE. C FUNCTION RAN2(IDUM) DIMENSION IR(97) DATA M/714025/ DATA IA/1366/ DATA IC/150889/ DATA IFF/0/ RM=1.0/FLOAT(M) IF (IDUM.LT.0.OR.IFF.EQ.0) THEN IFF=1 IDUM=MOD(IC-IDUM,M) DO 11 J=1,97 IDUM=MOD(IA*IDUM+IC,M) IR(J)=IDUM 11 CONTINUE IDUM=MOD(IA*IDUM+IC,M) IY=IDUM ENDIF J=1+(37*IY)/M IF (J.GT.97.OR.J.LT.1) PAUSE IY=IR(J) RAN2=IY*RM IDUM=MOD(IA*IDUM+IC,M) IR(J)=IDUM RETURN END 超级感谢！！！！！！！ |
回复此楼» 猜你喜欢
关于如何从代码看上不上会
已经有21人回复
-
chemdraw
已经有4人回复
-
大龄残疾硕士的一点执念
已经有24人回复
-
祈祷青基必中
已经有10人回复
-
b口会评
已经有6人回复
-
信息学部会评时间
已经有3人回复
-
在职考研
已经有7人回复
-
MSER送审了还被拒稿
已经有5人回复
-
硝基还原
已经有5人回复
-
化学口会评？应用基础研究和基础研究，希望能分开
已经有3人回复