newto divididas

function newton(x,b,xk)

n=length(xk);

for i=1:n

   f(i,1)=b(i);

end

for j=2:n

   k=j-1;

   for i=1:n+1-j

      f(i,j)=(f(i+1,k)-f(i,k))/(a(i+j-1)-a(i));

   end

end

prod=1;

sum=0;

for j=1:n

   sum=sum+f(1,j)*prod;

   prod=prod*(xk-x(j));

end

fprintf('\n        x           F[]\n');

fprintf('=============');

for j=1:n

   fprintf('=============');

end

fprintf('\n');

for i=1:n

   fprintf('%10.6f',x(i));

   for j=1:n+1-i

      fprintf('%10.6f',f(i,j));

   end

end

P=sum;

disp('EL RESULTADO DE LA INTERPOLACION ES:');

disp(P);