#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;
}