function films(aa,bb,cc)  %aa打粒子的次數(shù) bb一次打上粒子的個數(shù) cc整個循環(huán)的次數(shù)
Z=zeros(200);
Z(:)=100;
for ccc=1:cc
    
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
ad=adsorption_reflect(A,N,M)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
vibrate
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
cg_energy.m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
move.m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fanning.m
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    
DDD=100-DD;
M=DDD(12:211,12:211);
M=Z+M-100;
imshow(M,[min(min(M)),max(max(M))]);
Z=M;      %運行完一次后，把M記下，以便再次顯示時能方便疊加顯示
end