\[ \left \{x'(t)=x(t) (a (p x(t)+q y(t))+\alpha ),y'(t)=y(t) (b (p x(t)+q y(t))+\beta )\right \} \] ✗ Mathematica : cpu = 300.058 (sec), leaf count = 0 , timed out
$Aborted
✓ Maple : cpu = 10.175 (sec), leaf count = 147
\[ \left \{ [ \left \{ x \left ( t \right ) =0 \right \} , \left \{ y \left ( t \right ) ={\frac {\beta }{{{\rm e}^{-\beta \,t}}{\it \_C1}\,\beta -qb}} \right \} ],[ \left \{ x \left ( t \right ) ={\it ODESolStruc} \left ( {\it \_b} \left ( {\it \_a} \right ) ,[ \left \{ \left ( \left ( {\it \_b} \left ( {\it \_a} \right ) \right ) ^{{\frac {-a-b}{a}}}{\frac {\rm d}{{\rm d}{\it \_a}}}{\it \_b} \left ( {\it \_a} \right ) - \left ( {\it \_b} \left ( {\it \_a} \right ) \right ) ^{-{\frac {b}{a}}} \left ( ap{\it \_b} \left ( {\it \_a} \right ) +\alpha \right ) \right ) {{\rm e}^{-{\frac {{\it \_a}\, \left ( a\beta -\alpha \,b \right ) }{a}}}}+{\it \_C1}=0 \right \} , \left \{ {\it \_a}=t,{\it \_b} \left ( {\it \_a} \right ) =x \left ( t \right ) \right \} , \left \{ t={\it \_a},x \left ( t \right ) ={\it \_b} \left ( {\it \_a} \right ) \right \} ] \right ) \right \} , \left \{ y \left ( t \right ) ={\frac {- \left ( x \left ( t \right ) \right ) ^{2}ap-\alpha \,x \left ( t \right ) +{\frac {\rm d}{{\rm d}t}}x \left ( t \right ) }{qx \left ( t \right ) a}} \right \} ] \right \} \]