| 查看: 395 | 回复: 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 超级感谢!!!!!!! |













回复此楼