x=1:10;
y=cos(x);
n=length(x);
xi=linspace(x(1),x(n),100);
for i=1:length(xi)
    tmp=nevilleQ(x,y,xi(i));
    yi(i)=tmp(n,n);
end
plot(x,y,'o',xi,yi)
