xi = np.arange(x0, xn + dx / 2., dx)