//[]
clf;
N=1000;
U=[];

for k=1:N,
  u=4*rand()-2;
  v=0.4*rand();
    while v>sqrt(4-u^2)/(2*%pi), 
    u=4*rand()-2;
    v=0.4*rand();
    end
  U=[U,u];
end
arithm=1:N;
S=cumsum(U)./arithm;
aux=ones(1,N);
suitmin=-(1.96*aux./sqrt(arithm));
suitmax=(1.96*aux./sqrt(arithm));
x=[0,N]; y=[0,0];
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()