prob_sec_case function

Calculate the probability of at least one secondary infection