function ymany = baryzentrischauswertung(x,xmany,y,lambda)
% Wertet Interpolationspolynom in baryzentrischer Darstellung aus.
% Basierend auf einem Vorschlag von [Berrut/Trefethen, SIAM Review'04].

n = length(x)-1;
zaehler = zeros(size(xmany));
nenner  = zeros(size(xmany));
exakt   = zeros(size(xmany));
for j = 1:n+1
   xdiff = xmany-x(j);
   temp = lambda(j)./xdiff;
   zaehler = zaehler + temp*y(j);
   nenner = nenner + temp;
   exakt(xdiff==0) = j;
end
ymany = zaehler./nenner;
jj = find(exakt); ymany(jj) = y(exakt(jj));