The sampl.mcmc() function creates samples of the feasible region of a knapsack problem with both equalities and inequalities constraints.