# File problem3_7_5.ode
# Problem 3.7.5 part (a)

par v0=2, Vm1=9, Vm2=12, Vm3=15, Km1=1, Km2=0.4, Km3=3

s1'=v0 - Vm1*s1/(Km1+s1) 
s2'=Vm1*s1/(Km1+s1)-Vm2*s2/(Km2+s2) 
s3'=Vm2*s2/(Km2+s2)-Vm3*s3/(Km3+s3)


# initial condition
init s1=0.3, s2=0.2, s3=0.1

# final time
@ total=2
# set the plotting window size:
@ xlo=0, xhi=2, ylo=0, yhi=1

# run the simulation in the GUI by selecting (I)nitialconds|(G)o
# plot species s2 and s3 by selecting (G)raphic stuff|(A)dd curve and 
# assigning the y-axis

done