053       xllp1=<span class="keywordtype">real(l*(l+1)</span>,kind=dp)
<a name="l01054"></a>01054   
<a name="l01055"></a>01055       <span class="comment">! Initialising coefficients cp(m)= cplus(l-m).</span>
<a name="l01056"></a>01056   
<a name="l01057"></a>01057       <span class="keyword">do</span> m=1,l
<a name="l01058"></a>01058         xm=<span class="keywordtype">real</span>(l-m,kind=dp)
<a name="l01059"></a>01059         cpi(m)=2.d0/dsqrt(xllp1-xm*(xm+1))
<a name="l01060"></a>01060         cp(m)=1