/* a proper drand48 implementation, at least if you're win32 and RAND_MAX==32768 */
double
drand48()
{
        unsigned short rv[4];
        rv[0]  = rand();
        rv[1]  = rand();
        rv[2]  = rand();
        rv[3]  = 0x3FF0 | (rv[0] & 0xF);
        rv[0] &= 0xFFF0;

        return *((double *)rv) - 1.0;
}