 Solution

On an outward-oriented sphere of radius $R>0$, the electric field is normal with constant magnitude $A/R^2$. Thus its flux isIt is independent of $R$. There is no contradiction with solution: the classical divergence theorem cannot be applied to a full ball while pretending that $E$ is smooth at its center. On a punctured ball the outer flux and the inward-normal inner flux cancel, in agreement with zero ordinary divergence there.

The three-dimensional Dirac delta distribution ${\delta }^{(3)}(r)$ is supported at the origin and satisfies the sampling identityfor smooth compactly supported test functions $\phi$. It is not an ordinary function with a finite value at zero; its total integral is one, and translating it moves the point at which the test function is sampled.

The globally valid source equation isFor a direct weak verification, excise a ball of radius $\epsilon$ and integrate $-E\cdot \nabla \phi$ outside it. Since the ordinary divergence vanishes there, the result is $A{\int }_{S^2}\phi (\epsilon \hat{r})d\Omega \to 4\pi A\phi (0)$. This is precisely the displayed distributional divergence and restores the divergence theorem with the point charge included.