function F_prime=deRHS(t,y)
global K1 K2 c1 c2 

s = y(1);
u = y(2);
v = y(3);

Sp = (1-s) - c1*s*u/(K1+s);
Up = c1*s*u/(K1+s) - c2*u*v/(K2+u)-u;
Vp =  c2*u*v/(K2+u)-v;

F_prime = [Sp;Up;Vp];