Solving a nonhomogeneous linear system \(\mathbf{X}^{\prime}=\mathbf{A
X}+\mathbf{F}(t)\) by variation of parameters when \(\mathbf{A}\) is a \(3 \times
3\) (or larger) matrix is almost an impossible task to do by hand. Consider the
system
$$\mathbf{X}^{\prime}=\left(\begin{array}{rrrr}
2 & -2 & 2 & 1 \\
-1 & 3 & 0 & 3 \\
0 & 0 & 4 & -2 \\
0 & 0 & 2 & -1
\end{array}\right) \mathbf{X}+\left(\begin{array}{c}
t e^{t} \\
e^{-t} \\
e^{2 t} \\
1
\end{array}\right)$$
(a) Use a CAS or linear algebra software to find the eigenvalues and
eigenvectors of the coefficient matrix.
(b) Form a fundamental matrix \(\Phi(t)\) and use the computer to find
\(\Phi^{-1}(t)\)
(c) Use the computer to carry out the computations of:
\(\Phi^{-1}(t) \mathbf{F}(t), \quad \int \Phi^{-1}(t) \mathbf{F}(t) d t, \quad
\Phi(t) \int \Phi^{-1}(t) \mathbf{F}(t) d t\)
\(\Phi(t) \mathbf{C},\) and \(\Phi(t) \mathbf{C}+\int \Phi^{-1}(t) \mathbf{F}(t)
d t,\) where \(\mathbf{C}\) is a column
matrix of constants \(c_{1}, c_{2}, c_{3},\) and \(c_{4}\)
(d) Rewrite the computer output for the general solution of the system in the
form \(\mathbf{X}=\mathbf{X}_{c}+\mathbf{X}_{p},\) where \(\mathbf{X}_{c}=c_{1}
\mathbf{X}_{1}+c_{2} \mathbf{X}_{2}+c_{3} \mathbf{X}_{3}+c_{4} \mathbf{X}_{4}\)