Working version with subplots, need to test further
Working version with subplots and histograms

 - adds print statements and function to calculate
   num patients cured given a certain timestep