animaldemo

 demoscript for chemotaxis model over animal shaped domain

0001 % demoscript for chemotaxis model over animal shaped domain
0002 close all; clear all; 
0003 mypause('Trivial branch with small dlammax and lam=11.8 to find bifurcation points');  
0004 p=[];p=animalinit(p); p.smod=0; p=cont(p); 
0005 mypause('First bifurcating branch');  
0006 q=swibra('p','bp1','q',0.05); q.xi=0.2; q.nsteps=5; q.dlammax=1;q=cont(q); 
0007 q.ds=0.2; q.smod=25;q.nsteps=100; q=cont(q);
0008 mypause('Continue trivial branch to left for plotting');
0009 p0=[];p0=animalinit(p0);p0.smod=0;p0.ds=-0.5; p0.lammin=7.5; p0.dlammax=1;p0=cont(p0);
0010 mypause('Plot bifurcation diagram'); 
0011 figure(3);clf;
0012 plotbra(p,3,1,'ms',10,'lw',4,'cl','b');
0013 plotbra(p0,3,1,'ms',10,'lw',4,'cl','b');
0014 plotbra(q,3,1,'ms',10,'lw',4,'lab',[25]);
0015 xlabel('\lambda');ylabel('||u_1-1||_{L^1}/|\Omega|');
0016 mypause('Plot some solutions'); 
0017 plotsolf('q','p25',4,1,2); axis image; colormap copper
0018 plottauf('p','bp1',5,1,2);axis image; colormap copper;
0019 plottauf('p','bp2',6,1,2);axis image; colormap copper;
0020