tsirSim function

Function to simulate the stochastic TSIR