a>01044         log_first_row(i)=log_first_row(i-1) &amp;
<a name="l01045"></a>01045           +dlog(ratio)+dlog(dabs(t))
<a name="l01046"></a>01046         sign(i)=sign(i-1)*t/dabs(t)
<a name="l01047"></a>01047       <span class="keyword">end do</span>
<a name="l01048"></a>01048   
<a name="l01049"></a>01049       big=big_const
<a name="l01050"></a>01050       lbig=dlog(big)
<a name="l01051"></a>01051       bigi=1.d0/big_const
<a name="l01052"></a>01052       bigi2=1.d0/big_const**2
<a name="l01053"></a>01