\[ \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.066 (sec), leaf count = 0 , timed out
$Aborted
✓ Maple : cpu = 12. (sec), leaf count = 181
\[ \left \{ [ \left \{ x \left ( t \right ) =0 \right \} , \left \{ y \left ( t \right ) ={\frac {\beta }{{{\rm e}^{-\beta \,t}}{\it \_C1}\,\beta -bq}} \right \} ],[ \left \{ x \left ( t \right ) ={\it ODESolStruc} \left ( {\it \_b} \left ( {\it \_a} \right ) ,[ \left \{ \left ( {\it \_b} \left ( {\it \_a} \right ) \right ) ^{-{\frac {a+b}{a}}}{{\rm e}^{-{\frac {{\it \_a}\, \left ( a\beta -b\alpha \right ) }{a}}}}{\frac {\rm d}{{\rm d}{\it \_a}}}{\it \_b} \left ( {\it \_a} \right ) -{{\rm e}^{-{\frac {{\it \_a}\, \left ( a\beta -b\alpha \right ) }{a}}}}{\it \_b} \left ( {\it \_a} \right ) \left ( {\it \_b} \left ( {\it \_a} \right ) \right ) ^{-{\frac {b}{a}}}ap-{{\rm e}^{-{\frac {{\it \_a}\, \left ( a\beta -b\alpha \right ) }{a}}}} \left ( {\it \_b} \left ( {\it \_a} \right ) \right ) ^{-{\frac {b}{a}}}\alpha +{\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 \} \]