function y = intpolyval (x, xi)
%
% function y = intpolyval (x, xi)
%
% Valuta il polinomio interpolatore delle coppie [xi(:,1),xi(:,2)] sui
% punti x, mediante polinomi elementari di Lagrange.
n = size (xi, 1) - 1;
y = zeros (size (x));
for j = 1:n + 1
  y = y + xi(j, 2) * Li (j, x, xi(:,1));
end
%!demo
%! xi = rand (5, 2);
%! x = linspace (0, 1);
%! y = intpolyval (x, xi);
%! plot (x, y, xi(:, 1), xi(:, 2), 'o')
%!test
%! xi(:, 1) = linspace (0, 1, 5);
%! xi(:, 2) = sin (xi(:, 1));
%! x = linspace (0, 1);
%! y1 = intpolyval (x, xi);
%! V = vander (xi(:, 1));
%! a = V \ xi(:, 2);
%! y2 = polyval (a, x);
%! assert (y1, y2, 1e-12)
