gpcmds

 command templates for Gross-Pitaevski and vector GP; run cell-by-cell;

0001 % command templates for Gross-Pitaevski and vector GP; run cell-by-cell;
0002 close all; clear all; 
0003 %% first scalar GP
0004 p=[];p=gpinit(p);p=cont(p);
0005 %% branch plotting
0006 plotbra(p,3,3,'lab',[10,15]);
0007 %% solution plotting
0008 plotsolf('p','p1',4,1,5); % VF plot
0009 plotsolf('p','p10',5,1,5); % VF plot
0010 plotsolf('p','p10',6,10,3); % amp plot
0011 plotsolf('p','p15',7,10,3); % amp plot
0012 plotsolf('p','p20',8,13,2); % phase plot
0013 %% vector GP
0014 q=[];q=vgpinit(q);q=cont(q);
0015 %% branch plotting
0016 clf(3);
0017 plotbra(p,3,3,'lab',[1,10,15]);
0018 plotbra(q,3,3,'lab',[1,10,15],'cl','r');
0019 xlabel('\omega'); ylabel('p=max|u|/max|v| and p_1');
0020 %% solution plotting
0021 plotsolf('q','p1',10,1,5); % VF plot
0022 plotsolf('q','p10',11,1,5); % VF plot