Numerical integration rules for multivariate inversions