1 | c.......................................................................
   2 | c
   3 | c     tirage de nombres aleatoires en (fausse) double precision
   4 | c     
   5 | c     ce prog se trouve a l'identique sous un autre nom
   6 | c     unra.f
   7 | c
   8 | c.......................................................................
   9 |       subroutine rand_uniforme(xunra)
  10 |       implicit double precision (a-h,o-z)
  11 |       real urand
  12 | c.......................................................................
  13 |       common/comsee/iseed
  14 | c.......................................................................
  15 |       if (iseed.le.0) iseed=14396
  16 | 1000  continue
  17 |       xunra=dble(urand(iseed))
  18 |       if ((xunra.eq.0.d0).or.(xunra.eq.1.d0)) goto 1000
  19 | c.......................................................................
  20 |       return
  21 |       end
  22 | 
  23 | 
  24 | 
  25 | 
  26 | 
  27 | 
  28 | 
  29 | 
  30 | 
  31 |