How to treat $\nabla \cdot \left ( \frac{\vec{r}}{\left | \vec{r} \right |^{3}} \right )= 4\pi \delta \left ( \vec{r} \right )$ numerically?

Question

In electrodynamics we have $\nabla \cdot \left ( \frac{\vec{r}}{\left | \vec{r} \right |^{3}} \right )= 4\pi \delta \left ( \vec{r} \right )$, my question is how to use numerical method to treat this formula so that we can address the cases in which the singularity exists with the help of computer