function Lnf = lagrange_barycentric(nodes,f_nodes,x)
  n = length(nodes);
  m = length(x);
  Lnf = zeros(size(x));
  l = zeros(size(nodes));
  for i = 1:m
    ux = prod( x(i) - nodes([1:n]));
    for j = 1:n
      wi(j) = 1/prod( nodes(j) - nodes([1:j-1,j+1:n]));
    endfor
    Lnf(i) = ux* ((wi./(x(i) - nodes([1:n])))*f_nodes');

  endfor
end;