Integrate[1/(a^2+x^2+y^2)^(3/2),{x,0,a},{y,0,a}]