xi = np.pi/2