I was thinking of something similar. One has to start with a finite mass, say with a continuous mass distribution, and spherical, and use one of those theorems, Gauss or Green, to calculate the PE from its origin, to the surface, and then infinitely outward. IIRC, Brent once explained this, but I can't recall the details. AG