Output:\begin{align*}u_{1}\left(t\right) &= a^{2} u_{1}\left(0\right) + a b u_{2}\left(0\right) + c a u_{3}\left(0\right) - {\left(a^{2} u_{1}\left(0\right) + a b u_{2}\left(0\right) + c a u_{3}\left(0\right) - u_{1}\left(0\right)\right)} \cos\left(t\right)+{\left(c u_{2}\left(0\right)-b u_{3}\left(0\right)\right)} \sin\left(t\right)\\
u_{2}\left(t\right) &= a b u_{1}\left(0\right) + b^{2} u_{2}\left(0\right) + c b u_{3}\left(0\right) - {\left(a b u_{1}\left(0\right) + b^{2} u_{2}\left(0\right) + c b u_{3}\left(0\right) - u_{2}\left(0\right)\right)} \cos\left(t\right) + {\left(a u_{3}\left(0\right) - c u_{1}\left(0\right)\right)} \sin\left(t\right)\\
u_{3}\left(t\right) &= ca u_{1}\left(0\right) + bc u_{2}\left(0\right)+c^{2} u_{3}\left(0\right) - \left(ca u_{1}\left(0\right) + bc u_{2}\left(0\right)+c^{2} u_{3}\left(0\right)-u_{3}\left(0\right)\right) \cos\left(t\right) + {\left(b u_{1}\left(0\right) - a u_{2}\left(0\right)\right)} \sin\left(t\right)\end{align*}