anchor.diff:          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        write (6,*) wb,wa,result
  +*
  +*     [(c+gamma)!(c-gamma)!/(a+alpha)!(a-alpha)!(b+beta)!(b-beta)!]^1/2
   *
           m = 1
           n = 1
  
positive.code: C$Id: ecp_3j_prod.f,v 1.2 1996-10-11 16:47:45 d3e129 Exp $
  ************************************************************************
  *                                                                      *
        subroutine ecp_3j_prod (l1,l2,l3,m1,m2,result)
  
negative.code:         else
            wb = (l_max+m_mid)*(l_max+m_mid+1)/2
          end if
          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 = 2*l_max
          ll_min = 2*l_min
          ll_mid = 2*l_mid
          m = 1
  