keywordtype">real</span>(k+1,kind=dp))
<a name="l00162"></a>00162             <span class="keyword">do</span> i = 0, j - 1
<a name="l00163"></a>00163                sqrt_jmi = dsqrt(<span class="keywordtype">real</span>(j-i,kind=dp))
<a name="l00164"></a>00164                sqrt_ip1 = dsqrt(<span class="keywordtype">real</span>(i+1,kind=dp))
<a name="l00165"></a>00165                ddj = dd(i, k) / rj
<a name="l00166"></a>00166                dl(k - l, i - l) = dl(k - l, i - l) &amp;
<a name="l00167"></a>