利用ftont-tracking 方法在matlab中模拟水滴下落

%========================================================================
% code2.m:
% A very simple Navier-Stokes solver for a drop falling in a rectangular
% box. A forward in time centered in space discretization is used.
% The density is advected by a front tracking scheme and a stretched grid
% is used allowing us to concentrate the grid points in specific areas
%========================================================================
%domain size and physical variables
Lx=1.0;Ly=1.0;gx=0.0;gy=-100.0; rho1=1.0; rho2=2; m0=0.01; rro=rho1;
unorth=0;usouth=0;veast=0;vwest=0;time=0.0;
rad=0.15;xc=0.5;yc=0.7; % Initial drop size and location
% Numerical variables
nx=32;ny=32;dt=0.00125;nstep=200; maxit=200;maxError=0.001;beta=1.2;
% Zero various arrys
u=zeros（nx+1ny+2）; v=