 * l + 1, 0: 2 * l + 1) = 0.0_dp
<a name="l00379"></a>00379 
<a name="l00380"></a>00380          j = 2 * l - 1
<a name="l00381"></a>00381          rj = <span class="keywordtype">real</span>(j)
<a name="l00382"></a>00382          <span class="keyword">do</span> k = 0, j - 1
<a name="l00383"></a>00383             <span class="keyword">do</span> i = 0, l <span class="comment">!j/2 + 1</span>
<a name="l00384"></a>00384                dlj = dl(k - (l - 1), i - (l - 1)) / rj
<a name="l00385"></a>00385            