clf;
N=1000;
X=rand(1,N);
arithm=1:N;
S=cumsum(X)./arithm;
aux=ones(1,N);
suitmin=0.5*aux-(1.96*aux./sqrt(12*arithm));
suitmax=0.5*aux+(1.96*aux./sqrt(12*arithm));
x=[1,N]; y=[0.5,0.5];
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()
