function IteracniPolynomy()

A = mkd1d('sin',0.1);
Lambda = eig(full(A));
l1 = min(Lambda);
ln = max(Lambda);

l=l1:0.1:ln;
for k=1:10,
    i=1:k,
    xi=cos(pi*(i-1/2)./k);
    PR = (1-(2/(l1+ln))*l).^k; % Rich.
    chi=cos(k*acos((2*l-l1-ln)/(ln-l1)))./cos(k*acos((2*0-l1-ln)/(ln-l1))); % Ceb.
    plot(l,chi,'r',l,PR,'b');
    pause
end