# Brussellator
#

dx/dt = a - (b+1)*x + x*x*y
dy/dt = b*x - x*x*y

par a=1, b=1

done