x=0:5;
y=cos(x);
dy=-sin(x);
n=length(x);
xi=linspace(x(1),x(n),100);
yi=hermfun(x,y,dy,xi);
plot(x,y,'o',xi,yi)
