Help for this page
#define Drand01() ((rand() & 0x7FFF) / (double) ((unsigned long)1 << +15)) ^^^^
#define Drand01() drand48()