sim_n_persons function

Simulate n persons