plt.plot(xvec, yvec, 'b', label='Pchip')