norm(i))
<a name="l01116"></a>01116         <span class="keyword">do</span> j=1, i
<a name="l01117"></a>01117           dl(i - lp1, j - lp1)= dl(i - lp1, j - lp1)*renorm
<a name="l01118"></a>01118         <span class="keyword">end do</span>
<a name="l01119"></a>01119       <span class="keyword">end do</span>
<a name="l01120"></a>01120       <span class="keyword">do</span> i=l+2,2*l+1
<a name="l01121"></a>01121         expo=log_first_row(i)-lrenorm(i)
<a name="l01122"></a>01122         renorm=sign(i)*dexp(lo