thenon7.m

来自「计算动力学系统的分岔图」· M 代码 · 共 15 行

C=[-1.62114638486435 -1.562 -1.443  -1.0987856 -0.27959455684736  0.62285460449110 -0.48255078907068 -1.24433960857123 -1.51139945471518 -1.59072793547249;-1.62114638486435 -1.443   -1.0987856    -0.27959455684736   0.62285460449110 -0.48255078907068  -1.24433960857123 -1.51139945471518 -1.59072793547249 -1.61409645976328];
 
p=[-0.4;1.03;-0.1;0];ap=[2];%p=[a,b,r] 
disp('>> opt = contset; ');
opt = contset;
disp('>> [x0,v0]=init_HO_HO(@Henon,C, p, ap,J);');
[x0,v0]=init_HO_HO(@Ghmap,C, p, ap,1);
opt=contset; opt=contset(opt,'MaxNumpoints',35);
opt=contset(opt,'Singularities',1);
opt=contset(opt,'MinStepsize',0.001);
opt=contset(opt,'Backward',1);
%opt = contset(opt,'AutDerivative',1);
[xhom,vhom,shom,hhom,fhom]=cont(@homoclinic,x0,[],opt);
cpl(xhom,vhom,shom,[23 11])