% FEMLAB Model M-file
% Generated 06-Jun-2002 10:40:36 by FEMLAB 2.1.0.127.

flclear fem
% FEMLAB Version
clear vrsn;
vrsn.name='FEMLAB 2.1';
vrsn.major=0;
vrsn.build=127;
fem.version=vrsn;

% Space dimensions
fem.sdim={'x','y'};

% Geometry
clear s c p
R1=rect2(-0.80000000000000004,0.80000000000000004,-0.80000000000000004, ...
0.80000000000000004,0);
objs={R1};
names={'R1'};
s.objs=objs;
s.name=names;

objs={};
names={};
c.objs=objs;
c.name=names;

objs={};
names={};
p.objs=objs;
p.name=names;

drawstruct=struct('s',s,'c',c,'p',p);
fem.draw=drawstruct;
fem.geom=geomcsg(fem);

% Initialize mesh
fem.mesh=meshinit(fem,...
	'Out',    {'mesh'},...
	'jiggle', 'mean',...
	'Hcurve', 0.29999999999999999,...
	'Hgrad',  1.3,...
	'Hpnt',   {10,[]});

% Refine mesh
fem.mesh=meshrefine(fem,...
	'out',    {'mesh'},...
	'rmethod','regular');

% Refine mesh
fem.mesh=meshrefine(fem,...
	'out',    {'mesh'},...
	'rmethod','regular');

% Dimension
fem.dim={{'u'}};

% Boundary conditions
clear bnd
bnd.q={{{'0'}}};
bnd.g={{{'0'}}};
bnd.h={{{'0'}}};
bnd.r={{{'0'}}};
bnd.var={};
bnd.ind=ones(1,4);
fem.bnd=bnd;

% PDE coefficients
clear equ
equ.da={{{'1./(ux.^2+uy.^2).^0.5'}}};
equ.c={{{'1./(ux.^2+uy.^2).^0.5'}}};
equ.al={{{'0';'0'}}};
equ.ga={{{'0';'0'}}};
equ.be={{{'0';'0'}}};
equ.a={{{'0'}}};
equ.f={{{'0'}}};
equ.var={'absux','sqrt(ux.^2+uy.^2)','abscu1x','sqrt(cu1x.^2+cu1y.^2)'};
equ.ind=1;
fem.equ=equ;

% Evaluate initial condition
fem.init=asseminit(fem,...
	'context','local',...
	'init',   struct('sd',{{{{'(x.^2+y.^2).^0.5-0.6'}}}},'ind',{1}));

% Solve dynamic problem
fem.sol=femtime(fem,...
	'tlist',  0:0.01:0.18,...
	'atol',   0.001,...
	'rtol',   0.01,...
	'jacobian','equ',...
	'mass',   'full',...
	'ode',    'ode15s',...
	'odeopt', struct('InitialStep',{[]},'MaxOrder',{5},'MaxStep',{[]}),...
	'out',    'sol',...
	'stop',   'on',...
	'report', 'on',...
	'context','local',...
	'sd',     'off',...
	'Epoint', 'gauss2',...
	'Nullfun','flnullorth',...
	'Tpoint', 'gauss2',...
	'Solcomp',1);

figure(1);
% Plot solution
postplot(fem,...
	'context','local',...
	'contdata',{'u','cont','on'},...
	'contlevels',[0 0],...
	'contstyle','color',...
	'contlabel','off',...
	'contmaxmin','off',...
	'contbar','off',...
	'contmap','cool',...
	'geom',   'on',...
	'geomcol','bginv',...
	'view',   [0 90],...
	'title',  'Time=0',...
	'renderer','zbuffer',...
	'solnum', 1,...
	'axisvisible','on')

figure(2)                                                                      
% Plot solution                                                                 
postplot(fem,...                                                                
        'context','local',...                                                   
        'contdata',{'u','cont','on'},...                                        
        'contlevels',[0 0],...                                                  
        'contstyle','color',...                                                 
        'contlabel','off',...                                                   
        'contmaxmin','off',...                                                  
        'contbar','off',...                                                     
        'contmap','cool',...                                                    
        'geom',   'on',...                                                      
        'geomcol','bginv',...                                                   
        'view',   [0 90],...                                                    
        'title',  'Time=0.1',...                                                
        'renderer','zbuffer',...                                                
        'solnum', 11,...                                                         
        'axisvisible','on')                                                     
                    
figure(3);

% Plot solution                                                                 
postplot(fem,...                                                                
        'context','local',...                                                   
        'contdata',{'u','cont','on'},...                                        
        'contlevels',[0 0],...                                                  
        'contstyle','color',...                                                 
        'contlabel','off',...                                                   
        'contmaxmin','off',...                                                  
        'contbar','off',...                                                     
        'contmap','cool',...                                                    
        'geom',   'on',...                                                      
        'geomcol','bginv',...                                                   
        'view',   [0 90],...                                                    
        'title',  'Time=0.14',...                                                
        'renderer','zbuffer',...                                                
        'solnum', 15,...                                                       $
        'axisvisible','on')                     

figure(4)

postplot(fem,...                                                                
        'context','local',...                                                   
        'contdata',{'u','cont','on'},...                                        
        'contlevels',[0 0],...                                                  
        'contstyle','color',...                                                 
        'contlabel','off',...                                                   
        'contmaxmin','off',...                                                  
        'contbar','off',...                                                     
        'contmap','cool',...                                                    
        'geom',   'on',...                                                      
        'geomcol','bginv',...                                                   
        'view',   [0 90],...                                                    
        'title',  'Time=0.16',...                                              $
        'renderer','zbuffer',...                                                
        'solnum', 17,...                                                       $
        'axisvisible','on')