\[ \left \{4 x'(t)+2 x(t)+9 y'(t)+31 y(t)=e^t,3 x'(t)+x(t)+7 y'(t)+24 y(t)=3\right \} \] ✓ Mathematica : cpu = 0.165582 (sec), leaf count = 83
\[\left \{\left \{x(t)\to \frac {1}{442} e^{-4 t} \left (-442 \left (c_1+c_2\right ) \sin (t)+442 c_1 \cos (t)+31 e^{4 t} \left (17 e^t-78\right )\right ),y(t)\to \left (2 c_1+c_2\right ) e^{-4 t} \sin (t)+c_2 e^{-4 t} \cos (t)-\frac {2 e^t}{13}+\frac {6}{17}\right \}\right \}\]
✓ Maple : cpu = 0.077 (sec), leaf count = 62
\[ \left \{ \left \{ x \left ( t \right ) ={{\rm e}^{-4\,t}}\sin \left ( t \right ) {\it \_C2}+{{\rm e}^{-4\,t}}\cos \left ( t \right ) {\it \_C1}-{\frac {93}{17}}+{\frac {31\,{{\rm e}^{t}}}{26}},y \left ( t \right ) ={\frac { \left ( \left ( -221\,{\it \_C1}-221\,{\it \_C2} \right ) \cos \left ( t \right ) +221\,\sin \left ( t \right ) \left ( {\it \_C1}-{\it \_C2} \right ) \right ) {{\rm e}^{-4\,t}}}{221}}-{\frac {2\,{{\rm e}^{t}}}{13}}+{\frac {6}{17}} \right \} \right \} \]