global q u V k r 
q = 3;
u = .5;
V = 10;
k = 100;
r = 5;

tt = 0:0.1:10;
x0= [1,2];

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

plot(t,x)
xlabel('time')
ylabel('abundance')
title('Population and food source over time')
legend('Population','Food')