\[ \left \{\text {x1}'(t)=a \text {x2}(t)+b \text {x3}(t) \cos (c t)+b \text {x4}(t) \sin (c t),\text {x2}'(t)=-a \text {x1}(t)+b \text {x3}(t) \sin (c t)-b \text {x4}(t) \cos (c t),\text {x3}'(t)=a \text {x4}(t)-b \text {x1}(t) \cos (c t)-b \text {x2}(t) \sin (c t),\text {x4}'(t)=-a \text {x3}(t)-b \text {x1}(t) \sin (c t)+b \text {x2}(t) \cos (c t)\right \} \] ✓ Mathematica : cpu = 0.0099564 (sec), leaf count = 798
\[\left \{\left \{\text {x1}(t)\to c_3 \cos \left (\left (\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right )+c_1 \cos \left (\left (\frac {c}{2}+\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right )+c_4 \sin \left (\left (\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right )+c_2 \sin \left (\left (\frac {c}{2}+\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right ),\text {x2}(t)\to -c_4 \cos \left (\left (\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right )-c_2 \cos \left (\left (\frac {c}{2}+\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right )+c_3 \sin \left (\left (\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right )+c_1 \sin \left (\left (\frac {c}{2}+\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right ),\text {x3}(t)\to \frac {\left (a+\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) c_4 \cos \left (\frac {1}{2} \left (c+\sqrt {4 a^2+4 c a+4 b^2+c^2}\right ) t\right )}{b}+\frac {\left (a+\frac {c}{2}+\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) c_2 \cos \left (\left (\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right )}{b}+\frac {\left (a+\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) c_3 \sin \left (\frac {1}{2} \left (c+\sqrt {4 a^2+4 c a+4 b^2+c^2}\right ) t\right )}{b}+\frac {\left (a+\frac {c}{2}+\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) c_1 \sin \left (\left (\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right )}{b},\text {x4}(t)\to -\frac {\left (a+\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) c_3 \cos \left (\frac {1}{2} \left (c+\sqrt {4 a^2+4 c a+4 b^2+c^2}\right ) t\right )}{b}-\frac {\left (a+\frac {c}{2}+\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) c_1 \cos \left (\left (\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right )}{b}+\frac {\left (a+\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) c_4 \sin \left (\frac {1}{2} \left (c+\sqrt {4 a^2+4 c a+4 b^2+c^2}\right ) t\right )}{b}+\frac {\left (a+\frac {c}{2}+\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) c_2 \sin \left (\left (\frac {c}{2}-\frac {1}{2} \sqrt {4 b^2+(2 a+c)^2}\right ) t\right )}{b}\right \}\right \}\] ✓ Maple : cpu = 1.161 (sec), leaf count = 2788
\[\left \{\mathit {x1} \relax (t ) = c_{2}+c_{3} \sin \left (c t \right )+c_{4} \cos \left (c t \right ), \mathit {x2} \relax (t ) = -\cos \left (c t \right ) c_{3}+\sin \left (c t \right ) c_{4}+c_{1}, \mathit {x3} \relax (t ) = \frac {b \left (\cos \left (c t \right ) c_{1} a -\sin \left (c t \right ) c_{2} a -c_{3} \left (a +c \right )\right )}{a \left (a +c \right )}, \mathit {x4} \relax (t ) = \frac {b \left (\cos \left (c t \right ) c_{2} a +\sin \left (c t \right ) c_{1} a +c_{4} \left (a +c \right )\right )}{a \left (a +c \right )}\right \}\]