codstock13.m

计算动力学系统的分岔图

disp('>> global   opt cds fpmds  ');
global   opt cds fpmds 
disp('>> p=[200;0.5;0.5;0.5;1;2;3];ap=5; ');
p=[200;0.5;0.5;0.5;1;2;3];ap=5;
disp('>> opt = contset; ');
opt = contset;
disp('>> opt = contset(opt,''Multipliers'',1); ');
opt = contset(opt,'Multipliers',1);
opt = contset(opt,'Backward',1);
opt = contset(opt,'MaxNumPoints',1500);
opt = contset(opt,'Singularities',1);
disp('>>> curve of fixed point <<< ')
disp('>> [x0,v0]=init_FPm_FPm(@NAFstock,[72.8206;1.33341], p, ap,1);');
[x0,v0]=init_FPm_FPm(@NAFstock,[72.8206;1.33341], p, ap,1);
disp('>>[x1,v1,s1,h1,f1]=cont(@fixedpointmap,x0,[],opt); ');
[x1,v1,s1,h1,f1]=cont(@fixedpointmap,x0,[],opt);
cpl(x1,v1,s1,[3 1])