# standard nonlinear oscillator
R=x^2+y^2
x'=x*(a-R)-y*(1+q*R)
y'=y*(a-R)+x*(1+q*R)
par a=1,q=1
done