% Zuerst mal die Funktion f zum Plotten auswerten

XX = -5:0.05:5;
YY = f(XX);
PP = zeros(1,length(XX));
TT = zeros(1,length(XX));

% Dann Interpolation mit aequidistante Stuetzstellen

X = -5:5;
Y = f(X);
L  = lambda(X);

for i=1:length(XX)
  PP(i) = barycentric (XX(i), L, X, Y);
end

% Jetzt Tschebycheff Stuetzstellen

K = 0:10;
T = 5 * cos (pi*(10-K)/10);
Z = f(T);
L  = lambda(T);

for i=1:length(XX)
  TT(i) = barycentric (XX(i), L, T, Z);
end

% Und schliesslich alles plotten

plot (XX, YY, XX, PP, XX, TT, X, Y, "*", T, Z, "o")
grid, replot