%% Video at: http://personales.upv.es/asala/YT/V/wienermlab1.html
Ts=0.0001;
NPow=17;
T=(0:(2^17))*Ts;

A=1;B=sqrt(Ts);
sys=ss(A,B,1,0,Ts);
u=randn(length(T),1);
y=lsim(sys,u);
limxy=[min(T) max(T) min(y) max(y)];
T2=T;y2=y;
%plot(T,y) TODO.
for k=1:NPow
    Tstep=2^(NPow-k);
    elcolor=[0 0.8 0.9].^k;
    if(Tstep>=1)
        plot(T2(1:Tstep:end),y2(1:Tstep:end),'.',MarkerSize=2+NPow*2-2*k,Color=elcolor)%rand(1,3))
    else 
       plot(T,y,Color=elcolor)%rand(1,3))
    end
    T2(1:Tstep:end)=0;
    y2(1:Tstep:end)=0;
    grid on
    axis(limxy)
    hold on
    drawnow
    pause
end   
hold off
return 
%% 

%versión 2... convergencia de sumas... cada vez más "apretadas"...
pause
Nmax=16000;
T=0:(Nmax-1);
Nrealiz=60;
y=zeros(Nrealiz,Nmax);
oldn=0;
for n=1:(Nmax-1)
    y(:,n+1)=y(:,n)+randn(Nrealiz,1); %veremos varias realizaciones
    Tesc=T(1:(n+1))/n;
    yesc=y(:,1:(n+1))/sqrt(n);
    %no dibujo "siempre"
    if(n<=250) || (n>=1.005*oldn)
        plot(Tesc,yesc',LineWidth=2), grid on
        axis([0 1 -3 3])
        drawnow
        oldn=n;
    end
    %if(~mod(n,50))
    %    n
    %end
    %pause
end
