% generazione di variabili gaussiane
% applicazione dell'algoritmo di Box-Muller

clear;
clf;
figure(1);
clf;
n=10000;
x=sqrt(-2*log(rand(n/2,1)));
y=pi*(2*rand(n/2,1)-1);
z=x.*cos(y);
w=x.*sin(y);
zw=[z;w];
nbin=50;
b=5;
dx=2*b/nbin;
xgr=-b-dx/2:dx:b+dx/2;
nc=histc(zw,xgr)/(n*dx);
bar(xgr,nc);
y=(1/sqrt(2*pi))*exp(-((xgr+dx/2).^2)/2);
hold on
plot(xgr,y,'-r');
hold off







