function [c]=newton(x,y)
ff = y;
Matrix=ff;
m=length(x);
for j=2:m
    for i=1:(m-j+1)
        ff(i)=( ff(i+1) - ff(i) ) / ( x(j+i-1) - x(i));
    end
    Matrix=[Matrix ff];
end
for i=2:m
    for j=0:i-2
        Matrix(i,m-j)=0;
    end
end 
figure
uitable('Data',Matrix)
c=Matrix(1,:);
end