function [realizace]=clv_polynom(mi,sigma,n)
% funkce generuje cisla z normalniho rozlozeni pomoci CLV
% pouziva aproximaci polynomem
% syntaxe: realizace=clv_polynom(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
a1 = 3.949846138;
a3 = 0.252408784;
a5 = 0.076542912;
a7 = 0.008355968;
a9 = 0.029899776;
realizace=[];
for i=1:n
u=unifrnd(0,1,12,1);
x=sum(u)-6;
y=0.25*x;
u=a1*y+a3*y^3+a5*y^5+a7*y^7+a9*y^9;
realizace=[realizace;u];
end
sh=mi*ones(n,1);
realizace=sh+sigma*realizace;
hist(realizace)
