function [n,DD]=Newton_int_1(x,y)
% Inperpolación con polinomios de Newton
%  vectores de entrada x e y; 
%  salida: coeficientes en n y diferencias en DD 
% 
 N  = length(x)-1;
 DD = zeros(N+1,N+1); DD(1:N+1,1) = y(:);
 for k=2:N+1 
    for m=1:N+2-k
       DD(m,k)=(DD(m+1,k-1)-DD(m,k-1))/(x(m+k-1)-x(m)); % Diferencias divididas
    end
 end
 a = DD(1,:); n = a(N+1); 
 for k=N:-1:1
    n = [n a(k)] - [0 n*x(k)];                          % n(x)*(x-x(k-1))+a_k-1
 end
end
