My question is about distributional meanings in PDEs. My specific question is at the bottom, but I'd be interested in a bit of general theory (even if this is a link to a particularly good set of notes). In my lectures, not much is said about them and I don't exactly understand what's going on!

(I believe that $|x|^{-1}$ is $\| x \|^{-1}$, not $|x_1|^{-1}$; this makes sense in the latter part of the question.)
I think that, given $T \in \mathcal S'(\Bbb R^3)$ and letting $f:\Bbb R^3 \rightarrow \Bbb R$ via $f(x) = \| x \|^{-1}$, we obtain $$-T(\Delta f) = 4 \pi \delta_0 \ \iff T({2 \over {\|x\|^3}}(3\|x\| + 4) = 4 \pi \delta_0,$$ by taking the Laplacian $\Delta$ of $f$. Is this correct?
If it is, then I still don't know how to do the next bit. By taking the divergence of $R_\epsilon := B_R \setminus B_\epsilon$, we should get $0$ on the LHS for all $\epsilon > 0$, but then $4\pi$ in the limit; how do I take the divergence with the $T$ there? Do I use the integral notation and then try using one of Green's identities?