sample_marginal function

Exact independent samples from an approximate posterior distribution