function Int_New(xi,x,y)

n = columns(x);
M(n,n)=0;

for i = 1:n
   M(i,1) = y(i);
endfor

for i = 1:n
   for j = 1:n-i
       M(j,i+1) = (M(j+1,i)-M(j,i))/(x(j+i)-x(j));
   endfor
endfor

p  = 1;
yi = M(1,1);
for i = 1:n-1
   p  = p*(xi-x(i));
   yi = yi + p*M(1,i+1);
endfor
M
yi
end