function BP = Bernpol(f,n)
% BP = Bernpol(f,n)
% Bernsteinuv polynom stupne n funkce f
% f - inline funkce
if (n<1) || (round(n)~=n)
    error('Nespravny stupen polynomu')
end
BP=zeros(1,n+1);
koreny=ones(1,n);
znam=(-1)^n;
for k=0:n
    bk=nchoosek(n,k)*poly(koreny)*znam;
    BP=BP+f(k/n)*bk;
    znam=-znam;
    if k<n
        koreny(k+1)=0;
    end
end

end

