# GP oscillator

params eps=0.05,sq=1.58,rho=50,sig=1,del=20

x'= del*((sq*(1+rho*x*x)/(1+x*x +sig*y*y))-x)
y'= del*eps*sq*(1+rho*x*x)/(1+x*x) -y


@ TOTAL=30,DT=.05,xlo=0,xhi=100,ylo=0,yhi=100
@ xplot=x,yplot=y
init x=1,y=1

done