# file bistable_example1.ode
# Bistable network in Figure 4.6
#
par  k1=20, k3=5, k2=20, k4=5, n1=2, n2=2
#
S1'=k1/(1+S2^n1) - k3*S1
S2'=k2/(1+S1^n2) - k4*S2
#
init S1=1,S2=0.1
#
done