function dy = bacterial_nutrient_eqn(t,y)
global k alpha
dy = zeros(2,1); % need to be a column vector
dy(1) = k*y(2)*y(1);
dy(2) = -alpha*k*y(2)*y(1);