function y=myeval(xi,k,c,x)
y=c(k);
for i=k:-1:2
    y=y*(xi-x(i-1))+c(i-1);
end
