%move.m
h=rand(1);
if (0<=h&h<0.25)%空穴下一時刻向下移動；
    if(m~=1)
m=m-1;    
    end
end
if(0.25<=h&h<0.50)%空穴下一時刻向上移動；
    if(m~=M)
m=m+1;
    end
end
if(0.50<=h&h<0.75)%空穴下一時刻向左移動；
    if(n~=1)
n=n-1;
    else
        n=1;
    end
end
if(0.75<=h&h<1)%空穴下一時刻向右移動；
    if(n~=M)
n=n+1;    
    end
end
B=B+1;%空穴擴散一步，平均自由程減一；