def METHOD_NAME(m, dynamic=True, ntfe=None):
    if ntfe is not None and not dynamic:
        raise RuntimeError("Cannot provide ntfe to initialize steady model")
    elif dynamic and ntfe is None:
        ntfe = 10
    if dynamic:
        disc = pyo.TransformationFactory("dae.finite_difference")
        disc.apply_to(m, wrt=m.time, nfe=ntfe, scheme="BACKWARD")
    t0 = m.time.first()
    # Fix inputs
    m.conc_in[:, "A"].fix(5.0)
    m.conc_in[:, "B"].fix(0.01)
    m.flow_in[:].fix(1.0)
    m.flow_in[t0].fix(0.1)
    if dynamic:
        # Fix initial conditions if dynamic
        m.conc[t0, "A"].fix(1.0)
        m.conc[t0, "B"].fix(0.0)