hettestrun.asv

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

C=[-0.508776078,  0.5155,  .5354,  .5021,  .5573,  .4643,  .6167,  .3560,  .7689,  .0380,  1.0657,  -.8020,  .8098,  -1.0941,  .0027,  -.7955,   -.3299,   -.6318 ,  -.4505,   -.5641,  -.4949,  -.5377, -.5113, -.5280,  -.5178,   -0.50877460786069 ;0.53604733513342,.5354,.5021,.5572,.4641,.6167,.3560 ,.7689,.0380, 1.0657, -.8019 ,.8098,-1.0941,.0027,-0.7955,-.3299,-.6318, -.4505 ,-.5641 ,-.4950, -.5379, -.5113, -.5275, -.5178,-.5238,-0.50877460786069];
 p=[0.3;-1.03;-0.1];ap=[2 3];%p=[a,b,r]
 %0.53604733513342 %
 C,pause
 disp('>> opt = contset; ');
 opt = contset;
disp('>> [x0,v0]=init_HE_HE(@Henon,C, p, ap,J,eps0,eps1);');%(mapfile, C, p,ap,J, eps0,eps1)
[x0,v0]=init_HE_HE(@Henon,C, p, ap,0.020568,0.017528,1);
%   opt=contset; opt=contset(opt,'Maxstepsize',5); 
%   opt=contset; opt=contset(opt,'Minstepsize',4); 
opt=contset; opt=contset(opt,'MaxNumPoints',2);

opt=contset(opt,'Singularities',0);
opt=contset(opt,'Backward',0);
opt = contset(opt,'AutDerivative',0);

[x1,v1,s1,h1,f1]=cont(@heteroclinic,x0,[],opt);
%cpl(x1,v1,s1,[3 1])