<a name="l01018"></a>01018         <span class="keyword">end do</span>
<a name="l01019"></a>01019         return      <span class="comment">! ** Exit routine</span>
<a name="l01020"></a>01020       <span class="keyword">end if</span>
<a name="l01021"></a>01021   
<a name="l01022"></a>01022       lp1 = l + 1
<a name="l01023"></a>01023       big_const = 1.0d150
<a name="l01024"></a>01024 
<a name="l01025"></a>01025       ang=beta
<a name="l01026"></a>01026       c=dcos(ang)
<a name="l01027"></a>01027       s=