global a b epsiLon g
a = 2;
b = 0.5;
epsiLon = 1.1;
g = 5;

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

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