pleio_simulate function

Create simulations