I am trying to solve the following system of 2nd order differential equations ($\dot x$ denotes the derivative with respect to $t$, $x$ and $y$ are functions of $t$): $$\begin{cases} \dot{x}=3x-3y+5e^{2t} &(1)\\ \dot{y}=3x-3y &(2) \end{cases}$$
So far, I have done the following:
Differentiating both sides of $(1)$ with respect to $t$, $$\ddot{x}=3\dot x -3\dot y +10e^{2t}$$ Substituting $\dot y$, $$\ddot{x}=3\dot x -9x+9y+10e^{2t}$$ $$\ddot{x}-3\dot x +9x = 9y+10e^{2t}$$ Doing the same for $y$ gives $$\ddot y+3\dot y+9y = 9x+15e^{2t}$$
I don't know where to go further with this. I don't think I can use the method of undetermined coefficients here since the RHS is a function of both $t$ and another variable ($x$ or $y$).
Wolfram Alpha gives a neat (not messy) solution to the problem, so I assume we don't need some crazy methods to solve this.