Additionally, the authors studied system (REF ) in \(\bar{\Omega }=\lbrace \left(u,v\right)\in \mathbb {R}^2,~u>0,~v\ge 0\rbrace \)  and they proved that \(\varphi \)  (REF ) is a diffeomorphism which preserve the orientation of time [1]}, [2]}. The \(u\) -nullclines of system (REF ) are given by \(u=0\)  and \(v=\left(u+A\right)\left(1-u\right)\left(u-M\right)/Q\) , while the \(v\) -nullclines are given by \(v=0\)  and \(v=u\) . Hence, the equilibrium point of system (REF ) are \(\left(0,0\right)\) , \(\left(1,0\right)\) , and the positive equilibrium point(s) \(\left(u^*,v^*\right)\)  with \(v^*=u^*\)  and where \(u^*\)  is determined by the solution(s) of
\(u^3-T\left(A,M\right)u^2-L\left(A,M,Q\right)u+AM =0\,,\) 
