68 DOUBLE PRECISION FUNCTION dlaran( ISEED )
82 INTEGER M1, M2, M3, M4
83 parameter ( m1 = 494, m2 = 322, m3 = 2508, m4 = 2549 )
85 parameter ( one = 1.0d+0 )
88 parameter ( ipw2 = 4096, r = one / ipw2 )
91 INTEGER IT1, IT2, IT3, IT4
92 DOUBLE PRECISION RNDOUT
105 it3 = it3 + iseed( 3 )*m4 + iseed( 4 )*m3
108 it2 = it2 + iseed( 2 )*m4 + iseed( 3 )*m3 + iseed( 4 )*m2
111 it1 = it1 + iseed( 1 )*m4 + iseed( 2 )*m3 + iseed( 3 )*m2 +
113 it1 = mod( it1, ipw2 )
124 rndout = r*( dble( it1 )+r*( dble( it2 )+r*( dble( it3 )+r*
125 $ ( dble( it4 ) ) ) ) )
127 IF (rndout.EQ.1.0d+0)
THEN
double precision function dlaran(ISEED)
DLARAN