A continuation run of limit cycles has the standard form
[x,v,s,h,f]=cont(@limitcycle,x0,v0,opt);
Here x is a matrix in which each column corresponds to a computed orbit. The last element of each column contains the value of the active parameter.
CL/SMALL>_MATCONT provides only limited possibilities to visualise the output. Namely, it provides the routine plotcycle.m which is to be called
in the form
plotcycle(x,v,s,e)
where e is either a two- or a three-dimensional row vector.
- If e=[e1,e2] one must have
and a
2D plot will be generated in which the projections of the orbits on the
coordinate plane are shown.
The first and last curves, and the curves with singulaties are in red.
- If e=[e1,e2,e3] then e1 must be the index of the active parameter, i.e. e1=size(x,1) and
one must have
Then a
3D plot will be generated in which the projections of the orbits on the
coordinate plane are shown versus the
active parameter. The first and last curves, and the curves with singulaties are in red.