quadgk is fine for vector-valued functions (and in fact any integrand type supporting +, -, *real, and norm, and can also integrate over infinite intervals and contours in the complex plane.
But for multidimensional integrals you should use Cuba.jl or Cubature.jl or GSL.jl or similar.
