self.ensure(
            x1 >= self.map_ymin and x1 <= self.map_ymax,
            "The initial condition x1 should be in [%g, %g].",
            self.map_ymin,
            self.map_ymax,