global cds homTds homds xhom vhom shom

p=[-0.4;1.1;-0.1];
ap=[1 2];
x=xhom(1:end-1,shom(2).index);
p2=p;p2(homds.ActiveParams)=xhom(end,shom(2).index);
nu=homds.nu;
ns=homds.ns;
nphase=homds.nphase;


[x0,v0]=init_homT_homT(@Henon,x,nphase,nu,ns, p2, ap,1);


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

[xhomT,vhomT,shomT,homT,fhomT]=cont(@homoclinicT,x0,[],opt);

 cpl(xhomT,vhomT,shomT,[26 25])
 pause
opt=contset(opt,'Backward',0);

[xhomT,vhomT,shomT,homT,fhomT]=cont(@homoclinicT,x0,[],opt);

 cpldel(xhomT,vhomT,shomT,[26 25])