function p = newton(x,y)

n = length(x);
p=y;

for i=2:n
  for j=n:-1:i
    p(j) = (p(j)-p(j-1))/( x(j)-x(j-i+1) );
  end
end

return