92"></a>01092   
<a name="l01093"></a>01093       <span class="keyword">do</span> index= l+2,2*l
<a name="l01094"></a>01094         dl(index - lp1, 1 - lp1)=1.d0
<a name="l01095"></a>01095         lambda=(<span class="keywordtype">real</span>(l+1,kind=dp)*omc-<span class="keywordtype">real</span>(index,kind=dp)+c)/s
<a name="l01096"></a>01096         dl(index - lp1, 2 - lp1)=lambda*dl(index - lp1, 1 - lp1)*cpi(1)
<a name="l01097"></a>01097         <span class="keyword">if</span>(index.lt.2*l) <span class="k