anchor.diff:          else
             wb = (l_max+m_mid)*(l_max+m_mid+1)/2
           end if
  -C        write (6,*) wa,wb
           result = sqrt(wb)/wa
   *
   *   General code for l_min > 1
   *
         else
  +*
  +*     (a+b-c-1)!!(b+c-a-1)!!(c+a-b-1)!!/(a+b+c+1)!!
  +*
           l_sum = l1+l2+l3
  -        ll_max = l_sum
  +        ll_max = 2*l_max
           ll_min = 2*l_min
           ll_mid = 2*l_mid
           m = 1
  
positive.code: C$Id: ecp_3j_prod.f,v 1.4 1997-05-21 09:07:49 d3e129 Exp $
  ************************************************************************
  *                                                                      *
        subroutine ecp_3j_prod (l1,l2,l3,m1,m2,result)
  
negative.code:         do l = 3,l_sum-ll_mid-1,2
            m = m*l
          end do
          do l = 3,l_sum-ll_max-1,2
            m = m*l
          end do
          wb = m
          result = wb/wa
  *
  *     [(c+gamma)!(c-gamma)!/(a+alpha)!(a-alpha)!(b+beta)!(b-beta)!]^1/2
  *
          m = 1
          n = 1
  