%========================================================================
%?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=