This is to simulate the results based on a SIR model For the details of the specific model, please check the blog post here. Paramters for the simulation: