24小时热门版块排行榜

返回列表

dianananana

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9
帖子: 9
在线: 8.9小时
虫号: 4119036
注册: 2015-10-04
专业: 遥感机理与方法

[交流] 随机数生成

请教大神，这个生成随机数的代码里的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
超级感谢！！！！！！！

回复此楼

» 猜你喜欢

1楼 2018-03-22 11:39:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 dianananana 的主题更新

返回列表