###########################################################
#
#  BRUSELÁTOR - Chemický model oscilace koncentrací x a y
#
#  parametry a,b > 0 
#
#
#  Kritická hodnota superkritické Hopfovy bifurkace b=1+a^2
#  singulárního bodu [a,b/a]
#  AUTO spusťte pro b v intervalu <1,3>, zobrazujte y v <1,3>
#  l1=-3/4 v singulárním bodě [1,2] pro a=1, b=2
#
#
###########################################################
#
#       parameter values defining model
#
par a=1,b=1
#
#       definice funkcí
#
f1(x,y)=a-(b+1)*x+x^2*y
f2(x,y)=b*x-x^2*y
#
#      diferenciální rovnice
#
x'=f1(x,y)
y'=f2(x,y)
#
#      počáteční podmínky
#
x(0)=0
y(0)=0
#
#
# Numerické parametry:
#
#
@ total=100,xlo=0,xhi=2,ylo=0,yhi=4,dt=0.01,maxstor=20000,xp=x,yp=y
#
done
