proves the gauss theorem. First they prove that there has to be charge allocated in the origo of a function and that this is where you get the gauss formula and therefore you can integrate over any shape and obtain the same value. It is expressed as
$\nabla \cdot E= \frac{q_i}{\varepsilon_o}\delta(r)$ (208)
where
$\delta(r)$
Is a funcion that is 0 everywhere except at the origin where it is 1 to make the fact that it is 0 everywhere than in the origo as proved by equation (192)
$\nabla \cdot E=\frac{\partial E_x}{\partial x}+\frac{\partial E_y}{\partial y }+\frac{\partial E_z}{\partial z}=\frac{q}{4\pi\varepsilon_o}\frac{3r^2-3x^2-3y^2-3z^2}{r^5}=0$ (192)
and (192) is only undefined at the origo In (210) they extend this to many charges by the formula
$\nabla \cdot E= \sum_{i=1}^{N}\frac{q_i}{\varepsilon_o}\delta(r-r_i)$ (210)
My problem is that I don't get why the spatial arrangment of the additional charges other places then at the origo don't matter for the electrical dot product. They just extend it to (210) from (208) without a mathematically justification. I wanted a mathematically justification for this generalization. I wanted to ask this to an expert but I am not sure how to pose the question clever enough. Can you help me formulate a mathematical problem for this extention?