clear all;
close all;

global a1 a2 Km n;

a1 = 1;
a2 = 0.5;
%Km = 1;
n = 2;
%Go = 10;
% Initial conditions: Po, Gluo: We let Glu o = 10, when Glu = 3, then we
% stop growing and saturate...
i=1;
for Km = 1:3
[T,Y] = ode45('glu_growth_cons1', [0:0.1:200], [1 0.1]);
Glu(:,i) = Y(:,1);
N(:,i) = Y(:,2);
i=i+1;
end


subplot(1,2,1);plot(T,Glu);
title('[Glucose]');
legend('Th=1','Th=2', 'Th=3');
subplot(1,2,2);plot(T,N);
title('Population');