% tehtävä 5
function y = teht5(rounds)
xs = linspace(-2,2);
xsnorm = linspace(0,1);


plot(xs, f(xs), 'black')
hold on;

for n = 1:rounds
    % Bernsteiniä
    B = zeros(1,n+1);
    for k = 1:n
        B = B + fphi(k/n) * bernstein(n,k,n);
    end
    plot(xs, polyval(B, xsnorm), 'r');
    drawnow
    pause(0.2)
    % jos vielä saisi pyyhittyä edelliset pois...
end


hold off;


y = 0;

% tehtävän 5 funktio f
function y = f(x)
y = cos(x) + sqrt(abs(x));

% f . phi
function y = fphi(x)
y = f(phi(x, -2, 2));