span>
<a name="l00154"></a>00154 
<a name="l00155"></a>00155          dl(- l: l, - l: l) = 0.0_dp
<a name="l00156"></a>00156 
<a name="l00157"></a>00157          j = 2 * l
<a name="l00158"></a>00158          rj = <span class="keywordtype">real</span>(j)
<a name="l00159"></a>00159          <span class="keyword">do</span> k = 0, j - 1
<a name="l00160"></a>00160             sqrt_jmk = dsqrt(<span class="keywordtype">real</span>(j-k,kind=dp))
<a name="l00161"></a>00161             sqrt_kp1 = dsqrt(<span class="