def get_outbreaks(flowmat, incidence, R0=2.5, asymf=10, attenuate=1.0):
    """
    Calculate the probabilities of outbreak for all regions

    :param flowmat: Arriving passengers row -> column
    :param incidence: fraction of infectious in the populations
    :param R0: Basic reproduction number
    :param asymf: how many asymptomatics per reported case
    :param attenuate: Attenuation factor for flow
    :return:
    """
    # Adjusting arrivals by incidence

    inflows = (flowmat.T * attenuate) @ incidence

    probs = 1 - (1 / R0) ** (inflows * 8 * asymf)

    return probs