While doing this derivation, the the source coordinates are mentioned as "$s$" and the coordinate of the point at which field is to be calculated is mentioned as "$r$". Kindly follow this Wikipedia link and click on the "Outline proof" under "Derivation of Gauss Law from Coulomb's law".
Finally it comes out that $$\nabla\cdot E(r)= \frac{\rho(r)}{\epsilon_0}. $$But $\rho$ is actually defined for the "$s$" coordinates and $\rho(r)$, where $r$ is the point at which electric field is calculated is 0. Here I can not understand how the $\nabla\cdot E(r)$ is equal to $\frac{\rho(r)}{\epsilon_0}$.The information about $\rho(s)$ is totally lost in the final equation. What does the Dirac delta function actually do?