6"></a>00136             <span class="keyword">do</span> i = 0, j - 1
<a name="l00137"></a>00137                sqrt_jmi = dsqrt(<span class="keywordtype">real</span>(j-i,kind=dp))
<a name="l00138"></a>00138                sqrt_ip1 = dsqrt(<span class="keywordtype">real</span>(i+1,kind=dp))
<a name="l00139"></a>00139                dlj = dl(k - (l - 1), i - (l - 1)) / rj
<a name="l00140"></a>00140                dd(i, k) = dd(i, k) &amp;
<a name="l00141"></a>00141                     + sqrt_jmi * sqrt_jmk *