vkcmds

 command templates for von Karman plate; run cell-by-cell;

0001 % command templates for von Karman plate; run cell-by-cell;
0002 p=[]; p=vkinit(p); p=cont(p); % the trivial solution
0003 %% 2 primary and one secondary bifurcations
0004 q=swibra('p','bp1','q'); q.dlammax=1; q=cont(q);
0005 r=swibra('p','bp2','r'); r.dlammax=1; r=cont(r);
0006 w=swibra('r','bp1','w'); w.dlammax=1; w.smod=3;w=cont(w);
0007 %% branch plotting
0008 clf(3); cmp=2;
0009 plotbra(q,3,cmp,'ms',10,'lwst',2,'lwun',2,'cl','k');
0010 plotbra(r,3,cmp,'ms',10,'lwst',2,'lwun',2,'cl','b');
0011 plotbra(w,3,cmp,'ms',10,'lab',[3,6],'lw',5,'cl','r');
0012 axis([5 8 0 8]);xlabel('\lambda'); ylabel('||u_1||_2');
0013 %% solution plotting
0014 plotsolf('q','bp1',4,1,2);
0015 plotsolf('r','bp1',5,1,2);
0016 plotsolf('w','p3',6,1,2);
0017 plotsolf('w','p6',7,1,2);
0018 plotsolf('q','bp1',8,3,2);
0019 plotsolf('r','bp1',9,3,2);