function p = difdiv(x,y)
% Esta funcion calcula el polinomio interpolador mediante diferencias divididas
% y posterior simplificacion de la expresion del polinomio interpolador de Newton.
	n=length(y);
	d=y;
	p=d(1);
	for k=1:n-1
		for i=1:n-k
			d(i)=(d(i)-d(i+1))/(x(i)-x(i+k));
		end
		d(length(d))=[ ];
		p=[0,p]+d(1)*poly(x(1:k));
	end
end