      lambda*cpi(m)*dl(index - lp1, m - lp1)-cp2(m) &amp;
<a name="l01080"></a>01080               *dl(index - lp1, m-1 - lp1)
<a name="l01081"></a>01081             <span class="keyword">if</span>(dl(index - lp1, m+1 - lp1).gt.big) <span class="keyword">then</span>
<a name="l01082"></a>01082               lrenorm(index)=lrenorm(index)-lbig
<a name="l01083"></a>01083               <span class="keyword">do</span> im=1,m+1
<a name="l01084"></a>01084                 dl(index - lp1, im - lp1)=dl(index - lp1, im