24小时热门版块排行榜

返回列表

jihe01

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 2867.8
帖子: 336
在线: 291.5小时
虫号: 361438
注册: 2007-05-04
性别: GG
专业: 土壤生态学

[交流] 【求助】程序是产生随机数,请看以下问题，谢谢！

程序是产生随机数的,f90，我用fortran 6.6运行出现 Access Violation 好像是IDUM=1的问题

FUNCTION RAN1(IDUM)

   DIMENSION R(97)
   PARAMETER(M1=259200,IA1=7141,IC1=54773,RM1=1./M1)
   PARAMETER(M2=134456,IA2=8121,IC2=28411,RM2=1./M2)
   PARAMETER(M3=243000,IA3=4561,IC3=51349)

   DATA IFF/0/

   IF(IDUM.LT.0.OR.IFF.EQ.0) THEN
      IFF=1
      IX1=MOD(IC1-IDUM,M1)
      IX1=MOD(IA1*IX1+IC1,M1)
      IX2=MOD(IX1,M2)
      IX1=MOD(IA1*IX1+IC1,M1)
      IX3=MOD(IX1,M3)

      DO 11 J=1,97
      IX1=MOD(IA1*IX1+IC1,M1)
      IX2=MOD(IA2*IX2+IC2,M2)
      R(J)=(real(IX1)+real(IX2)*RM2)*RM1
11 CONTINUE
      IDUM=1.0
   ENDIF

   IX1=MOD(IA1*IX1+IC1,M1)
   IX2=MOD(IA2*IX2+IC2,M2)
   IX3=MOD(IA3*IX3+IC3,M3)
   J=1+(97*IX3)/M3
   IF(J.GT.97.OR.J.LT.1)PAUSE
   RAN1=R(J)
   R(J)=(real(IX1)+real(IX2)*RM2)*RM1

   RETURN
   END function

回复此楼

» 猜你喜欢

各位大神，目前国内有哪些比较好用的逆合成软件？已经有14人回复
国社科系统bug了，是不是要放榜了？已经有10人回复
现代”学阀”该如何界定已经有9人回复
课题组招2027级博士上海工程技术大学激光智能制造方向已经有4人回复
上海工程技术大学激光智能制造课题组 2027级博士研究生招生已经有3人回复
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告已经有11人回复
申博发邮件已经有11人回复
我的奶奶已经有3人回复

1楼 2011-04-03 17:43:18

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

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

返回列表