Question: Integration involving dot products of 3-d vectors

I am trying to get Maple to understand, and evaluate the following integral analytically AND numerically.  
The original integrand to be evaluated is given by 

(w²  + m²)^(1/4) (p² + mu² )^(1/4) / ( (w² +a²)² (p² +a²) (q² + b2 ²)² ((p - q + w)² + b₁²) )

which is to be integrated simultaneously with respect to dw dp dq (where w, p, q are 3-d vectors).  
Limits of integration are (in cartesian co-ordinates) from -Infinity to + Infinity

All variables contained in the integrand are Real and >0.
a, m, mu, b₁, b₂ are constants and >0.  
dp = d³p, dq = d³q, dw = d³w
p²=p², w²=w², q²=q²

It is of course possible to help Maple along by noting that 
(p - q + w)² = p² + q² + w² - (2 p q Cos(\theta_pq) ) + (2 p w Cos(\theta_pw) ) - (2 q w Cos(\theta_wq) ) + b₁² 
where \theta_pq is the angle between vectors p and q, theta_pw is the angle between vectors p and w, etc. 

However, using the expanded form of (p - q + w)² would imply that the integral would now need to be written out in tems of 
dw_radial, dq_radial, dp_radial, d\theta_w, d\theta_q, d\theta_p, and d\phi_w, d\phi_q, d\phi_p.
Now, the integral in spherical co-ordinates would be given by the pseudocode below (note that now there are no longer any vector symbols, and the limits of integration have changed)

int ( Sin(\theta_w) d\theta_w Limits {0, pi}) int ( w² dw_radial Limits:{0, Infinity} )
int ( Sin(\theta_p) d\theta_p Limits {0, pi}) int ( p² dp_radial Limits:{0, Infinity} )
int ( Sin(\theta_q) d\theta_q Limits {0, pi}) int ( q² dq_radial Limits:{0, Infinity} )

where we have omitted the integration over \phi_w, \phi_p, \phi_q, assuming that there is no dependence on the \phi angles, resulting in all the \phi angles being integrated over, with the result being a constant multiplied to the rest of the remaining integrals. No additional information on the vectors is available, e.g. the angle between p and q, p and w, q and w, are not specified, the magnitude of the vectors is not specified. 

Any ideas if Maple can understand and execute the integral with p, q, w are input as vectors in the integrand directly, as shown earlier?
If it's easier for Maple to evaluate the integral in spherical co-ordinates, any suggestions on setting up the input, i.e. set up that the variables and constants are real, setup the limits of the theta and radial components etc. would also be appreciated.