function interpolacion = interpnewton(puntos,valores)
% **********************************************************************
% ********************** AUTOR: José Miguel Gómez **********************
%
% Esta función calcula el polinomio interpolador de una función usando
% el método de Netown.
% 
% Las variables de entrada son:
% puntos: Puntos usados para interpolar
% valores: Valores de la función en dichos puntos.
%
% **********************************************************************
% **********************************************************************

n=size(puntos,2);
m=size(puntos,2);
if n~=m
   disp('ERROR:los vectores de datos tienen que ser de la misma dimensión');
else
   aux=valores;
	interpolacion=[valores(1)];
	A=[1];
	for i=2:n
   	for j=i:n
      	aux(j)=(valores(j)-valores(j-1))/(puntos(j)-puntos(j-i+1));
   	end;
   	valores([i:n])=aux([i:n]);
   	A=conv(A,[1,-puntos(i-1)]);
   	interpolacion=polisuma(interpolacion,valores(i)*A);
	end;
end;