c----------------------------------------------------------------------
      function ANGINT(l1,m1,l2,m2)
c  evaluates the integral over the solid angle of a product
c  of three spherical harmonics via clebsch gordon coefficients:
c
c  for the calculation the relation:
c  \int d\omega Y_l1^m1 \cos\theta Y_l2^m2 = \sqrt{\frac{4\pi}{3}} 
c  \int d\omega Y_l1^m1 Y_1^0 Y_l2^m2 = (-1)^m2 \sqrt{\frac{2 l1 +1}
c  {2 l2+1}} <l1 1 0 0|l2 0> <l1 1 m1 0|l2 -m2>
c  is used.
c  author: Heiko Appel,  5.18.99
c      in: l1,m1,l2,m2 - quantum numbers for the spherical harmonics in
c                        the integral
c     out: return value of function: integral value for integral over 
c                                    solid angle
      implicit real*8(a-h,o-z)
      real*8 l1,l2,m1,m2