emonstration, first we need some includes to access the negative
          binomial distribution (and the binomial, beta and gamma distributions too).
        </p>
<p>
          To demonstrate the use with a high precision User-defined floating-point
          type <code class="computeroutput"><span class="identifier">cpp_dec_float</span></code> we also
          need an include from Boost.Multiprecision.
        </p>
<pre class="programlisting"><span class="preprocessor">#include</span> <span class="specia