function obsah = montecarlo_reseni(n)
% obsah = montecarlo_reseni(n)
%
% funkce vypocte obsah plochy ohranicene grafem funkce y = x^2, osou x a primkou x = 1 
%
% n ... pocet pokusu

suma = 0;

for i = 1 : n
    bod = rand([1,2]); % generuje bod ve ctverci <0,1> x <0,1>
    x = bod(1); % x-ova souradnice bodu
    y = bod(2); % y-ova souradnice bodu
    
    if y <= x^2
        suma = suma + 1; % bod lezi v plose ohranicene grafem funkce y = x^2, osou x a primkou x = 1
    end
end

obsah = suma / n;

% vypis vysledku
disp(['Obsah plochy ohranicene grafem funkce y = x^2, osou x a primkou x = 1: ', num2str(obsah, '%.8f')])

% konec


