/* 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)); }