quantile
Computes the quantile (inverse CDF) for the given probability p.
Uses ln1p for numerical stability when p is close to zero.
Return
the value x at which cdf(x) = p. Returns Double.POSITIVE_INFINITY for p = 1.
Parameters
p
the cumulative probability, must be in [0, 1].