def METHOD_NAME(self, mc_method):
    evol = MCSolver.avail_integrators()[mc_method](self.se_system, {})
    state = qutip.basis(2,0).data
    evol.set_state(0, state)
    t = 0
    for i in range(1, 21):
        t_target = i * 0.05
        while t < t_target:
            t_old, y_old = evol.get_state()
            t, state = evol.mcstep(t_target)
            assert t <= t_target
            assert t > t_old
            assert_allclose(self._analytical_se(t),
                            state.to_array()[0, 0], atol=2e-5)
            t_back = (t + t_old) / 2
            t_got, bstate = evol.mcstep(t_back)
            assert t_back == t_got
            assert_allclose(self._analytical_se(t),
                            state.to_array()[0, 0], atol=2e-5)
            t, state = evol.mcstep(t)