xo = np.hstack((xo, x1 + xo[-1]))