xn = np.linspace(xnmin, xnmax, n)