%%
%% TP1 Exercice 4
%% Simulation de la loi normale par
%% l'algorithme de Box-Muller
%%

N=input('Entrer la taille de l''échantillon N : ');
m=input('Préciser la valeur de la moyenne m : ');   
sigma2=input('Préciser la valeur de la variance : ');   
X=m*ones(N,1)+sqrt(sigma2)*sqrt(-2*log(rand(N,1))).*cos(2*pi*rand(N,1));
[E,C]=histo(X,sqrt(N),0,1,'b');
hold on
title('Simulation de la loi normale');
xlabel('Valeurs');
ylabel('Effectifs');
plot(C,dnorm(C,m,sqrt(sigma)),'r--');
legend('Empirique','Theorique');
hold off