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".