#include <stdlib.h>

double rand_unif()
{
    return (1.0 * random()) / RAND_MAX;
}