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

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

%initial state for us in DE solver

X0 = [s0;p0];

V = 2;
K = 1;

[t,Xt] = ode45('nataliewilliamsDE',(t0:.01:t1),X0);

plot(t,Xt(:,1), t,Xt(:,2),'--')
title('Enzymatic Utlization and Productivity');
xlabel('Time');
ylabel('Enzyme Function');
legend('s', 'p')