| 查看: 169 | 回复: 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 |
回复此楼» 猜你喜欢
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
-
300求调剂，材料科学英一数二
已经有8人回复
-
招收生物学/细胞生物学调剂
已经有5人回复
-
070305高分子化学与物理 304分求调剂
已经有7人回复
-
289求调剂
已经有13人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有9人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
调剂求院校招收
已经有5人回复
-
调剂310
已经有8人回复
