function p=newton(x,y)
T=Tabladifer(x,y)
n=length(x);
p=0;
for k=1:n
    p=sumapoly(p,T(k,k)*poly(x(1:k-1)));
end