Working version with subplots and histograms

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