sim_prevalence function

Estimate prevalence using Monte Carlo simulation.