Given a finite Borel measure μ on R n and basic semi-algebraic sets Ω\_i ⊂ R n , i = 1,. .. , p, we provide a systematic numerical scheme to approximate as closely as desired μ(\cup\_i Ω\_i), when all moments of μ are available (and finite). More precisely , we provide a hierarchy of semidefinite programs whose associated sequence of optimal values is monotone and converges to the desired value from above. The same methodology applied to the complement R n \ (\cup\_i Ω\_i) provides a monotone sequence that converges to the desired value from below. When μ is the Lebesgue measure we assume that Ω := \cup\_i Ω\_i is compact and contained in a known box B and in this case the complement is taken to be B \ Ω. In fact, not only μ(Ω) but also every finite vector of moments of μ\_Ω (the restriction of μ on Ω) can be approximated as closely as desired, and so permits to approximate the integral on Ω of any given polynomial.