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

回复此楼

» 猜你喜欢

假如你的研究生提出不合理要求已经有6人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
论文终于录用啦！满足毕业条件了已经有28人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复
实验室接单子已经有3人回复
磺酰氟产物，毕不了业了！已经有8人回复

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

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

智能机器人

Robot (super robot)

我们都爱小木虫

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

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

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

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

返回列表