% plot the localization properties of the Bernstein basis polynomials
nn=[6,20,50]; % degrees to choose from

x=linspace(0,1,100);
[xx,yy]=meshgrid(x,x);
for k=1:length(nn)
    figure
    n=nn(k);
    z=bernstein(n,x);
    [xx1,yy1]=meshgrid(x,[0:n]);
    mesh(xx1,yy1,z)
    title 'Bernstein Grundpolynome (durchgezogene Linien)'
    xlabel 'x'
    ylabel '0\leq k\leq n'

    figure
    zz=(n+1)*(z'*z);  % kernel of the Bernstein Durrmeyer operator
    mesh(xx,yy,zz)
    title 'Kern des Bernstein-Durrmeyer Integraloperators'
    xlabel 'x'
    ylabel 'y'
end