| 查看: 386 | 回复: 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 超级感谢！！！！！！！ |
回复此楼» 猜你喜欢
276求调剂。有半年电池和半年高分子实习经历
已经有4人回复
-
284求调剂
已经有6人回复
-
306求调剂
已经有7人回复
-
291求调剂
已经有5人回复
-
070300化学求调剂
已经有4人回复
-
070300，一志愿北航320求调剂
已经有5人回复
-
085600材料与化工
已经有9人回复
-
0703 调剂
已经有6人回复
-
招08考数学
已经有11人回复
-
276求调剂
已经有4人回复
