%
%
%       Biofisica 2006 - FCEyN - UBA
%       Random Walks
%       Jacobo Sitt
%

function peli(n)

ps=200; % pasos


for i=1:n
  
w(i,:)=rw(.9,ps);   % hace una tabla con n random walk 
                     % la fila i tiene el i-esimo RW
end

for i=1:ps
    m(i)=mean(w(:,i)); % calcula la media para la columna i ( paso i de n RW's)
    v(i)=var(w(:,i));  % calcula la d.est para la columna i ( paso i de n RW's)
end

x=-ps:.2:ps; % arma un vector para todos los posibles x


 M = moviein(n);  % inicializa la matriz de imagenes
 
       for j=1:ps
        
        y = 1./sqrt(2*3.14*v(j)).*exp(-(x-m(j)).^2./(2*v(j)))*n; % calcula la 
            % distrucion gaussiana para todos los x con la media y la var
        hist(w(:,j),100) % plotea el histograma de los rw
        hold on
        plot(x,y) % plotea la distribucion
                hold off
       
        axis([-ps ps 0 n])
        

        M(:,j) = getframe; %captura la imagen para la peli
        
       end
       movie(M,10) % muestra la peli 10 veces