# bruss.ode
# the brusselator
u'=a-(b+1)*u+v*u^2
v'=b*u-v*u^2
par b=1.5,a=1
init u=1,v=1.5
done