for real=1:10

iter=3000;

N=100;
rho=zeros(1,N);

k=50;
store=[k];
for i=1:iter
    
    p=rand(1);
    if p>0.5
        if k<N
        k=k+1;
        else
        k=1;
        %k=k;
        end
    else
        if k>1
        k=k-1;
        else
        k=N;
        %k=k;
        end
    end
    store=[store k];
end

plot(store,'-o')
hold on
axis([0 iter 0 100])
end