SimTsir function

Function to simulate the stochastic TSIR