Hi Everyone,

I had some trouble when going through the PyCont documentation, and wondered if anyone would be able to help me out.

I have a system of ODEs and I'm looking to map out (as a function of two parameters) where the system has stable oscillations.

Following the catalytic oscillator example, I am able to get the EQ1 curve that has (for one parameter) the Hopf bifurcation at the upper and lower bound of the first parameter.

1 PCargs = args(name='EQ1', type='EP-C') # as many arguments can be supplied here as desired 2 PCargs.freepars = ['q2'] # rest of the arguments created like a struct 3 PCargs.StepSize = 1e-3 4 PCargs.MaxNumPoints = 175 5 PCargs.MaxStepSize = 1e-2 6 PCargs.LocBifPoints = ['LP','H'] 7 8 PyCont.newCurve(PCargs) 9 PyCont['EQ1'].forward()

I then try to free the second parameter and compute the hopf curve, but I don't see where I could get a plot (x is first parameter, y is second) of where the solution has oscillations.

[similar to:

1 PCargs = args(name='HO1', type='H-C2') 2 PCargs.initpoint = 'EQ1:H2' 3 PCargs.freepars = ['q2', 'k'] 4 PCargs.MaxStepSize = 1e-2 5 PCargs.LocBifPoints = ['GH', 'BT', 'ZH'] 6 PCargs.MaxNumPoints = 100 7 8 PyCont.newCurve(PCargs) 9 PyCont['HO1'].forward()

Is this curve even the right thing to be calculating? Is there an easy way to see or get the data with respect to the two parameters?

Thanks!

Peter