function xdot=diseaseecology(t,x);
global a b c q g sigma

%linear tradeoff
xdot=[(a*x(3)+b).*x(1)-q*(x(1)+x(2)).*x(1)-x(3).*x(1).*x(2);x(3).*x(1).*x(2)-g*x(2);sigma*(a-x(2))];

