% this is the lac operon in KS
function F_prime=deRHS(t,y)
global mu lambda kappa eps m m0


ee = y(1);
p = y(2);
s = y(3);

Ep = m0+p^m/(kappa^m+p^m)-eps*ee;
Pp = mu*ee*(s/(s+1)-lambda*p/(p+1));
Sp = -ee*s/(s+1);

F_prime = [Ep;Pp;Sp];