@classmethod
def fit(cls, expc, alpha, beta, scale=1e-05):
    """
        makes a model fitted to experimental concentrations expc (dict) with flux distribution given by alpha and beta
        :input expc: dict with experimental concentrations in M
        :input alpha: fraction of trioses exported as GAP
        :input beta: fraction of carbons exported as hexoses
        :input scale: optional parameter scaling the flux. Default: 1e-5, found to work.
        :returns: model with steady-state concentrations and fluxes fitted
        """
    m = cls()
    vref = m.standardCBBFluxByPartition(alpha, beta, scale=1e-05)
    res = m.stabilize(vref, expc)
    if res.success != True:
        print('Warning: Optimization not successful!')
        print(res)
    elif res.fun >= 0:
        print('Warning: Optimization did not find stable stationary state!')
        print(res)
    return m
