| 查看: 190 | 回复: 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 |
回复此楼» 猜你喜欢
国社科又开始会评了，不知道这次命运如何
已经有13人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
30