function dxdt = michaelis_menten_lineweaver_burk(t,x)
% USAGE: dxdt = michaelis_menten_lineweaver_burk(t,x)
global K Vmax
% initialize the time derivative for the model
dxdt = zeros(size(x));
% rename substrate (s) and product (p) for ease of reading
s = x(1);
p = x(2);
% the time derivatives of substrate (1) and product (2) are
% negatives of each other and defined by the mm model
dxdt(1) = -Vmax*s/(K+s);
dxdt(2) = Vmax*s/(K+s);
% differential equation values computed. function completed.
return