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