# Model from Hervagault et al., in Dynamics of Biochemical Systems (Ricard # and Cornish-Bowden, Eds.), Plenum: New York, 1984, pp. 157-169. ds/dt = s0-s - lambda*exp(-kappa*s)*s/(s+1) param s0=4, lambda=0, kappa=1 s(0)=4 done