# file NfkB_signaling.m
# model of Nf-kB signaling pathway
# from Krishna et al. (2006) PNAS 103, pp. 10840-10845
# problem 7.8.15

#parameter values
par A=0.007
par B=954.5
par C=0.035
par del=0.029
par eps=2e-5

#model equations
Nn' = A*(1-Nn)/(eps+I) - B*I*Nn/(del+Nn)
Im' = Nn^2-Im
I' = Im-C*(1-Nn)*I/(eps+I)

# initial condition
init Nn=0.0014, Im=0.0207, I=0.0127

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

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

done