function obsah = montecarlo(n)
% obsah = montecarlo(n)
%
% vypocte obsah plochy ohranicene shora grafem funkce x^2 na intervalu <0,1>
%
% input:
% n ... pocet pokusu

if n <= 0 | n - floor(n) > eps
    error('n neni prirozene cislo')
end

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, osou x a primkami x=0, x=1
    end
end

obsah = suma/n;

disp(['Obsah plochy pod grafem funkce x^2 na intervalu <0,1>: ', num2str(obsah)])

% konec


