I am trying to formally learn electrodynamics on my own (I only took an introductory course). I have come across the differential form of Gauss's Law.
$$ \nabla \cdot \mathbf E = \frac {\rho}{\epsilon_0}.$$
That's fine and all, but I run into what I believe to be a conceptual misunderstanding when evaluating this for a point charge.
I know the math looks better in spherical coordinates, but I will be using Cartesian.
So when I calculate the divergence I obtain:
$$ \nabla \cdot \mathbf E = \nabla \cdot kQ\langle\frac{x}{(x^2+y^2+z^2)^{\frac{3}{2}}},\frac{y}{(x^2+y^2+z^2)^{\frac{3}{2}}},\frac{z}{(x^2+y^2+z^2)^{\frac{3}{2}}}\rangle = \frac{-3(x^2+y^2+z^2)}{(x^2+y^2+z^2)^{\frac{5}{2}}}+\frac{3}{(x^2+y^2+z^2)^{\frac{3}{2}}}.$$
This can further be simplified:
$$\frac{-3(x^2+y^2+z^2)}{(x^2+y^2+z^2)^{\frac{5}{2}}}+\frac{3}{(x^2+y^2+z^2)^{\frac{3}{2}}} = \frac{3}{(x^2+y^2+z^2)^{\frac{3}{2}}}-\frac{3}{(x^2+y^2+z^2)^{\frac{3}{2}}} = \frac{3-3}{(x^2+y^2+z^2)^{\frac{3}{2}}}.$$
Now instinctively I would say that 3-3 is zero and then the while thing is zero everywhere. I am confused as to why (purely mathematically) this expression is not equal to zero at the origin. I completely understand why it physically has to be that way. And I also understand that it is modeled with the delta dirac function. But what (again, mathematically) is stopping me from saying that equation is just zero even at the origin?