calculate_epidemic_probability function

Main function for calculating epidemic probability