%原子團遷移
DD=zeros(222);
EE=zeros(222);
EE(:)=100;
r=ceil((100-E)/10); %搜索半徑
for i=12:211
    for j=12:211
        a=r(i,j);
        if a<2    %若半徑過小則不執行
            continue
        end
        for k=2:a
            for ii=i-k:i+k
                for jj=j-k:j+k
                    if ii==i&jj==j   %判斷自身點的時候跳過執行程序
                        continue
                    else if E(ii,jj)<100
                            m=fix(i+((ii-i)*N(ii,jj))/(N(i,j)+N(ii,jj)));
                            n=fix(j+((jj-j)*N(ii,jj))/(N(i,j)+N(ii,jj)));%m,n為遷移后的坐標
                            DD(m,n)=D(i,j)+D(ii,jj);%遷移后的粒子個數
                            EE(m,n)=E(i,j)-D(ii,jj);%遷移后的能量
                        end
                    end
                end
            end
            break
        end
    end
end