###########################################################
#
#  Model glykolýzy
#
#  parametry a,b > 0 
#
#
#  Kritická hodnota superkritické Hopfovy bifurkace 
#  singulárního bodu [b,b/(a+b^2)]
#  AUTO spusťte pro b v intervalu <0,1>, zobrazujte y v <0,6>
#  
#
###########################################################
#
#       parameter values defining model
#
par a=0.2,b=0.8
#
#       definice funkcí
#
f1(x,y)=-x+a*y+x^2*y
f2(x,y)=b-a*y-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=4,ylo=0,yhi=2,dt=0.01,maxstor=20000,xp=x,yp=y
#
done
