sim_abundance function

Simulate basic population dynamics model