fy0 = np.interp(y, x, fx)