">-</span><span class="mi">1</span><span class="p">)</span>
        <span class="n">pdf</span> <span class="o">=</span> <span class="p">(</span><span class="n">f_rolled</span> <span class="o">-</span> <span class="n">f</span><span class="p">)[</span><span class="mi">0</span><span class="p">:</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span><span class="o">/</span><span class="n">dx</span>  <span class="c1"># PDF is the derivative of CDF</span>

        <span class="k">return</s