LKanalikow = 60;
pozycja = [];
pozycja1 = [];
przesun = LKanalikow/2+1;
Kanaliki = zeros(1,2*przesun);
figure(99);clf
xvec = -przesun+1:przesun;
hbar = bar(xvec,Kanaliki);hold on
kulki = plot(0,LKanalikow,'ko');
set(kulki,'markerfacecolor','blue')
axis([-przesun, przesun, 0, 250]);
for i = 1:1000
    pozycja = [0,LKanalikow*0.3; pozycja]
    % wykres
    set(hbar,'ydata',Kanaliki);
    set(kulki,'xdata',pozycja(:,1),'ydata',pozycja(:,2)*10);

    drawnow

    % add constant -y and random x
    [m,n]=size(pozycja);
    stepx = sign(randn(m,1))/2 -1.2;
    stepy = -1*ones(m,1);
    pozycja = pozycja + [stepx, stepy];

%     check bottom of stack for zero y crossing...
%     if so, remove and record x position (an integer)
    if pozycja(end,2) <= 0
        xfinal = pozycja(end,1);
        kanalik = round(xfinal+przesun);
        Kanaliki(kanalik) = Kanaliki(kanalik) + 1;
        pozycja(end,:) = [];
    end
    pozycja1 = [0,LKanalikow*0.3; pozycja1]
    % wykres
    set(hbar,'ydata',Kanaliki);
    set(kulki,'xdata',pozycja1(:,1),'ydata',pozycja1(:,2)*10);

    drawnow

    % add constant -y and random x
    [m,n]=size(pozycja1);
    stepx1 = sign(randn(m,1))/2 + 0.5;
    stepy1 = -1*ones(m,1);
    pozycja1 = pozycja1 + [stepx1, stepy1];

%     check bottom of stack for zero y crossing...
%     if so, remove and record x position (an integer)
    if pozycja1(end,2) <= 0
        xfinal = pozycja1(end,1);
        kanalik = round(xfinal+przesun);
        Kanaliki(kanalik) = Kanaliki(kanalik) + 1;
        pozycja1(end,:) = [];
    end
end