The following equation is correct for all non-negative real numbers:
$$4\pi\delta^{(3)}(\mathbf{r})=\nabla\cdot\frac{\mathbf{r}}{r^{3}},$$ $$r\in[0,+\infty)$$
especially, when r=0, both sides give infinity.
However, the right handside can further give:
$$\nabla\cdot\frac{\mathbf{r}}{r^{3}}=\frac{\nabla\cdot\mathbf{r}}{r^{3}}+\mathbf{r}\cdot\nabla\frac{1}{r^{3}}=\frac{3}{r^{3}}-\frac{3}{r^{3}}=0$$
which is very strange, since $$4\pi\delta^{(3)}(\mathbf{r})=0,$$
I have read the proof for $r\neq0$ at Laplacians and Dirac delta functionsenter link description here, but there, this issue is not addressed.
Can anyone explain what happens here?