dd_sim function

Function to simulate the diversity-dependent diversification process