def update_toestand(A, x, resistent=[]):
    """
    Voor een gegeven verbindingsmatrix A en een binaire toestandsvector x
    (0: vatbaar/resitent, 1: geinfecteerd) op tijdstip t, bereken de toestandsvector op
    tijdstip t+1. Optioneel kan je een lijst met indices van resistente individuen
    geven.
    """
    # update x
    x[:] = (x + A @ x) > 0
    # resitente worden niet ziek
    if resistent:
        x[resistent] = 0
    return x