(</span><span class="n">rp</span><span class="p">)</span>  <span class="kd">::</span> <span class="n">p</span>
    <span class="kt">REAL</span><span class="p">(</span><span class="n">rp</span><span class="p">)</span>  <span class="kd">::</span> <span class="n">h_j</span>

    <span class="n">gam</span><span class="o">=</span><span class="mi">1</span><span class="o">/</span><span class="nb">sqrt</span><span class="p">(</span><span class="mi">1</span><span class="o">-</span><span class="n">v</span><span class