expdev.for

来自「Numerical Recipes一书中例子的源码所用到的函数集,William」· FOR 代码 · 共 11 行

      FUNCTION expdev(idum)
      INTEGER idum
      REAL expdev
CU    USES ran1
      REAL dum,ran1
1     dum=ran1(idum)
      if(dum.eq.0.)goto 1
      expdev=-log(dum)
      return
      END