global r1 rmin1 r2 rmin2 r3

t0 = 0; %initial time
t1 = 5; %final time
a0 = 2; %initial alphaketo
b0 = 1; %initial glutamate
c0 = 0; %initial gluatmine

%initial state for use in DE solver

X0 = [a0;b0;c0];

r1 = 1;
rmin1 = 1;
r2 = 2;
rmin2 = 1;
r3 = 1;

[t,Xt] = ode45('presentation_ka_nw', [t0,t1],X0);

plot(t,Xt(:,1), t,Xt(:,2),'rd-', t,Xt(:,3),'g^');
title('Protein Productivity');
xlabel('Time');
ylabel('Protein');
legend('a','b','c');