calc_epidemic_onset function

Calculate epidemic onset