 of linear interpolation, fRE_H is evaluation of input Hollman energy</span>
    <span class="c">! distribution PDF at gamma supplied to function</span>

    <span class="n">rm</span><span class="o">=</span><span class="nb">sqrt</span><span class="p">((</span><span class="n">R</span><span class="o">-</span><span class="n">R0</span><span class="p">)</span><span class="o">**</span><span class="mi">2</span><span class="o">+</span><span class="p">(</span><span class="n">Z</span><span class="o">-</span><span cla