24小时热门版块排行榜    

CyRhmU.jpeg
查看: 359  |  回复: 0

dianananana

新虫 (初入文坛)

[交流] 随机数生成

请教大神,这个生成随机数的代码里的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
超级感谢!!!!!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dianananana 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见