A flexible numerical integration method is proposed for the computation of moments of a multivariate posterior density with different tail properties in different directions. The method (called mixed integration) amounts to a combination of classical numerical integration and Monte Carlo integration. Mixed integration is parsimonious in the sense that the method makes use of the same parameters as the more restrictive multivariate normal importance function.