usr/src/lib/libF77/rand_.c

/*
Uniform random number generator.  Code courtesy of Bob Morris.
Linear congruential generator, suitable for 32 bit machines;
multiplication is mod 2**31
*/

static  long    randx = 1;

srand_(x)       /* subroutine to set seed */
long *x;
{
randx = *x;
}


double rand_()
{
double ldexp();
return(ldexp((double)(((randx = randx*1103515245 + 12345)>>7) & 077777777), -24));
}