int( (a^2+z^2)^(-3/2),z)