plt.plot(xvec, yvec3, 'g', label='CHS')