
function F_prime=deRHS(t,y)
global vm vp km k3 k1 k2 J K A

M = y(1);
P = y(2);
q = 2/(1+sqrt(1+8*K*P));
P1 = q*P;
P2 = ((1-q)/2)*P;

Xp = vm/(1+(P2/A)^2) - km*M;

Yp = vp*M-(k1*P1+2*k2*P2)/(J+P)-k3*P;

F_prime = [Xp;Yp];
