ac6cmds

 command templates for AC over hexagon; run cell-by-cell;

0001 % command templates for AC over hexagon; run cell-by-cell;
0002 close all; clear all; 
0003 p=[];p=ac6init(p);p.nsteps=40;p=cont(p); % AC with x-dependent-DBC
0004 %%
0005 q=swibra('p','bp1','q',-0.1);q=cont(q);
0006 r1=swibra('p','bp2','r1',-0.1);r1=cont(r1);
0007 r2=loadp('r1','p1','r2'); r2.ds=-r1.ds; r2=cont(r2);
0008 %% plot first bif-diagram
0009 figure(3);clf; cmp=2; 
0010 plotbra(p,3,cmp,'ms',8,'lw',4,'cl','b','lab',40); 
0011 plotbra(q,3,cmp,'ms',8,'lab',[10],'lw',4);
0012 plotbra(r1,3,cmp,'ms',8,'lw',4,'cl','r');
0013 plotbra(r2,3,cmp,'ms',8,'lw',4,'cl','r');
0014 xlabel('\lambda');ylabel('||u||_2');
0015 %%
0016 figure(4);clf(4); cmp=1; 
0017 plotbra(p,4,cmp,'ms',8,'lab',[10],'lw',4, 'cl', 'b');
0018 plotbra(r1,4,cmp,'ms',8,'lab',[5],'lw',4,'cl','r');
0019 plotbra(r2,4,cmp,'ms',8,'lab',[5],'lw',4,'cl','r');
0020 axis([0.8 1.5 1.1 1.4]);
0021 xlabel('\lambda');ylabel('max(u)');
0022 %%
0023 plotsolf('p','p1',5,1,4);
0024 plotsolf('q','p10',6,1,1);view(-20,60);
0025 plotsolf('p','p30',7,1,1); view(-20,60);
0026 %% with homogenenous DBC
0027 close all;clear all;
0028 ph=[];ph=ac6init(ph);ph.bcf=@ac6bcf;ph=cont(ph);