% this is the lac operon in KS edition 2
function F_prime=deRHS(t,y)
global aA gA aB gB aP gP aM gM aL gL bA bL K K1 KA KL KL1 KLe Le

% variable are

M = y(1);
P = aP*M/gP;
B = aB *M / gB ;
A = y(2);

Mp = aM*(1+K1*A^2) /( K + K1*A^2) - gM *M;
Apr = aL *P *Le/( KLe + Le) - bA *B*A/( KA + A) - gA *A;


F_prime = [Mp;Apr];
