ref: 020b9bdd65d06b73e92264c4cfe96dd73be94aed
dir: /sys/src/libc/port/frand.c/
#include <u.h>
#include <libc.h>
#define MASK 0x7fffffffL
#define NORM (1.0/(1.0+MASK))
double
frand(void)
{
double x;
do {
x = lrand() * NORM;
x = (x + lrand()) * NORM;
} while(x >= 1);
return x;
}