xi = np.linspace(0, 1, 45)