move.m

来自「1． 对薄膜形成与生长中相关物理过程及现象进行分析和建模。 2． 利用蒙特卡罗」· M 代码 · 共 27 行

%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;%空穴扩散一步，平均自由程减一；