s***e 发帖数: 911 | 2
我把程序直接给你? 很短. 这个程序很好用. 其中有些注释句很多余,是以前我调程序
干的,没删掉. 是单精度实行.反正是random number, 没必要双精度. 主程序掉用强行
转一下就成.
FUNCTION ran2(idum)
c REAL FUNCTION ran2(idum)
INTEGER idum,IM1,IM2,IMM1,IA1,IA2,IQ1,IQ2,IR1,IR2,NTAB,NDIV
REAL ran2,AM,EPS,RNMX
* REAL*8 FUNCTION ran2(idum)
c REAL FUNCTION ran2(idum)
PARAMETER (IM1=2147483563,IM2=2147483399,AM=1./IM1,IMM1=IM1-1,
*IA1=40014,IA2=40692,IQ1=53668,IQ2=52774,IR1=12211,IR2=3791,
*NTAB=32,NDIV=1+IMM1/NTAB,EPS=1.2e-7,RNMX |
|