# File two_component_pathway.ode
#Model of two-component signalling pathway 
# Figure 6.3

par k1=5, km1=1, k2=6, k3=2, RT=2, PT=8

#model equations
R=RT-RL
Pstar=PT-P

#set time-varying profile for ligand level
LT=heav(t-1.001) - heav(t-3.001)

RL'=k1*R*LT - km1*RL
P'=-k2*RL*P+k3*Pstar

# initial condition
init RL=0, P=0

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

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

done