///////////////////////////////////////////////////////////////////////////////////
//                  ECHANTILLONNAGE D'IMPORTANCE          (EXO2)             //////
///////////////////////////////////////////////////////////////////////////////////


stacksize(100000000);

//Nombre de simulation (à chaque expérience)
n=10000;

//Nombre de répétitions de l'expérience
nb=10;

//Simulation des gaussiennes centrées réduites
x=grand(n,nb,'nor',0,1);

//g-Valeur de la sortie g(Y)=g(f(X)), avec f(x)=x et g(x)=(x-2)_+
y=max((x-2),0);

//Les estimateurs empiriques de E(g(Y))=E((Y-2)_+)
my=cumsum(y,'r')./((1:n)'*ones(1,nb));


if nb<51 then

clf(0);xset("window",0);show_window();
plot2d(my);
xtitle(['Nombre d''expériences : "+string(nb)+" , Nombre de simulations par expérience : "+string(n)+" ',
'estimation empirique de E((Y-2)_+) pour Y loi N(0,1)']);



end