function [sr_vred,std_dev] = crtaj_raspodelu (gaus)
subplot(2,1,1)
plot(gaus)
[g,gbin]=hist(gaus,100);
dxg=gbin(2)-gbin(1);
n=length(gaus);
fg=g/n/dxg;
subplot(2,1,2)
bar(gbin,fg)
sr_vred=mean(gaus)
std_dev=std(gaus)
endfunction