% C=[-1.69449184777905 -1.6837643 -1.65949674 -1.581150109 -1.337009429 -0.659661307 0.547365455 0.062065530 -1.009338894 -1.480864536 -1.632164287;-1.69449184777905 -1.65948946 -1.581161931 -1.337017143 -0.659669162 0.547301706 0.062074908 -1.009383240 -1.480879089 -1.632160728 -1.676628792];
%  p=[-0.4;1.1;-0.1];ap=[2];%p=[a,b,r] *perfect
%  C,pause
% % C=[-0.3144706727040 -1.04920163606 -0.640540809357 0.392687493818 0.324301394255 -0.718447548185 -1.02365494748  -0.588597191062 0.437429350429 0.234771666962 -0.711565274823 -0.929086704063 -0.403516106938 0.547353044603 -0.0659223493914 -0.759037581779 -0.713901008241 0.010375914355 0.528812106379  -0.490774367408 -0.954295085512 -0.656923612169     0.279142457533  0.410479099965 -0.664676612665 -1.03320009562 -0.658206686637 0.352622134037 0.370238275605 -0.709806378525 -1.05518838715 -0.6643179710 ;0.755102259927  -0.640540809357 0.392687493818 0.324301394255 -0.718447548185 -1.02365494748 -0.588597191062 0.437429350429 0.234771666962 -0.711565274823 -0.929086704063 -0.403516106938 0.547353044603 -0.0659223493914 -0.759037581779 -0.713901008241 0.010375914355 0.528812106379 -0.490774367408 -0.954295085512 -0.656923612169 0.279142457533 0.410479099965 -0.664676612665  -1.03320009562 -0.658206686637 0.352622134037 0.370238275605 -0.709806378525 -1.05518838715 -0.66431797106 0.364875727378 ];
% % C=[-0.3144706727040 -1.04920163606  -1.02365494748   -0.929086704063  -0.713901008241  0.279142457533   0.352622134037  ;0.755102259927  -0.640540809357  -0.588597191062  -0.403516106938  0.010375914355   0.410479099965  0.370238275605  ];
%  
% 
%  disp('>> opt = contset; ');
%  opt = contset;
% disp('>> [x0,v0]=init_HO_HO(@Henon,C, p, ap,J);');
% [x0,v0]=init_HO_HO(@Henon,C, p, ap,1);
% opt=contset; opt=contset(opt,'MaxNumpoints',41);
% 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);
% 
% for i=1:1
%     %cpl(xhom,vhom,shom,[25 2*(i-1)+1])
%     cpl(xhom,vhom,shom,[25 11])
% end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%46 NOV

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];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(@Henon,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);

for i=1:1
    %cpl(xhom,vhom,shom,[25 2*(i-1)+1])
    cpl(xhom,vhom,shom,[23 11])
end