Question: why divergence of E field is not zero?

r := sqrt(x^2+y^2+z^2);
divE := diff(e*x/(4*Pi*r^3), x)+diff(e*x/(4*Pi*r^3), y)+diff(e*x/(4*Pi*r^3),z);
simplify(divE);

Please Wait...