% Diese Funktion wertet das Lagrange'sche Interpolationspolynom aus
% in der baryzentrischen Darstellung aus.

function y = barycentric (x, Lambda, X, Y)

  U = (x-X);
  [u_min, i] = min (abs(U));
  
  % Falls x sehr nahe an einer Stuetzstelle liegt, berechnen wir nichts,
  % sondern geben den dazugehoerigen y-Wert aus.

  if u_min < 100*eps
    y = Y(i);
  else
    W = Lambda ./ U;
    y = sum (Y .* W) / sum (W);
  endif