I want to use matcont for a simple example:
p=; ap=1;
opt=contset;
opt=contset(opt,'MaxNumPoints',200);
opt=contset(opt,'Singularities',2);
opt=contset(opt,'Backward',1);
=init_EP_EP(@ab,,p,ap);
=cont(@equilibrium,x0,,opt);
I can see a LP point, now I want to do two bifurcation analysis. this is my code:
I want to use matcont for a simple example:
p=; ap=1;
opt=contset;
opt=contset(opt,'MaxNumPoints',200);
opt=contset(opt,'Singularities',2);
opt=contset(opt,'Backward',1);
=init_EP_EP(@ab,,p,ap);
=cont(@equilibrium,x0,,opt);
I can see a LP point, now I want to do two bifurcation analysis. this is my code:
x1 = x(1:end-1,s(2).index);
n1 = x(end,s(2).index);
p=;
ap=;
= init_LP_LP(@ab, x1, p, ap)
=cont(@equilibrium,x0,v0,opt);
I can see an error which I can not solve it:
??? Error using ==> vertcat
CAT arguments dimensions are not consistent.
Error in ==> cjac at 12
fullx = ;
Error in ==> init_LP_LP at 71
jac = cjac(lpds.func,lpds.Jacobian,x,p,lpds.ActiveParams);
Error in ==> ab4 at 38
= init_LP_LP(@ab, x1, p, ap)
try replacing the last line (after initializing the LP continuation) for:
=cont(@limitpoint,x0,v0,opt);
regards,
gustavo