ing around and is unable to converge on the integral.
        The problem is actually a lack of exponent range: if we simply swap type
        double for something with a greater exponent range (an 80-bit long double
        or a quad precision type), then we can get to at least p = 0.99. If we want
        to go beyond that, or stick with type double, then we have to get smart.
      </p>
<p>
        The easiest method is to notice that for small x, then <code class="literal">tan(x) &#8773; x</code>,
     