function [ dldt ] = LaurenMageeDE(t,l)
global V s K p

dldt = zeros(size(l));

s = l(1); % substrate
p = l(2); % product

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

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

end