 dlj * coshb
<a name="l00142"></a>00142                dd(i + 1, k) = dd(i + 1, k) &amp;
<a name="l00143"></a>00143                     - sqrt_ip1 * sqrt_jmk * dlj * sinhb
<a name="l00144"></a>00144                dd(i, k + 1) = dd(i, k + 1) &amp;
<a name="l00145"></a>00145                     + sqrt_jmi * sqrt_kp1 * dlj * sinhb
<a name="l00146"></a>00146                dd(i + 1, k + 1) = dd(i + 1, k + 1) &amp;
<a name="l00147"></a>00147                     + sqrt_ip1 * sqrt_kp1 * dlj * coshb
<a name="l0014