Estamos dado que el $\vec F$ es continuamente diferenciable, es decir, $\vec F \in C^1(D, \Bbb R^3)$.
Luego, si no existe $p \in D$ con
$\nabla \cdot \vec F(p) \ne 0, \tag 1$
debemos tener bien
$\nabla \cdot \vec F(p) > 0 \tag 2$
o
$\nabla \cdot \vec F(p) < 0; \tag 3$
si (2) se une, se sigue de la continuidad de $\nabla \cdot \vec F$ que
existe alguna $\epsilon > 0$ tal que
$\nabla \cdot \vec F(x) > 0 \tag 4$
para
$x \in B(p, \epsilon) = \{x \in D \mid \Vert x - p \Vert < \epsilon \}; \tag 5$
entonces
$\displaystyle \int_{B(p, \epsilon)} \nabla \cdot \vec F(x) \; dV > 0; \tag 6$
por el teorema de la divergencia (6) implica
$\displaystyle \int_{\partial B(p, \epsilon)} \vec F \cdot d \vec S > 0, \tag 7$
lo que se contradice con el dado hipótesis. Por lo tanto podemos descartar (2); del mismo modo, si (3) se mantiene, podemos por un argumento similar a la conclusión de que
$\displaystyle \int_{\partial B(p, \epsilon)} \vec F \cdot d \vec S < 0, \tag 8$
lo cual también contradice la suposición de que
$\displaystyle \int_{\partial B(p, \epsilon)} \vec F \cdot d \vec S = 0; \tag 9$
por lo tanto debemos tener
$\nabla \cdot \vec F(p) = 0 \tag{10}$
para todos los $p \in D$.