Nt=100;
Np=30;
RN=2*randi(2,Np,Nt)-3;
RW=cumsum(RN,2);
meanpos=mean(RW,1);
stdpos=std(RW,'',1);
plot(stdpos.^2,'o')

