simulate_flu_fit function

Fitting a SIR-type model to flu data