 k + 1) = dd(i + 1, k + 1) &amp;
<a name="l00392"></a>00392                     + sqrt_tbl(i+1) * sqrt_tbl(k+1) * dlj * coshb
<a name="l00393"></a>00393             <span class="keyword">end do</span>
<a name="l00394"></a>00394          <span class="keyword">end do</span>
<a name="l00395"></a>00395 
<a name="l00396"></a>00396          <span class="comment">! Having constructed the d^(l+1/2) matrix in dd, do the second</span>
<a name="l00397"></a>00397          <span class="comment">! half-step recursion fro