adding the script for fitting a harmonic sir model to flu data