Hi Bill
Sorry, this is my fault. I got the initial question wrong (see my second post).
The mathematica code would be (from the corrected integral).
in[1] := $Assumptions = {r1 > 0 && r2 > 0 && gamma > 0};
Integrate[ Integrate[ s1 s2 Exp[-gamma Sqrt[s1*s1 + s2*s2]] BesselJ[0, s1 r1]...