function [realizace]=BM_transformace(mi,sigma,n)
% funkce generuje cisla z normalniho rozlozeni pomoci BM transformace
% syntaxe: realizace=BM_transformace(mi,sigma, n)
% vystupni parametr:
% realizace ... vektor realizaci nahodne veliciny s rozlozenim N(mi, sigma^2)
% vstupni parametry:
% mi ... stredni hodnota
% sigma ... smerodatna odchylka
% n ... pocet realizaci (musi byt sude cislo)
realizace=[];
for i=1:n/2
u1=unifrnd(0,1,1,1);
u2=unifrnd(0,1,1,1);
x1=sqrt(-2*log(u1))*cos(2*pi*u2);
x2=sqrt(-2*log(u1))*sin(2*pi*u2);
realizace=[realizace;x1;x2];
end
sh=mi*ones(n,1);
realizace=sh+sigma*realizace;
hist(realizace)
