t1.m

1． 对薄膜形成与生长中相关物理过程及现象进行分析和建模。 2． 利用蒙特卡罗（Monte Carlo）方法和分子动力学（Molecular Dynamics）方法对薄膜形成与生长过程进行计算机模拟

%t1.m
%判断是否为黑
if(A(m,n)==2)%A(m,n)==2表示黑色；
fill([n-1,n,n,n-1,n-1],[m-1,m-1,m,m,m-1],'k');%若为红，就涂黑；
if(m~=M&A(m+1,n)==0)%若为黑，则上面格子涂红，以防黑横着走；
    fill([n-1,n,n,n-1,n-1],[m,m,m+1,m+1,m],'r');
elseif(m~=M&n~=1&A(m,n-1)==0)%若为黑，则左面格子涂红；
     fill([n-2,n-1,n-1,n-2,n-2],[m-1,m-1,m,m,m-1],'r')
elseif(n~=M&A(m,n+1)==0)%若为黑，则右面格子涂红；
    fill([n,n+1,n+1,n,n],[m-1,m-1,m,m,m-1],'r')
elseif(m~=1&A(m-1,n)==0)%若为黑，则下面格子涂红；
    fill([n-1,n,n,n-1,n-1],[m-2,m-2,m-1,m-1,m-2],'r')
else
    m=m;n=n;
end
end
%判断是否为红
if(A(m,n)==1)%A(m,n)==1表示红色；
fill([n-1,n,n,n-1,n-1],[m-1,m-1,m,m,m-1],'k');%若为红，表示空穴以进入活性区，下时刻一定变黑；
if(m~=M&A(m+1,n)==0)%红上边下时刻将变为红（前提为白）；
    fill([n-1,n,n,n-1,n-1],[m,m,m+1,m+1,m],'r');
elseif(n~=1&A(m,n-1)==0)%红左边下时刻将变为红；
     fill([n-2,n-1,n-1,n-2,n-2],[m-1,m-1,m,m,m-1],'r')
elseif(n~=M&A(m,n+1)==0)%红右边下时刻将变为红；
    fill([n,n+1,n+1,n,n],[m-1,m-1,m,m,m-1],'r')
elseif(m~=1&A(m-1,n)==0)%红下边下时刻将变为红（前提为白）；
    fill([n-1,n,n,n-1,n-1],[m-2,m-2,m-1,m-1,m-2],'r')
else
    m=m;n=n;
end
end
%判断是否邻红
if(A(m,n)==0)%A(m,n)==0表示白色；
    up;
    down;
    left;
    right;
    rightup;
    rightdown;
    leftup;
    leftdown;
end