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

%linear tradeoff
r=a*x(1)+b;
xeq=g/x(1);
yeq=(r-q*g/x(1))/(x(1)+q);
xdot=sigma*(a*-yeq);



