anchor.diff:           if(l.eq.1.and.m.eq.0)ylm_real=dsqrt(3.d0/fourpi)*zchap
            if(l.eq.1.and.m.eq.-1)ylm_real=dsqrt(3.d0/fourpi)*ychap
   
  -         if(l.eq.2.and.m.eq.2)ylm_real=dsqrt(15.d0/16.d0/pi)*(xchap**2-ychap**2)
  +         if(l.eq.2.and.m.eq.2)ylm_real=dsqrt(15.d0/16.d0/pi)*(xchap*xchap-ychap*ychap)
            if(l.eq.2.and.m.eq.1)ylm_real=dsqrt(15.d0/fourpi)*xchap*zchap
  -         if(l.eq.2.and.m.eq.0)ylm_real=dsqrt(5.d0/16.d0/pi)*(-xchap**2-ychap**2+2.d0*zchap**2)
  +         if(l.eq.2.and.m.eq.0)ylm_real=dsqrt(5.d0/16.d0/pi)*(-xchap*xchap-ychap*ychap+2.d0*zchap*zchap)
            if(l.eq.2.and.m.eq.-1)ylm_real=dsqrt(15.d0/fourpi)*ychap*zchap
            if(l.eq.2.and.m.eq.-2)ylm_real=dsqrt(15.d0/fourpi)*xchap*ychap
   
  
positive.code:         DIMENSION PM(0:100,0:100)
          MM=100
          pi=dacos(-1.d0)
          iabs_m=iabs(m)
          if(iabs_m.gt.l)stop 'm must be between -l and l'
          factor= dsqrt( ((2*l+1)*fact(l-iabs_m))/(4.d0*pi*fact(l+iabs_m)) )
          if(dabs(x).gt.1.d0)then
           print*,'pb. in ylm_no'
           print*,'x=',x
  
negative.code:         two_qkmp1 = two_qkmp1-2.d0
          qk = qk-1.d0 
        enddo
        inverses(q) = a_over_b_square/(dble(q+n+q+n+3) * dble(q+1)) 
  !      do k=0,q
  !        sum=sum+s_q_k
  !        s_q_k = a_over_b_square * ( dble(2*(q-k+m)+1)*dble(q-k)/(dble(2*(k+n)+3) * dble(k+1)) ) * s_q_k
  