# Activator-inhibitor system
#
dx/dt = (a+b*x*x)/(1+x*x+r*y) -x
dy/dt = eps*(c*x+yo-y)

par a=1, b=5, c=4, r=1, yo=0, eps=.1

done