clear all % set up global variables for use in the DE routine. global r1 rn1 r2 rn2 r3 t0 = 0; %initial time t1 = 10; % final time c0 = 0; % initial ?? conc b0 = 0; % initial nitrogen conc a0 = 5; % initial yeast concentration % initial state for use in DE solver. r1 = 1; rn1 = 1; r2 = 1; rn2 =1; r3 = 1; S0 = [a0;b0;c0]; [t,St] = ode45('Lsquareddynamics',[t0,t1],S0); figure %ylim ([]) plot(t,St(:,1)) hold on plot(t,St(:,2)) hold on plot(t,St(:,3)),title('Reaction Models') xlabel ('time (s)') ylabel ('concentrations (mM)')