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

t0  = 0;     %initial time
t1  = 10;    % final time
s0  = 8;     %initial substrate
p0  = 0;     %initial product

% initial state for use in DE solver.

S0  = [s0;p0];

V=2;
K=1;



[t,St] = ode45('AlyssaGomesDE',(t0:0.1:t1),S0);
    
plot(t,St(:,1), t,St(:,2),'-o')
title('Enzyme Substrate Versus Product Amounts');
xlabel('Time');
ylabel('Substrate Activity');


legend('s=substrate activity','p=product activity')