function kvalitaInterpolace()

a = -1; b = 1;
xx = a:0.01:b;

for n=1:20
    xLagr = a:((b-a)/n):b;
    kvalitaLagr = xx; kvalitaLagr(:) = 1;
    i = 0:n;
    xCeb = (b+a - 2*cos((2*i+1)*pi/2/(n+1))) ./ (b-a);
    kvalitaCeb = xx; kvalitaCeb(:) = 1;
    for i=0:n
        kvalitaLagr = kvalitaLagr .* (xx-xLagr(i+1));
        kvalitaCeb = kvalitaCeb .* (xx-xCeb(i+1));
    end
    plot(xx,kvalitaLagr,'b',xx,kvalitaCeb,'r')
    [max(abs(kvalitaLagr)),max(abs(kvalitaCeb))]
    pause
end