| 查看: 199 | 回复: 0 | |||
| 当前主题已经存档。 | |||
[交流]
随机数发生器（Fortran77编写的，测试可用）
|
|||
|
real function ranw(idum) Parameter (Mbig=2**30-2, Xinvers=1./Mbig) data ibit/ 1/ Integer IX(55) save if (ibit.ne.0) then ibit=0 cc cc fill up the vector ix with some random integers, which are cc not all even cc if (idum.eq.0) pause 'use nonzero value of idum' idum=abs(mod(idum,Mbig)) ibit=0 Ix(1)=871871 Do i=2,55 Ix(i)=mod(Ix(i-1)+idum,Ix(i-1)) Ix(i)=max(mod(Ix(i),Mbig),idum) end do j=24 k=55 cc cc warm up the generator cc do i=1,1258 Ix(k)=mod(Ix(k)+Ix(j),Mbig) j=j-1 if (j.eq.0) j=55 k=k-1 if (k.eq.0) k=55 end do end if cc cc this is where execution usually starts: cc Ix(k)=mod(Ix(k)+Ix(j),Mbig) j=j-1 if (j.eq.0) j=55 k=k-1 if (k.eq.0) k=55 ranw=Ix(k)*Xinvers end |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有13人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
申请科研项目纯靠实力能成功吗？
已经有11人回复
-
极限拉扯，如何抉择
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
200