Generates uniformly distributed pseudorandom-number sequences.
    Returns 48-bit nonnegative, double-precision, floating-point
    values.

    Format

      #include  <stdlib.h>

      double erand48  (unsigned short int xsubi[3]);