"></a>00252          coshb = - cos(beta / 2.0) <span class="comment">!-q</span>
<a name="l00253"></a>00253 
<a name="l00254"></a>00254          <span class="comment">! Initialise the plane of the dl-matrix to 0.0 for the recursion</span>
<a name="l00255"></a>00255          <span class="comment">! from l - 1 to l - 1/2.</span>
<a name="l00256"></a>00256 
<a name="l00257"></a>00257          dd(0: 2 * l + 1, 0: 2 * l + 1) = 0.0_dp
<a name="l00258"></a>00258 
<a name="l00259"></a>00259          j = 2 * l - 1
<a