function dx = nuPop(t,x)
global q u U v V k K R 

dx = zeros(size(x));

dx(1) = q*u - q*x(1) - (x(3)*v)*(x(1)/(k+x(1)))*(x(2)/(K+x(2)));
dx(2) = q*U - q*x(2) - (x(3)*V)*(x(1)/(k+x(1)))*(x(2)/(K+x(2)));
dx(3) = (x(3)*R)*(x(1)/(k+x(1)))*(x(2)/(K+x(2))) - q*x(3);

end


global q u U v V k K R 
q = 0.15;
u = 29;
U = 100;
v = 110;
V = 180;
k = 2.5; 
K = 10;
R = 15; 

tt = [0:0.01:10];
x0 = [118,20,6.67];


[t,x] = ode45('funfive',tt,x0);

plot(t,x)
xlabel('time')
ylabel('abundance')
title('u1=118')
legend('Ammonium Concentration','Glucose Concentration','Yeast Pop')