In Problems 25 – 28, use the elimination method to find a general solution for the given system of three equations in the three unknown functions x(t), y(t), z(t).

$$\begin{gathered} \mathbf{x}\mathbf{\text{'}}\mathbf{=}\mathbf{3}\mathbf{x}\mathbf{+}\mathbf{y}\mathbf{-}\mathbf{z}\mathbf{,} \\ \mathbf{y}\mathbf{\text{'}}\mathbf{=}\mathbf{x}\mathbf{+}\mathbf{2}\mathbf{y}\mathbf{-}\mathbf{z}\mathbf{,} \\ \mathbf{z}\mathbf{\text{'}}\mathbf{=}\mathbf{3}\mathbf{x}\mathbf{+}\mathbf{3}\mathbf{y}\mathbf{-}\mathbf{z} \end{gathered}$$

Elimination Procedure for 2 × 2 Systems:

To find a general solution for the system

$$\begin{gathered} {\mathbf{L}}_{\mathbf{1}}\left[\mathbf{x}\right]\mathbf{+}{\mathbf{L}}_{\mathbf{2}}\left[\mathbf{y}\right]\mathbf{=}{\mathbf{f}}_{\mathbf{1}}\mathbf{,} \\ {\mathbf{L}}_{\mathbf{3}}\left[\mathbf{x}\right]\mathbf{+}{\mathbf{L}}_{\mathbf{4}}\left[\mathbf{y}\right]\mathbf{=}{\mathbf{f}}_{\mathbf{2}}\mathbf{,} \end{gathered}$$

Where${\mathbf{L}}_{\mathbf{1}}\mathbf{,}{\mathbf{L}}_{\mathbf{2}}\mathbf{,}{\mathbf{L}}_{\mathbf{3}}\mathbf{,}$ and L₄are polynomials in $D=\frac{d}{dt}$:

1. Make sure that the system is written in operator form.
2. Eliminate one of the variables, say, y, and solve the resulting equation for x(t). If the system is degenerating, stop! A separate analysis is required to determine whether or not there are solutions.
3. (Shortcut) If possible, use the system to derive an equation that involves y(t) but not its derivatives. [Otherwise, go to step (d).] Substitute the found expression for x(t) into this equation to get a formula for y(t). The expressions for x(t), and y(t) give the desired general solution.
4. Eliminate x from the system and solve for y(t). [Solving for y(t) gives more constants----in fact, twice as many as needed.]
5. Remove the extra constants by substituting the expressions for x(t) and y(t) into one or both of the equations in the system. Write the expressions for x(t) and y(t) in terms of the remaining constants.

Given that,$x\text{'}=3x+y-z ......\mathbf{\left(}\mathbf{1}\mathbf{\right)}$

$$\begin{gathered} y\text{'}=x+2y-z ......\mathbf{\left(}\mathbf{2}\mathbf{\right)} \\ z\text{'}=3x+3y-z ......\mathbf{\left(}\mathbf{3}\mathbf{\right)} \end{gathered}$$

Let us rewrite the given system of equations into operator form.

$$\begin{gathered} \left(D-3\right)\left[x\right]-\left[y\right]+\left[z\right]=0 ......\mathbf{\left(}\mathbf{4}\mathbf{\right)} \\ \left[x\right]-\left(D-2\right)\left[y\right]-\left[z\right]=0 ......\mathbf{\left(}\mathbf{5}\mathbf{\right)} \\ 3\left[x\right]+3\left[y\right]-\left(D+1\right)\left[z\right]=0 ......\mathbf{\left(}\mathbf{6}\mathbf{\right)} \end{gathered}$$

Add the equation (4) and (5) together to eliminate z(t),

$$\begin{gathered} \left(D-3\right)\left[x\right]-\left[y\right]+\left[z\right]+\left[x\right]-\left(D-2\right)\left[y\right]-\left[z\right]=0 \\ \left(D-3+1\right)\left[x\right]-\left(D-2+1\right)\left[y\right]=0 \\ \left(D-2\right)\left[x\right]-\left(D-1\right)\left[y\right]=0 \\ \left(D-2\right)\left[x\right]-\left(D-1\right)\left[y\right]=0 ......\mathbf{\left(}\mathbf{7}\mathbf{\right)} \end{gathered}$$

Multiply D+1 on equation (5). Then, subtract the equation with equation (6).

$$\begin{gathered} \left(D+1\right)\left[x\right]-\left(D+1\right)\left(D-2\right)\left[y\right]-\left(D+1\right)\left[z\right]-3\left[x\right]-3\left[y\right]+\left(D+1\right)\left[z\right]=0 \\ \left(D+1-3\right)\left[x\right]-\left(D^2-D-2+3\right)\left[y\right]=0 \\ \left(D-2\right)\left[x\right]-\left(D^2-D+1\right)\left[y\right]=0 \\ \left(D-2\right)\left[x\right]-\left(D^2-D+1\right)\left[y\right]=0 ......\mathbf{\left(}\mathbf{8}\mathbf{\right)} \end{gathered}$$

Subtract equations (7) and (8) to get

$$\begin{gathered} \left(D-2\right)\left[x\right]-\left(D-1\right)\left[y\right]-\left(D-2\right)\left[x\right]+\left(D^2-D+1\right)\left[y\right]=0 \\ \left(D^2-D+1-D+1\right)\left[y\right]=0 \\ \left(D^2-2D+2\right)\left[y\right]=0 \end{gathered}$$

Since the auxiliary equation to the corresponding homogeneous equation is $r^2-2r+2=0$. The roots are$r=1+i$and $r=1-i$.

Then, the general solution of y is

$y\left(t\right)=A{e}^t\sin t+B{e}^t\cos t ......\mathbf{\left(}\mathbf{9}\mathbf{\right)}$

Now derivate the equation (9)

$$\begin{gathered} D\left[y\left(t\right)\right]=A{e}^t\sin t+A{e}^t\cos t+B{e}^t\cos t-B{e}^t\sin t \\ =(A-B)e^t\sin t+(A+B)e^t\cos t \end{gathered}$$

Substitute the derivation in equation (7).

$$\begin{gathered} \left(D-2\right)\left[x\right]-\left(D-1\right)\left[y\right]=0 \\ \left(D-2\right)\left[x\right]=\left(D-1\right)\left[A{e}^t\sin t+B{e}^t\cos t\right] \\ =A{e}^t\sin t+A{e}^t\cos t+B{e}^t\cos t-B{e}^t\sin t-A{e}^t\sin t-B{e}^t\cos t \\ =A{e}^t\cos t-B{e}^t\sin t \end{gathered}$$

$\left(D-2\right)\left[x\right]=A{e}^t\cos t-B{e}^t\sin t ......\mathbf{\left(}\mathbf{10}\mathbf{\right)}$

So, the complementary solution of the differential equation is$x_h\left(t\right)=C{e}^{2t} ......\mathbf{\left(}\mathbf{11}\mathbf{\right)}$

Let us assume that function:$x_p\left(t\right)=D{e}^t\sin t+E{e}^t\cos t ......\mathbf{\left(}\mathbf{12}\mathbf{\right)}$

Find the derivation of equation (12).

$D\left[x_p\left(t\right)\right]=D{e}^t\sin t+D{e}^t\cos t+E{e}^t\cos t-E{e}^t\sin t$

Use the derivation in equations (10) to get,

$$\begin{gathered} \left(D-2\right)\left[x\right]=A{e}^t\cos t-B{e}^t\sin t \\ \left(D-2\right)\left[D{e}^t\sin t+E{e}^t\cos t\right]=A{e}^t\cos t-B{e}^t\sin t \\ D{e}^t\sin t+D{e}^t\cos t-E{e}^t\sin t+E{e}^t\cos t-2D{e}^t\sin t-2E{e}^t\cos t=A{e}^t\cos t-B{e}^t\sin t \\ D{e}^t\cos t-E{e}^t\sin t-D{e}^t\sin t-E{e}^t\cos t=A{e}^t\cos t-B{e}^t\sin t \\ \left(D-E\right)e^t\cos t-\left(D+E\right)e^t\sin t=A{e}^t\cos t-B{e}^t\sin t \end{gathered}$$

Now, equalize the like terms.

$$\begin{gathered} D-E=A \\ D+E=B \end{gathered}$$

Then,

$$\begin{gathered} 2D=A+B \\ D=\frac{A+B}{2} \end{gathered}$$

$$\begin{gathered} \frac{A+B}{2}+E=B \\ E=B-\frac{A+B}{2} \\ =\frac{B-A}{2} \end{gathered}$$

So, $x_p\left(t\right)=\frac{1}{2}\left(A+B\right)e^t\sin t+\frac{1}{2}\left(B-A\right)e^t\cos t$.

Then,$x\left(t\right)=C{e}^{2t}+\frac{1}{2}\left(A+B\right)e^t\sin t+\frac{1}{2}\left(B-A\right)e^t\cos t ......\mathbf{\left(}\mathbf{13}\mathbf{\right)}$

Now substitute the equation (9) and (13) in equation (4)

role="math" $$\begin{gathered} \left(D-3\right)\left[x\right]-\left[y\right]+\left[z\right]=0 \\ z=\left[y\right]-\left(D-3\right)\left[x\right] \\ =A{e}^t\sin t+B{e}^t\cos t-\left(D-3\right)\left[C{e}^{2t}+\frac{1}{2}\left(A+B\right)e^t\sin t+\frac{1}{2}\left(B-A\right)e^t\cos t\right] \\ =A{e}^t\sin t+B{e}^t\cos t-2C{e}^{2t}+\frac{1}{2}\left(A+B\right)e^t\sin t+\frac{1}{2}\left(A+B\right)e^t\cos t+\frac{1}{2}\left(B-A\right)e^t\cos t-\frac{1}{2}\left(B-A\right)e^t\sin t \\ +3C{e}^{2t}+\frac{3}{2}\left(A+B\right)e^t\sin t+\frac{3}{2}\left(B-A\right)e^t\cos t \\ z\left(t\right)=\frac{3}{2}\left(A+B\right)e^t\sin t+\frac{3}{2}\left(A-B\right)e^t\cos t+C{e}^{2t} \end{gathered}$$

So, the solution is founded.