adding seasonality to betas in four-serotype model closes #2