1 2 3 4 5 6 7 8 9 10 11
#include "internal.h" /* Returns a random float64 between 0 and 1 */ double rng·exponential(double lambda) { double f; f = rng·random(); return -log(1 - f)/lambda; }