x1 = dx * np.arange(1, np.ceil((max_x - xo[-1]) / dx) + 1)