fo = np.interp(xi, xo, fo)