The Dirac-d (d=delta) function in spherical
coordinates is
d(x - x_0) = (1/p^2) d(p - p_0)
That is, there is a 1/p^2 multiplying the usual delta functions. I've
left off the angular pieces since we assume symmetry in infinite space.
To determine G let G = H/p and solve the resulting equation for H; this
will produce the waves.