function [ dXdt ] = TessaAMorrisDE(t,X)
global K V 

dXdt = zeros(size(X));

s = X(1); % substrate

dsdt = -V*s/(K+s);
dpdt = V*s/(K+s);

dXdt(1) = dsdt;
dXdt(2) = dpdt;

end