###########################################################
#
#  Model Michaelis-Mentenové
#
#                k2 > k1
#
#
#      s = koncentrace substrátu
#      c = koncentrace komplexu enzym-substrát
#
#      e = e0-c koncentrace enzymu
#
#
###########################################################
#
#       parametry 
#
par k1=1, kr1=1, k2=4, e0=0.01
#
#       definice funkcí
#
f1(s,c)=kr1*c - k1*s*(e0-c)
f2(s,c)=k1*s*(e0-c) - (k2 + kr1)*c 
#
#      diferenciální rovnice
#
s'=f1(s,c)
c'=f2(s,c)
#
#      počáteční podmínky
#
s(0)=1
c(0)=0
#
#
# další nastavení
#
@ total=40, xp=t, yp=c, xlo=0,xhi=40,ylo=0,yhi=4,dt=0.05,maxstor=20000
#
done
