% Lagrange-Interpolation von sin auf dem Intervall [0,2*pi]

a = 0; b = 2*pi;
xmany = linspace(a,b,500);
for n = 2:2:20,
   x = linspace(a,b,n+1);
   y = sin(x);
   ymany = sin(xmany);
   lambda = baryzentrisch(x);
   p = baryzentrischauswertung(x,xmany,y,lambda);
   subplot(2,1,1);
   plot(xmany,ymany,'b-',xmany,p,'r-');
   makelarge
   title(sprintf('Interpolation, Polynomgrad %i',n));
   subplot(2,1,2);
   semilogy(xmany,abs(ymany-p),'r-');
   makelarge
   title('Interpolationsfehler');
   pause
end