Two parameters bifurcation analysis

chan
2011-06-28
2013-01-23
  • chan
    chan
    2011-06-28

    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)

     
  • Gustavo
    Gustavo
    2011-06-28

    try replacing the last line (after initializing the LP continuation) for:

    =cont(@limitpoint,x0,v0,opt);

    regards,
    gustavo