元胞自动机 Matlab实现表面张力

clf
clear all
%build the GUI
%define the plot button
plotbutton=uicontrol('style','pushbutton',...
    'string','Run',...
    'fontsize',12,...
    'position',[100,400,50,20],...
    'callback','run=1;');


%define the stop button
erasebutton=uicontrol('style','pushbutton',...
    'string','Stop',...
    'fontsize',12,...
    'position',[200,400,50,20],...
    'callback','freeze=1;');


%define the Quit button
quitbutton=uicontrol('style','pushbutton',...
    'string','Quit',...
    'fontsize',12,...
    'position',[300,400,50,20],...
    'callback','stop=1;close;');


number=uicontrol('style','text',...
    'string','1',...
    'fontsize',12,...
    'position',[20,400,50,20]);
%=================================================
%CA setup
n=128;
%initialize the arrays
z=zeros(n,n);
cells=z;
sum=z;
%set a few cells to one
cells(n/2,.25*n:.75*n)=1;
cells(.25*n:.75*n,n/2)=1;
cells=(rand(n,n))<.5;
%how long for each case to stability or simple oscillato