clear all
% set up global variables for use in the DE routine.
global V K

V = 2;
K = 1;

t = (0:.1:10); %increments of time
s0 = 8; %initial substrate
p0 = 0; %initial product

% initial state for use in DE solver.

X0  = [s0;p0];


[t,Xt] = ode45('KaraDismukeDE',t,X0);
    
plot(t,Xt)
title('Enzyme Substrate and Product Changes')
xlabel('Time (s)');
ylabel('Enzyme Role');
legend('s','p')