function [aprox_pi, sigma] = ukol3_reseni(itermax)
% [aprox_pi, sigma] = ukol3_reseni(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);

for i = 1:itermax
    x = rand(); % generuje x-ovou souradnici bodu v intervalu <0,1>
    s(i) = sqrt(1.0 - x^2); % funkcni hodnota funkce sqrt(1 - x^2)
end

aprox_pi = 4.0 * sum(s) / double(itermax);


% vypis vysledku
disp(['Aproximace pi: ', num2str(aprox_pi, '%.15f')])
disp(['Absolutni chyba: ', num2str(abs(aprox_pi - pi), '%.15f')])
disp(['Relativni chyba: ', num2str(100*abs(aprox_pi - pi)/pi, '%.3f'),'%'])

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

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

% konec

