# Gene expression
 
dx/dt=1+a*x*x/(1+x*x+y*y)-b*x
dy/dt=eps*(1+a*x*x/(1+x*x+y*y))-c*y

p a=50, b=8, c=.4, eps=.1

done