o fill columns down to diagonals.</span>
<a name="l01069"></a>01069   
<a name="l01070"></a>01070       <span class="keyword">do</span> index= 2,l+1
<a name="l01071"></a>01071         dl(index - lp1, 1 - lp1)=1.d0
<a name="l01072"></a>01072         lambda=(<span class="keywordtype">real</span>(l+1,kind=dp)*omc-<span class="keywordtype">real</span>(index,kind=dp)+c)*si
<a name="l01073"></a>01073         dl(index - lp1, 2 - lp1)=lambda*dl(index - lp1, 1 - lp1)*cpi(1)
<a name="l01074"></a>01074         <span c