clear all global V K V=2; K=1; t0=0; t1=10; s0=8; %initial substrate p0=0; %initial product X0=[s0;p0]; [t,Xt] = ode45('KarinaAlvarezDE',(t0:.1:t1),X0); plot(t,Xt(:,1),t,Xt(:,2)); title('Substrate and Product Vs. Time'); xlabel('Time') ylabel('Quantity') legend('Substrate','Product')