def model(
    time: float,
    variables: float, 
):
    e_coli, c_gluta, enterobactin = variables
    mu_e = 0.4
    mu_c = 0.3
    a_e = 6
    K_e = 0.5
    K_c = 0.5
    theta = 0.001
    r_prod = 0.2
    r_cons_e = 1
    r_cons_c = 1
    a_c = (10) - (a_e)
    uptake_E_growth = (a_e * enterobactin) / (K_e + enterobactin)
    uptake_C_growth = (a_c * enterobactin) / (K_c + enterobactin)
    cons_term_E = mu_e * (a_e * enterobactin) / (K_e + a_e * enterobactin) * e_coli
    cons_term_C = mu_c * (a_c * enterobactin) / (K_c + a_c * enterobactin) * c_gluta
    dEdt = mu_e * uptake_E_growth * e_coli
    dCdt = (mu_c * uptake_C_growth * c_gluta) - (theta * c_gluta * c_gluta)
    dBdt = ((r_prod * enterobactin) - (r_cons_e * cons_term_E)) - (r_cons_c * cons_term_C)
    de_colidt = +dEdt
    dc_glutadt = +dCdt
    denterobactindt = +dBdt
    return [de_colidt, dc_glutadt, denterobactindt]

y0 = {"e_coli": 5, "c_gluta": 5, "enterobactin": 1}