A general and rigorous method to deal with singularities at the origin of a polar coordinate system is presented. Its power derives from a clear distinction between the radial distance and the radial coordinate variable, which makes that all delta functions and their derivatives are automatically generated, and ensures that the Gauss theorem is…