R3d = InterpolatingFunction.Spline(r,R3dVal) -- calculate Slater integrals for 3d shell Fk = GetSlaterIntegrals({"3d"},{R3d}) F0dd = Fk["3d 3d 3d 3d"][0]*EnergyUnits.Ha.value F2dd = Fk["3d 3d 3d 3d"][2]*EnergyUnits.Ha.value F4dd = Fk["3d 3d 3d 3d"][4]*EnergyUnits.Ha.value