plt.plot(xvec, yvec2, 'k', label='interp1d')