function[output] = simulace_jablka(n, k)
% n = pocet pokusu
% k = pocet jablek v kosiku
vysledky = 1:n;
% vytvorim vektor s n prkvy kam budu ukladat vysledky
for i = 1:n
    vysledky(i) = jablka(k);
    % na i-tou pozici vektoru 'vysledky' ulozim vysledek pokusu
end

subplot(1,2,1)
    histogram(vysledky, 'facecolor', 'r', 'binwidth', 1)
    % sirku sloupce nastavim na 1 abychom videli absolutni pocet pokusu ktere dopadly stejne
    title('vysledky jednotlivych pokusu')
    xlabel('pocty deti')
    ylabel('cetnosti')
    
subplot(1,2,2)
    pocet_jablek = k;
    % do tohoto vektoru budu ukladat klesajici pocet jablek v prubehujednoho pokusu
    % prvnim elementem bude pocet jablek, se kterym zaciname, tedy k
    snezena_jablka = 0;
    % vytvorim vektor do ktereho budu zaznamenavat snezena jablka
    i = 1;
    % promenna i bude znacit, na jakou pozici vektoru 'pocet_jablek' chci ukladat novou hodnotu
    while snezena_jablka <= k
        identifikator = randi([0, 1], 1);
        % zda dite snedlo jablko, nebo nesnedlo
        pocet_jablek(i + 1) = pocet_jablek(i) - identifikator;
        % do vektoru 'pocet_jablek' ulozime novy stav poctu jablek
        snezena_jablka = snezena_jablka + identifikator;
        % zaznamename, kolik jablek bylo snezeno
        i = i + 1;
        % posuneme se na dalsi pozici
    stairs(pocet_jablek, 'r')
    ylim([0, k + 1])
    title('trajektorie jednoho pokusu')
    xlabel('pocet deti')
    ylabel('pocet jablek')
    end
end