//[]
clf;
N=1000;
U=sqrt(rand(1,N));


arithm=1:N;
S=cumsum(U)./arithm;
aux=ones(1,N);
suitmin=(2/3)*aux-(1.96*aux./(3*sqrt(2*arithm)));
suitmax=(2/3)*aux+(1.96*aux./(3*sqrt(2*arithm)));
x=[0,N]; y=[2/3,2/3];
plot2d(arithm,S)
plot2d(arithm,suitmin,4)
plot2d(arithm,suitmax,6)
plot2d(x,y,5)
legend("S_n","borne inferieure du TCL a 95%","borne superieure du TCL a 95%","y=m")
xtitle("Visualisation du TCL","n")
xselect()