acgccmds

 command templates for AC with global coupling; run cell-by-cell

0001 % command templates for AC with global coupling; run cell-by-cell
0002 close all; clear all; 
0003 %% first branch
0004 p=[];p=acgcinit(p,1);p=meshref(p);p=cont(p); 
0005 %% 2nd branch
0006 q=[]; q=acgcinit(q,2);q.nsteps=40;q=cont(q); 
0007 %% postprocessing
0008 figure(3);clf;
0009 plotbra(p,3,1,'lab',[1,20],'lw',5) 
0010 plotbra(q,3,1,'lab',[1,10],'cl','b','lw',5); 
0011 axis([-0.5 1 -4 18]); xlabel('\lambda');ylabel('<u>');
0012 %%
0013 plotsolf('p','p1',4,1,3);plotsolf('p','p20',5,1,3);
0014 %%
0015 plotsolf('q','p1',6,1,3);
0016 plotsolf('q','p10',7,1,3);
0017 %%
0018 plotsolf('p','p20',9,1,3);
0019 %%
0020 plotbra(p,10,2,'lab',[1,15,20],'lw',5)