function [aprox_pi, sigma] = ukol3(itermax)
% [aprox_pi, sigma] = ukol3(itermax)
%
% funkce vypocte metodou Crude Monte Carlo 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

s = 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(s) / double(itermax);

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

% konec

