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

V=2; 
K=1;
t=(0:.1:10); %time vector

s0=8; %initial substrate
p0=0; %initial product
X0=[s0;p0];

[t,Xt] = ode45('KristenHorstmannDE',t,X0);
    
plot(t,Xt)
title ('Substrate and Product Change over Time')
xlabel('Time')
ylabel('Enzyme Function')
legend('Substrate','Product')