I have derived the Green's function for the 3D wave equation as $$G(x,y,t,\tau)=\frac{\delta\left(|x-y|-c(t-\tau)\right)}{4\pi c|x-y|}$$ and I'm trying to use this to solve $$u_{tt}-c^2\nabla^2u=0 \hspace{10pt}u(x,0)=0\hspace{10pt} u_t(x,0)=f(x)$$ but I'm not sure how to proceed. I tried converting the problem to a forced wave equation with homogenous boundary conditions by setting $u'=u-tf(x)$ so $u'_{tt}-c^2\nabla^2u'=c^2t\nabla^2f(x)$ and then integrating to find $$u=f(x)t+\int_0^\infty\int_{\mathbb{R}^3}G(x,y,t,\tau)c^2\tau\nabla^2f(y)$$ but when I evaluate this integral I don't seem to recover the causal structure of the wave equation. I'm meant to find that $u=\langle f\rangle t$ where $\langle f\rangle$ is the average value of $f$ on a sphere of radius $ct$ around $x$.
Firstly, is this the right approach to using Green's functions here? And if so how can I recover the causal structure of the problem?