How to check a particular solution of System of Linear ODEs?

In summary, to check if a given solution works for a system of inhomogeneous linear ODEs, one can use symbolic differentiation in Mathematica to verify that the proposed solution satisfies the ODEs.
  • #1
Hall
351
87
If I have been given a system of inhomogeneous linear ODEs,
$$
\vec{x'} =
\begin{bmatrix}
4 & -1 \\
5 & -2 \\
\end{bmatrix}
\vec{x}
+
\begin{bmatrix}
18e^{2t} \\
30e^{2t}\\
\end{bmatrix}
$$

I have found its particular solution to be:
$$
1/4
\begin{bmatrix}
-31e^{2t} - 25e^{6t} \\
85e^{2t} - 25e^{6t} \\
\end{bmatrix}
$$

But this answer doesn't match with the answer given in the book. Can someone tell me how to check if this solutions works by writing some code in Mathematica? I know, I can use DSolve for solving them, but I'm asking a reverse of that.

Please guide me step by step, I'm new to Mathematica and I don't have any background in programming.
 
Last edited:
Physics news on Phys.org
  • #2
Hall said:
If I have been given a system of inhomogeneous linear ODEs,
$$
\vec{x'} =
\begin{bmatrix}
4 & -1 \\
5 & -2 \\
\end{bmatrix}
\vec{x}
+
\begin{bmatrix}
18e^{2t} \\
30e^{2t}\\
\end{bmatrix}
$$

I have found its particular solution to be:
$$
1/4
\begin{bmatrix}
-31e^{2t} - 25e^{6t} \\
85e^{2t} - 25e^{6t} \\
\end{bmatrix}
$$

But this answer doesn't match with the answer given in the book.

I look at your solution, and I must ask myself: Where does [itex]e^{6t}[/itex] come from? [itex]e^{2t}[/itex] is an eigenfunction of the derivative operator: [itex](e^{2t})' = 2e^{2t}[/itex]. So I would expect the particular function to be [itex]ae^{2t}[/itex] for some constant vector [itex]a[/itex], which can be determined by substituting this into the ODE. That your answer is not of this form, or of this form plus a complementary function (6 is not an eigenvalue of the matrix, so a multiple of [itex]e^{6t}[/itex] is not a complementary function), leads me to suspect that you have made an error, but since you haven't shown your working I can't tell you what it is.

Can someone tell me how to check if this solutions works by writing some code in Mathematica? I know, I can use DSolve for solving them, but I'm asking a reverse of that.

Please guide me step by step, I'm new to Mathematica and I don't have any background in programming.

The answer to this question is "symbolic differentiation". Define your proposed solution as a function, and check to see that
Code:
f'[t] - {{4, -1},{5, -2}} . f[t] - {18*exp[2*t], 30*exp[2*t]}
is zero. Have a look at the examples at https://reference.wolfram.com/language/ref/Derivative.html.
 
Last edited:
  • Like
Likes Hall

Similar threads

  • Differential Equations
Replies
2
Views
1K
  • Precalculus Mathematics Homework Help
Replies
19
Views
1K
  • Calculus and Beyond Homework Help
Replies
3
Views
589
  • Precalculus Mathematics Homework Help
Replies
21
Views
1K
  • Differential Equations
Replies
2
Views
2K
  • Calculus and Beyond Homework Help
Replies
2
Views
557
  • Engineering and Comp Sci Homework Help
Replies
2
Views
1K
  • MATLAB, Maple, Mathematica, LaTeX
Replies
0
Views
1K
  • Calculus and Beyond Homework Help
Replies
8
Views
1K
  • Calculus and Beyond Homework Help
Replies
4
Views
962
Back
Top