function [ dxdt ] = GendronFirstMatLabProjectFunc(t,x)
% this function defines the chemostat differential equation
global r1 rm1 r2 rm2 r3

dxdt = zeros(size(x));

a = x(1); %
b = x(2); %
c = x(3);

dadt  = -r1*a+rm1*b-r3*a*c;
dbdt  = r1*a-rm1*b+rm2*c-r2*b+2*r3*a*c;
dcdt  = -rm2*c+r2*b-r3*a*c

dxdt(1) = dadt;
dxdt(2) = dbdt;
dxdt(3) = dcdt;
end