classical inverse-distance potential in 3D in spherical coordinates