# Divergence of inverse square

Hello,

I tried to calculate the divergence of 1/r^2*e_r, but I got zero, where it should have been 4*pi*dirac_delta(r*e_r)

```
E.<r,th,ph>=EuclideanSpace(coordinates="spherical")
v=E.vector_field([1/r^2,0,0])
v.div().display()
```

Results is

```
E^3 --> R
(r, th, ph) |--> 0
```

Can this be considered as a bug?