def parameters_adp_atp_sum_dep(sum_atp_adp):
    # Cloutier steady state
    glc_ss = 0.16
    atp_ss = 2.2
    adp_ss = sum_atp_adp - atp_ss
    gly_ss = 0.05
    pyr_ss = 0.12
    lac_ss = 0.89


    f_glc = 0.048
    lac_degr = 0.0969


    b_0 = sum_atp_adp / 2
    c_0 = sum_atp_adp / 2


    # How we compute the reaction rate from the ODE and Cloutier's steady state:
    K_hxk = f_glc / (glc_ss * atp_ss ** 2)
    K_pyrk = 2 * f_glc / (adp_ss ** 2 * gly_ss)
    K_ldh = lac_ss * lac_degr / pyr_ss
    K_Mito = (2 * f_glc - lac_degr * lac_ss) / (pyr_ss * adp_ss ** 28)
    K_act = (58 * f_glc - 28 * lac_degr * lac_ss) / atp_ss

    return(b_0, c_0, K_hxk, K_pyrk, K_ldh, K_Mito, K_act, f_glc, lac_degr)