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