For the program pointc alternatively, the radial matrix elements can be given:
#--------------------------------------------------------- # Radial Matrix Elements (e.g. Abragam Bleaney 1971 p 399) #--------------------------------------------------------- #<r^2> in units of a0^2 a0=0.5292 Angstroem R2=1.114 #<r^4> in units of a0^4 a0=0.5292 Angstroem R4=2.91 #<r^6> in units of a0^6 a0=0.5292 Angstroem R6=15.03