Given a semialgebraic set S.Rn, we provide a numerical approximation procedure that provides upper and lower bounds on .(S), for measures . that satisfy some given moment conditions. The bounds are obtained as solutions of positive semidefinite programs that can be solved via standard software packages.