#Figure 4.11
#.ode file for plotting the functions in Figure () where
#x is the dummy variable for plotting the auxiliary functions.
#x is the log of the calcium concentration
#initial value of x
x(0)=-80
#the equation
dx/dt=1
#parameters
par IP1=0.5 ,IP2=0.65  ,IP3=0.8
par d1=0.027, d2=0.94, d3=0.2
par k85=0.21
#the equilibrium activation and inactivation 
#functions to plot
aux c=exp(x) # this makes c the calcium concentration
aux yinf1=(d1+k85*IP1)/(d1+k85*IP1+d2*d4*c+IP1*d4*c)
aux yinf2=(d1+k85*IP2)/(d1+k85*IP2+d2*d4*c+IP2*d4*c)
aux yinf3=(d1+k85*IP3)/(d1+k85*IP3+d2*d4*c+IP3*d4*c)

aux x01=8000*c*IP1/((52+400*IP1)*(1.64+20*c))
aux x02=8000*c*IP2/((52+400*IP2)*(1.64+20*c))
aux x03=8000*c*IP3/((52+400*IP3)*(1.64+20*c))
done