Случайное блуждание

n = 100;
x = n/2; y = n/2;
plot(x,y,'.','EraseMode','none')
axis([0,n,0,n]); axis square
hold on

for k = 1:300
    X = [x-1,x,x+1];
    Y = [y-1,y,y+1];
    x = X( round(1+2*rand) ); y = Y( round(1+2*rand) );
    plot(x,y,'.','EraseMode','none'); %     pause(.2)
    M(k) = getframe;
end

movie2avi(M,'proba.avi')