Most of the books about electromagnetism prove Gauss' law for a point charge in vacuum:
$$ \Phi = \int_{\Sigma} \mathbf{E} \centerdot d \mathbf{S} = q/\epsilon_0 $$
and then simply state that for a continuous charge distribution the charge is $$ q= \int_{V'} \rho (\mathbf{r'}) dV' $$ and thus the application of the divergence theorem gives the differential form of Gauss' law: $$ {\rm div} (\mathbf{E}) = \rho/\epsilon_0 .$$
But it is always true that given any (integrable) charge density $ \rho(\mathbf{r'}) $ distributed over an arbitrary volume $V'$, such that the produced electric field is:
$$ \mathbf{E}(\mathbf{r})=1/{4\pi \epsilon_0}\int_{V'} \rho(\mathbf{r'})\frac{(\mathbf{r}-\mathbf{r'})}{\mid \mathbf{r} - \mathbf{r'}\mid^3}dV', $$
that $$ {\rm div} (\mathbf{E})=\frac{\rho}{\epsilon_0} ~?$$
How could this be proved rigorously?