# brusselator.ode
# 17 Oct, 2009

x' = a + x^2*y - (b+1)*x
y' = b*x - x^2*y

par a=1, b=1
init x=1, y=1

# Settings.  You need a space after the @

# total time of integration
@ total=60

# allocate memory, set bounds and time step
@ maxstor=500000
@ bounds=1e7
@ dt=.05

# set up the initial screen
@ xp=t, yp=x, xlo=0, xhi=60, ylo=0, yhi=5

done