probability_epidemic function

Calculate the probability a disease will cause an outbreak based on R, k and initial cases