\[ \left \{x''(t)=\text {a1} x(t)+\text {b1} y(t)+\text {c1},y''(t)=\text {a2} x(t)+\text {b2} y(t)+\text {c2}\right \} \] ✓ Mathematica : cpu = 19.8316 (sec), leaf count = 15664
DSolve[{Derivative[2][x][t] == c1 + a1*x[t] + b1*y[t], Derivative[2][y][t] == c2 + a2*x[t] + b2*y[t]},{x[t], y[t]},t]
\[ \text {Too large to display} \] ✓ Maple : cpu = 0.149 (sec), leaf count = 457
dsolve({diff(diff(x(t),t),t) = a1*x(t)+b1*y(t)+c1, diff(diff(y(t),t),t) = a2*x(t)+b2*y(t)+c2})
\[\left \{x \left (t \right ) = c_{4} {\mathrm e}^{\frac {\sqrt {2 \sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}+2 \operatorname {a1} +2 \operatorname {b2}}\, t}{2}}+c_{3} {\mathrm e}^{-\frac {\sqrt {2 \sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}+2 \operatorname {a1} +2 \operatorname {b2}}\, t}{2}}+c_{2} {\mathrm e}^{\frac {\sqrt {2 \operatorname {a1} +2 \operatorname {b2} -2 \sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}}\, t}{2}}+c_{1} {\mathrm e}^{-\frac {\sqrt {2 \operatorname {a1} +2 \operatorname {b2} -2 \sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}}\, t}{2}}+\frac {\operatorname {b1} \operatorname {c2}}{\operatorname {a1} \operatorname {b2} -\operatorname {a2} \operatorname {b1}}-\frac {\operatorname {b2} \operatorname {c1}}{\operatorname {a1} \operatorname {b2} -\operatorname {a2} \operatorname {b1}}, y \left (t \right ) = \frac {-c_{1} \left (\operatorname {a1} \operatorname {b2} -\operatorname {a2} \operatorname {b1} \right ) \left (\sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}+\operatorname {a1} -\operatorname {b2} \right ) {\mathrm e}^{-\frac {\sqrt {2 \operatorname {a1} +2 \operatorname {b2} -2 \sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}}\, t}{2}}-c_{2} \left (\operatorname {a1} \operatorname {b2} -\operatorname {a2} \operatorname {b1} \right ) \left (\sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}+\operatorname {a1} -\operatorname {b2} \right ) {\mathrm e}^{\frac {\sqrt {2 \operatorname {a1} +2 \operatorname {b2} -2 \sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}}\, t}{2}}+c_{3} \left (\operatorname {a1} \operatorname {b2} -\operatorname {a2} \operatorname {b1} \right ) \left (\sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}-\operatorname {a1} +\operatorname {b2} \right ) {\mathrm e}^{-\frac {\sqrt {2 \sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}+2 \operatorname {a1} +2 \operatorname {b2}}\, t}{2}}+c_{4} \left (\operatorname {a1} \operatorname {b2} -\operatorname {a2} \operatorname {b1} \right ) \left (\sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}-\operatorname {a1} +\operatorname {b2} \right ) {\mathrm e}^{\frac {\sqrt {2 \sqrt {\operatorname {a1}^{2}-2 \operatorname {a1} \operatorname {b2} +4 \operatorname {a2} \operatorname {b1} +\operatorname {b2}^{2}}+2 \operatorname {a1} +2 \operatorname {b2}}\, t}{2}}-2 \operatorname {b1} \left (\operatorname {a1} \operatorname {c2} -\operatorname {a2} \operatorname {c1} \right )}{2 \operatorname {b1} \left (\operatorname {a1} \operatorname {b2} -\operatorname {a2} \operatorname {b1} \right )}\right \}\]