def METHOD_NAME():
    model = get_model(io.StringIO(par))
    toas = pint.simulation.make_fake_toas_uniform(
        54000, 56000, 100, model=model, obs="gbt"
    )
    dmxmodel = get_model(io.StringIO(par + dmxpar))
    indices = model.add_DMX_ranges([54500, 55500], [55000, 56000], dmxs=[2, 3])
    assert np.all(np.array(indices) == np.array([2, 3]))
    assert np.all(
        model.components["DispersionDMX"].dmx_dm(toas)
        == dmxmodel.components["DispersionDMX"].dmx_dm(toas)
    )