global V s K p t0 = 0; %initial time t1 = 10; % final time s0 = 8; % initial substrate conc p0 = 0; % initial prouct conc l0 = [s0;p0]; V = 2; K = 1; [t,lt] = ode45('LaurenMageeDE',[t0,t1],l0); figure plot(t,lt(:,1),'r') hold on plot(t,lt(:,2)),title('Substrate and Product vs. Time')