function [aprox_pi, sigma] = ukol2(itermax, delta)
% [aprox_pi, sigma] = ukol2(itermax)
%
% funkce vypocte metodou Monte Carlo s body generovanymi pomoci nahodne prochazky obsah plochy ohranicene grafem funkce sqrt(1 - x^2), 
% kladnym smerem osy x a osou y a smerodatnou odchylku sigma 
% pomoci zmineneho obsahu aproximuje cislo pi
%
% itermax ... pocet pokusu
% delta ... maximalni delka hodu

x = rand(); % generuje x-ovou souradnici bodu ve ctverci <0,1> x <0,1>
y = rand(); % generuje y-ovou souradnici bodu ve ctverci <0,1> x <0,1>

q = zeros(1,itermax);

% doplnte kod
%
%
%
%
%
% doplnte kod

% vypis vysledku
disp(['Aproximace pi: ', num2str(aprox_pi, '%.15f')])
%disp(['Absolutni chyba: ', num2str()])
%disp(['Relativni chyba: ', num2str(),'%'])

mu = sum(q) / double(itermax);

sigma2 = 1/itermax * (q - mu) * (q - mu)';
sigma = 2 * sqrt(sigma2)/sqrt(itermax);

% konec
