ran0.f90

来自「FORTRANvisualfortran常用数值算法集及源码」· F90 代码 · 共 14 行

FUNCTION ran0(idum)
INTEGER idum,IA,IM,IQ,IR,MASK
REAL ran0,AM
PARAMETER (IA=16807,IM=2147483647,AM=1./IM,&
                   IQ=127773,IR=2836,MASK=123459876)
INTEGER k
idum=ieor(idum,MASK)
k=idum/IQ
idum=IA*(idum-k*IQ)-IR*k
if (idum<0) idum=idum+IM
ran0=AM*idum
idum=ieor(idum,MASK)
END FUNCTION ran0