chemcmds

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

0001 % command templates for chemotaxis; run cell-by-cell;
0002 close all; clear all; p=[]; p=cheminit(p); p=cont(p);
0003 %% the first 5 bifurcating branches
0004 q1=swibra('p','bp1','q1',-0.2);q1.tol=1e-6; 
0005 q1.bifchecksw=0; q1.spcalcsw=0; q1=cont(q1); 
0006 %%
0007 q2=swibra('p','bp2','q2',0.2);q2.tol=1e-6; 
0008 q2.bifchecksw=0; q2.spcalcsw=0; q2=cont(q2);
0009 %%
0010 q3=swibra('p','bp3','q3',0.5);q3.tol=1e-6; 
0011 q3.bifchecksw=0; q3.spcalcsw=0; q3=cont(q3);
0012 %%
0013 q4=swibra('p','bp4','q4',0.1);q4.tol=1e-6; 
0014 q4.bifchecksw=0; q4.spcalcsw=0; q4=cont(q4);
0015 %%
0016 q5=swibra('p','bp5','q5',-0.1);q5.tol=1e-8; 
0017 q5.bifchecksw=0; q5.spcalcsw=0; q5=cont(q5);
0018 %% plot bifurcation diagram
0019 figure(4);clf(4);cmp=1;
0020 plotbra(p,4,cmp,'ms',10,'lw',4,'cl','b');
0021 plotbra(q1,4,cmp,'ms',10,'lw',4);
0022 plotbra(q2,4,cmp,'ms',10,'lw',4);
0023 plotbra(q3,4,cmp,'ms',10,'lw',4);
0024 plotbra(q4,4,cmp,'ms',10,'lw',4,'lab',[10,25],'cl','r');
0025 plotbra(q5,4,cmp,'ms',10,'lw',4);
0026 axis([10 22 0 0.8]);xlabel('\lambda');ylabel('||u_1-1||_{L^1}/|\Omega|');
0027 %% stuff that might be useful
0028 p=loadp('p','p25','p');
0029 q1=loadp('q1','p21','q1');
0030 q2=loadp('q2','p21','q2');
0031 q3=loadp('q3','p21','q3');
0032 q4=loadp('q4','p33','q4');
0033 q5=loadp('q5','p20','q5');