Solving a system of third order homogenous ODEs I am trying to solve this third order system of homogenous ODEs. 


*

*$x'''=2x+y$

*$y'''=x+2y$
Initial conditions are given as well.
Higher order systems weren't covered in the lectures hence I am a bit lost. Nor can I find any similar questions asked on this site. I have tried to solve the equations in two ways so far. Firstly express $y$  in terms of $x'''$ and $x$ from the first equation differentiate the equation three times giving: $x^{(6)}-2x=y'''$ and plugging this into the second equation finally obtaining:
$x^{(6)}-4x'''+3x=0$
A higher order homogeneous ODE with constant coefficients that I can solve with the substitution $x = e^{\lambda x}$. The solution can then be differentiated three times and plugged into $x'''=2x+y$ to obtain $y$. Using the given initial conditions we can then find the specific solution. However, this proves very tedious, especially solving the 6x6 system of equations. 
Is this approach correct and applicable to other similar systems? Or is there a completely different approach to this type of problem I am not aware of? Could you try to solve the system $\vec{u}''' = A \vec{u}$ by diagonalizing A and look for a solution that way, I've tried this as well but I cannot seem to find a solution. 
Any help is much appreciated. 

Because of the symmetric cyclic/circulant nature of the system matrix, you can decouple the system by Fourier transform. Set $u=x+y$, $v=x-y$ to get the system
$$
u'''=3u\\
v'''=v
$$
which can be both solved with standard methods for linear ODE with constant coefficients.
