This method uses a combination of Newton-Raphson iteration and
bisection with the given guess as a starting point. The upper and
lower bounds specify the interval in which the probability
distribution reaches the value p.