rand_unit function

Pseudo-random numbers on [0, 1) based on a linear congruential generator