Zero divergence of Electric field I'm trying to rigorously derive the integral form of Gauss's law from Coulomb's law and the divergence theorem. Arrive at 
 $$ \oint\limits_{\partial V} E\cdot da = \begin{cases}
   \frac q\epsilon_{o} & \text{if $\partial V$ encloses $q$ }  \\
   0       & \text{if $\partial V$ does not enclose $q$ }
  \end{cases} 
$$
using $ E=\frac{q\hat r}{4\pi\epsilon_{o}r^2}  $ and $\oint \limits_{\partial V} E\cdot da = \int \limits_V \nabla \cdot E \space d\tau $
Proof attempt: For the case where closed surface $\partial V $ doesn't enclose any charge, I'm trying to show that $\int \limits_V \nabla \cdot E d\tau $ is 0 and then use divergence theorem to complete the proof for this case. Used the fact that $\nabla \cdot E$ is 0 everywhere except at $r=0$, since I have $\frac1{r^2}$ radial field but how do I show that $\nabla \cdot E $ is 0 through out the entire volume $V$ ?
I know this very basic, but please point out what's it, that I'm missing.
EDIT:
Now I understand that $\nabla\cdot E$ is 0 as $V$ doesn't contain $r=0$(the location of charge). But why does it turn out so, even though the field lines appear to "diverge" at each point. I'm looking for some kind of interpretation for this. Is it because "if we take a region small enough within which $E$ is almost constant, $\nabla\cdot E$ is 0".
 A: $\oint  E\cdot dS = \frac{1}{\epsilon}\int\limits_V \rho dV= \int \limits_V \nabla \cdot E \space dV $
if $\frac{1}{\epsilon}\int\limits_V \rho dV= \int \limits_V \nabla \cdot E \space dV$ then $\frac{\rho}{\epsilon} = \nabla \cdot E$
if $\rho=0$ then $\frac{\rho}{\epsilon} = 0 = \nabla \cdot E$
Is this what your looking for?
$\rho$ would be zero say, outside of a point charge.
One of the paradoxes you'll find when considering a point charge is that the divergence is zero for the field created a point charge, except at the origin in which case it is undefined. This paradox is resolved using the dirac delta function like in this excellent website which I recommend. http://farside.ph.utexas.edu/teaching/em/lectures/node30.html
In case you didn't know the divergence theorem comes from summing up the flux in infinitesimally small cubes that constitute a closed surface. The only fluxes that "survive" are the exterior ones. If $\nabla \cdot E = 0$ then this means that the flux on an infinitesimally small cube is zero. Which means that if the field lines were an actual fluid then the density of that small cube is not changing. (Because the same water flowing in flows out).
