lesliegower23.asv

计算动力学系统的分岔图

disp('>>>> global opt cds fpmds ');
global opt cds fpmds 
disp('>> opt = contset; ');
opt = contset;
disp('>> opt = contset(opt,''Multipliers'',1); ');
opt = contset(opt,'Multipliers',1);
disp('>> ap=13; p=[20;18;0.72;0.23;0.29;0.98;0.36;0.55;0.18;0.26;0.3;0.23;0.3;.08];');
p = [20;18;0.72;0.23;0.29;0.98;0.36;0.55;0.18;0.26;0.3;0.23;0.3;.08];ap=13;

disp('>> [x0,vO]=init_FP_FPn(@LeslieGower,[11.17052012;11.94625068;20.99523834;15.21083594], p, ap,1);');
[x0,vO]=init_FPm_FPm(@LeslieGower,[11.17052012;11.94625068;20.99523834;15.21083594], p, ap,1);
disp('>> opt=contset;opt=contset(opt,''MaxNumPoints'',900);');
opt=contset;opt=contset(opt,'MaxNumPoints',900); 
disp('>> opt=contset(opt,''Singularities'',1);');
opt=contset(opt,'Singularities',1);
disp('>> [x231,v231,s231,h231,f231]=cont(@@fixedpointmap,x0,[],opt);');
[x231,v231,s231,h231,f231]=cont(@fixedPointmap,x0,[],opt,1);
disp('>> opt=contset(opt,''Backward'',1);');
opt=contset(opt,'Backward',1);
disp('>> [x232,v232,s232,h232,f232]=cont(@@fixedpointmap,x0,[],opt);');
[x232,v232,s232,h232,f232]=cont(@fixedPointmap,x0,[],opt,1);