function v = interpolante(x,y,u)
% calcula la interpolante polinomial por puntos (x,y),
% y la evalua en u
  n=length(x);
  v=zeros(size(u));
  for k=1:n
     w=ones(size(u));
     for j=[1:k-1 k+1:n]
        w=(u-x(j))./(x(k)-x(j)).*w;
     endfor
     v=v+w*y(k);
  endfor
endfunction