animalcmds

 command templates for animal; run cell-by-cell;

0001 % command templates for animal; run cell-by-cell;
0002 close all; clear all; 
0003 %% trivial branch with small dlammax and lam=11.8 to find bif-points
0004 p=[];p=animalinit(p); p.smod=0; p=cont(p); 
0005 %% first bif. 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 %% 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 %% bif-diagram
0011 clf(3);plotbra(p,3,1,'ms',10,'lw',4,'cl','b');
0012 plotbra(p0,3,1,'ms',10,'lw',4,'cl','b');
0013 plotbra(q,3,1,'ms',10,'lw',4,'lab',[25]);
0014 xlabel('\lambda');ylabel('||u_1-1||_{L^1}/|\Omega|');
0015 %% plots
0016 plotsolf('q','p25',4,1,2); axis image; colormap copper
0017 plottauf('p','bp1',5,1,2);axis image; colormap copper;
0018 plottauf('p','bp2',6,1,2);axis image; colormap copper;
0019