acgcdemo

 demo script for AC with global coupling

0001 % demo script for AC with global coupling
0002 close all; clear all; 
0003 fprintf('\nInit Allen-Cahn with global coupling, first branch:\n'); 
0004 fprintf('(here branch selected by initial guess via icsel-argument of acgcinit\n');
0005 p=[];p=acgcinit(p,1);
0006 mypause('mesh refine starting point of branch and run continuation'); 
0007 p=meshref(p);p=cont(p); 
0008 mypause('\nInit and run Allen-Cahn with global coupling, second branch:'); 
0009 q=[]; q=acgcinit(q,2);q.nsteps=40;q=cont(q); 
0010 mypause('Plot bifurcation diagram');
0011 figure(3);clf;
0012 plotbra(p,3,1,'lab',[1,20],'lw',5) 
0013 plotbra(q,3,1,'lab',[1,10],'cl','b','lw',5); 
0014 axis([-0.5 1 -4 18]); xlabel('\lambda');ylabel('<u>');
0015 mypause('Plot some solutions'); 
0016 plotsolf('p','p1',4,1,3);plotsolf('p','p20',5,1,3);
0017 plotsolf('q','p1',6,1,3);plotsolf('q','p10',7,1,3);