function d = diffdiv1(xi)
%
% function d = diffdiv1 (xi)
%
% Calcolo delle differenze divise dei punti [xi(:,1),xi(:,2)]. I nodi
% xi(:,1) non devono essere ripetuti.
yi = xi(:, 2);
xi = xi(:, 1);
n = length (xi) - 1;
d = yi;
for i = 1:n + 1
  for j = 1:i - 1
    d(i) = (d(i) - d(j)) / (xi(i) - xi(j));
  end
end
