This is about finding the Green's function for the Laplacian operator under infinite boundary condition. Basically, 9.173 gives the relation between the Green's function and the unit excitation. In 3D, the unit excitation is assumed located at the center of a sphere. If you write out the gradient operator in 9.173 in the spherical coordinate on the radial component, and perform surface integration, you will get 9.174. 9.175 is the indefinite integral of 9.174.
Haberman's book on PDE Section9.5.6 has the explicit derivation.