 dd(0: 2 * l + 1, 0: 2 * l + 1) = 0.0_dp
<a name="l00130"></a>00130 
<a name="l00131"></a>00131          j = 2 * l - 1
<a name="l00132"></a>00132          rj = <span class="keywordtype">real</span>(j)
<a name="l00133"></a>00133          <span class="keyword">do</span> k = 0, j - 1
<a name="l00134"></a>00134             sqrt_jmk = dsqrt(<span class="keywordtype">real</span>(j-k,kind=dp))
<a name="l00135"></a>00135             sqrt_kp1 = dsqrt(<span class="keywordtype">real</span>(k+1,kind=dp))
<a name="l0013