string: To further investigate the behaviour of the stability of the system while multiple parameters were varied, we constructed the characteristic polynomial of the Jacobian evaluated at each steady state (J ) and determined the stability of each steady state using the Routh-Hurwitz criterion, a preferred method for nonlinear ODE systems, which uses the coefficients of J and hence bypasses solving the roots of each polynomial [45,46].