sim_cf_dens function

Simulate counterfactual densities