Computes the quantile (inverse CDF) for the given probability p.
Uses Newton's method seeded with the distribution mean as the initial guess.
the value x >= 0 such that cdf(x) = p.
cdf(x) = p
the cumulative probability, must be in 0, 1.