1      Subroutine eumac(alpha,r,w,nr,rm)
2
3C$Id$
4      implicit none
5
6      double precision r(*) ! grid pts coord [output]
7      double precision w(*) ! grid pts weights [output]
8      double precision alpha !   [input]
9      double precision rm    !   [input]
10c
11      integer nr,i
12      double precision fmn,qi,ri,wi
13c
14c Definition of Euler-Maclaurin numerical quadrature points and weights
15c for radial integrals.
16c
17c***************************************************************************
18      nr = nr+1
19
20      fmn = dble(rm)/dble(nr)
21
22      do 10 i = 1,nr-1
23
24       qi = dble(i)/dble(nr)
25       ri = alpha*(qi/(1.D0 - qi))**rm
26       wi = fmn*(ri*ri*ri)/(qi*(1.D0 - qi))
27
28       r(i) = ri
29       w(i) = wi
30
31   10 continue
32
33      nr = nr-1
34
35      end
36