| 查看: 170 | 回复: 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 |
回复此楼» 猜你喜欢
0856求调剂
已经有13人回复
-
0703化学
已经有11人回复
-
329求调剂
已经有7人回复
-
321求调剂
已经有8人回复
-
0856材料化工调剂 总分330
已经有12人回复
-
309求调剂
已经有3人回复
-
070300化学354求调剂
已经有4人回复
-
330分求调剂
已经有3人回复
-
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
