The sampl.mcmc() function creates samples of the feasible region of a knapsack problem with both equalities and inequalities constraints.
conda install r::r-knapsacksampling