# bruss2.ode
# two Brusselator equations
f(u,v)=a-(b+1)*u+v*u^2
g(u,v)=b*u-v*u^2
u1'=f(u1,v1)+du*(u2-u1)
v1'=g(u1,v1)+dv*(v2-v1)
u2'=f(u2,v2)+du*(u1-u2)
v2'=g(u2,v2)+dv*(v1-v2)
par b=1.5,a=1,du=0,dv=.2
init u1=1,u2=1,v1=1.5,v2=1.5
@ ntst=30,dsmin=1e-5,dsmax=0.1,parmin=0,parmax=5,autoxmin=0
@ autoxmax=5,autoymin=0,autoymax=7
done