fun_stoch_sim function

Simulate for one station