function [p,er]=diffdiv(x,y)
p=[];
er='';
nx=length(x);
ny=length(y);
%% Verifica dei dati in ingresso				
if ny ~=nx
  er='x e y devono avere lo steso numero di elementi';
  return;
end
%% Costruzione della matrice delle differenze divise
A=y;
for k=2:nx
  A(k:nx)=(A(k:nx)-A(k-1:nx-1))./(x(k:nx)-x(1:nx-k+1));
end
%% Costruzione del polinomio interpolatore
p=A(1);
q=[1];
for k=2:nx
  p=[0,p];
  q=conv(q,[1,-x(k-1)]);
  p=p+A(k)*q;
end