Vamos a comenzar por escribir la integral como
I=∫R3∫R3f(x+y)δ(|x|−r)δ(|y|−R)d3xd3y
donde δ es de Dirac de la función delta. El cambio de las variables a (z,y)=(x+y,y) tenemos
I=∫R3∫R3f(z)δ(|z−y|−r)δ(|y|−R)d3yd3z=∫R3f(z)(∫R3δ(|z−y|−r)δ(|y|−R)d3y)d3z
Para 0≤r≤R, se puede calcular que
Para calcular el ∫R3δ(|z−y|−r)δ(|y|−R)d3y que nos permiten utilizar coordenadas esféricas para y donde el ángulo de θ se mide desde la dirección de z, por lo que
|y|=ρ
|z−y|=√(|z|−ρcosθ)2+(ρsinθ)2=√ρ2−2ρ|z|cosθ+|z|2
Tenemos a continuación
∫R3δ(|z−y|−r)δ(|y|−R)d3y==∫∞0dρ∫π0dθ∫2π0dϕρ2sinθδ(√ρ2−2ρ|z|cosθ+|z|2−r)δ(ρ−R)==2πR2∫π0dθsinθδ(√R2−2R|z|cosθ+|z|2−r)=u:=√R2−2R|z|cosθ+|z|2=2πR2∫R+|z||R−|z||duuR|z|δ(u−r)==2πR|z|∫∞−∞duuθ((R+|z|)−u)θ(u−|R−|z||)δ(u−r)==2πRr|z|θ((R+|z|)−r)θ(r−|R−|z||)==2πRr|z|θ((R+|z|)−r)θ(r−(R−|z|))θ(r−(|z|−R))==2πRr|z|θ(|z|+(R−r))θ(|z|−(R−r))θ(R+r−|z|)=since |z|>0,R−r>0=2πRr|z|θ(|z|−(R−r))θ(R+r−|z|)
donde θ es de Heaviside de la theta de la función.
En el extremo usted consigue
I=2πRr∫|z|∈[R−r,R+r]f(z)|z|d3z