% 
% This script includes matlab commands to store parameters in global memory
%      for access in functions, specification of parameter values, initial 
%      conditions, and time periods, and calls to the differential equation
%      solver ode45 in order to simulate a simple enymatic reaction.
%
global K Vmax

% MM parameters

% Times at which we find the solution to the differential equation
t = 0:0.1:100;

% specify rate parameters
K = 1.192448063;
Vmax = 0.100264699

% specify initial substrate and product concentrations
z0 = [4.8;0];

% call ode45 for the mm subs-prod model
% zz has two columns for the substrate and product
[tz,zz] = ode45('michaelis_menten_lineweaver_burk',t,z0);

% open a new figure and plot the substrate and products for mm and full
% model
figure
plot(tz,[zz(:,1),zz(:,2)]),legend('MM substrate','MM product');

Vmax = 0.084210526
K = 2.006568421
[ttreat,ztreat] = ode45('michaelis_menten_lineweaver_burk',t,z0);
hold on
plot(ttreat,ztreat,'*')