sim_div function

Simulation of dividends