dsin(ang)
<a name="l01028"></a>01028       si=1.d0/s
<a name="l01029"></a>01029       t=dtan(-ang/2.d0)
<a name="l01030"></a>01030       c2=dcos(ang/2.d0)
<a name="l01031"></a>01031       omc=1.d0-c
<a name="l01032"></a>01032   
<a name="l01033"></a>01033       <span class="keyword">do</span> i=1,2*l+1
<a name="l01034"></a>01034         lrenorm(i)=0.d0
<a name="l01035"></a>01035       <span class="keyword">end do</span>
<a name="l01036"></a>01036   
<a name="l01037"></a>01037       <span class="comment">! C