function dy=glu_growth_cons(t,y)

global a1 a2 Km n;

N = y(2); %Number of cells
Glu = y(1); %Levels of glucose
%Do stability analysis

fhglu = (Glu.^n)./ (Glu.^n + Km.^n);

%dGludt = -a2*N;
dGludt = -a2*Glu*N;
dNdt = a1*N*fhglu; %cell growth depends on the amount of glucose present in the medium. Stops growing when glucose is used up.
 %glucose is bein gused up as a function of the population of cells

dy = [dGludt; dNdt];