Answer on Question #64773 – Math – Differential Equations
Question
Solve ( D 3 − 7 D + 12 ) y = x 3 − 5 x 2 (D3 - 7D + 12)y = x3 - 5x2 ( D 3 − 7 D + 12 ) y = x 3 − 5 x 2
Solution
To solve a nonhomogeneous linear differential equation (DE)
( D 3 − 7 D + 12 ) y = x 3 − 5 x 2 , (D^3 - 7D + 12)y = x^3 - 5x^2, ( D 3 − 7 D + 12 ) y = x 3 − 5 x 2 ,
that is,
y ′ ′ ′ − 7 y ′ + 12 y = x 3 − 5 x 2 y''' - 7y' + 12y = x^3 - 5x^2 y ′′′ − 7 y ′ + 12 y = x 3 − 5 x 2
we must
1) find the complementary function y c y_c y c that is the general solution of the associated homogeneous DE
y ′ ′ ′ − 7 y ′ + 12 y = 0 ; y''' - 7y' + 12y = 0; y ′′′ − 7 y ′ + 12 y = 0 ;
2) find any particular solution y p y_p y p of the nonhomogeneous equation
y ′ ′ ′ − 7 y ′ + 12 y = x 3 − 5 x 2 . y''' - 7y' + 12y = x^3 - 5x^2. y ′′′ − 7 y ′ + 12 y = x 3 − 5 x 2 .
The general solution of the equation (1) is y = y c + y p y = y_c + y_p y = y c + y p .
1) To solve the associated homogeneous DE
y ′ ′ ′ − 7 y ′ + 12 y = 0 y''' - 7y' + 12y = 0 y ′′′ − 7 y ′ + 12 y = 0
we substitute into equation the solution as exponential function y = e m x y = e^{mx} y = e m x . We get
m 3 e m x − 7 m e m x + 5 e m x = 0 m^3 e^{mx} - 7m e^{mx} + 5e^{mx} = 0 m 3 e m x − 7 m e m x + 5 e m x = 0
or
( m 3 − 7 m + 12 ) e m x = 0. (m^3 - 7m + 12)e^{mx} = 0. ( m 3 − 7 m + 12 ) e m x = 0.
Since e m x e^{mx} e m x is never zero for real values of x x x , the last equation is satisfied only when m m m is a solution or root of the third-degree polynomial equation
m 3 − 7 m + 12 = 0. m^3 - 7m + 12 = 0. m 3 − 7 m + 12 = 0.
This equation is called the auxiliary equation of the differential equation (2).
Analytical solution of this equation is a difficult task. Approximate numerical solution is m 1 ≈ 1.63 − 1.00 i m_1 \approx 1.63 - 1.00i m 1 ≈ 1.63 − 1.00 i , m 2 ≈ 1.63 + 1.00 i m_2 \approx 1.63 + 1.00i m 2 ≈ 1.63 + 1.00 i , m 3 ≈ − 3.27 m_3 \approx -3.27 m 3 ≈ − 3.27 .
Thus we have three different roots of the auxiliary equation and the general solution of the associated homogeneous DE is
y c = C 1 e ( 1.63 − 1.00 i ) x + C 2 e ( 1.63 + 1.00 i ) x + C 3 e − 3.27 x . y_c = C_1 e^{(1.63 - 1.00i)x} + C_2 e^{(1.63 + 1.00i)x} + C_3 e^{-3.27x}. y c = C 1 e ( 1.63 − 1.00 i ) x + C 2 e ( 1.63 + 1.00 i ) x + C 3 e − 3.27 x .
Now we use Euler's formula e i β = cos β + i sin β e^{i\beta} = \cos\beta + i\sin\beta e i β = cos β + i sin β . Then we get
y c = c 1 e 1.63 x cos x + c 2 e 1.63 x sin x + c 3 e − 3.27 x y_c = c_1 e^{1.63x} \cos x + c_2 e^{1.63x} \sin x + c_3 e^{-3.27x} y c = c 1 e 1.63 x cos x + c 2 e 1.63 x sin x + c 3 e − 3.27 x
where c 1 , c 2 , c 3 c_1, c_2, c_3 c 1 , c 2 , c 3 are constants.
2) The particular solution of the nonhomogeneous equation will be found in the following form:
y p = A x 3 + B x 2 + C x + D , y_p = A x^3 + B x^2 + C x + D, y p = A x 3 + B x 2 + C x + D ,
where A , B , C , D A, B, C, D A , B , C , D are the undetermined coefficients.
Substituting this solution into the equation (1) and equating the coefficients of corresponding terms, we find A , B , C , D A, B, C, D A , B , C , D . Find derivatives y p ′ , y p ′ ′ , y p ′ ′ ′ y'_p, y''_p, y'''_p y p ′ , y p ′′ , y p ′′′
y p ′ = 3 A x 2 + 2 B x + C y'_p = 3A x^2 + 2B x + C y p ′ = 3 A x 2 + 2 B x + C y p ′ ′ = 6 A x + 2 B y''_p = 6A x + 2B y p ′′ = 6 A x + 2 B y p ′ ′ ′ = 6 A y'''_p = 6A y p ′′′ = 6 A
Substituting the derivatives into the equation (1)
6 A − 7 ( 3 A x 2 + 2 B x + C ) + 12 ( A x 3 + B x 2 + C x + D ) = x 3 − 5 x 2 6A - 7(3Ax^2 + 2Bx + C) + 12(Ax^3 + Bx^2 + Cx + D) = x^3 - 5x^2 6 A − 7 ( 3 A x 2 + 2 B x + C ) + 12 ( A x 3 + B x 2 + C x + D ) = x 3 − 5 x 2
Equating the coefficients of corresponding terms, we find
{ 12 A = 1 , − 21 A + 12 B = − 5 − 14 B + 12 C = 0 6 A − 7 C + 12 D = 0 \begin{cases}
12A = 1, \\
-21A + 12B = -5 \\
-14B + 12C = 0 \\
6A - 7C + 12D = 0
\end{cases} ⎩ ⎨ ⎧ 12 A = 1 , − 21 A + 12 B = − 5 − 14 B + 12 C = 0 6 A − 7 C + 12 D = 0 { A = 1 12 ≈ 0.083 , − 21 ⋅ 0.083 + 12 B = − 5 − 14 B + 12 C = 0 6 ⋅ 0.083 − 7 C + 12 D = 0 \begin{cases}
A = \frac{1}{12} \approx 0.083, \\
-21 \cdot 0.083 + 12B = -5 \\
-14B + 12C = 0 \\
6 \cdot 0.083 - 7C + 12D = 0
\end{cases} ⎩ ⎨ ⎧ A = 12 1 ≈ 0.083 , − 21 ⋅ 0.083 + 12 B = − 5 − 14 B + 12 C = 0 6 ⋅ 0.083 − 7 C + 12 D = 0 { 12 B = − 5 + 1.743 = − 3 , 257 − 14 B + 12 C = 0 0.5 − 7 C + 12 D = 0 \begin{cases}
12B = -5 + 1.743 = -3,257 \\
-14B + 12C = 0 \\
0.5 - 7C + 12D = 0
\end{cases} ⎩ ⎨ ⎧ 12 B = − 5 + 1.743 = − 3 , 257 − 14 B + 12 C = 0 0.5 − 7 C + 12 D = 0 { B = − 39 144 ≈ − 0.271 14 ⋅ 0.271 + 12 C = 0 1 2 − 7 C + 12 D = 0 \begin{cases}
B = -\frac{39}{144} \approx -0.271 \\
14 \cdot 0.271 + 12C = 0 \\
\frac{1}{2} - 7C + 12D = 0
\end{cases} ⎩ ⎨ ⎧ B = − 144 39 ≈ − 0.271 14 ⋅ 0.271 + 12 C = 0 2 1 − 7 C + 12 D = 0 { C ≈ − 0.316 1 2 + 7 ⋅ 0.316 + 12 D = 0 \begin{cases}
C \approx -0.316 \\
\frac{1}{2} + 7 \cdot 0.316 + 12D = 0
\end{cases} { C ≈ − 0.316 2 1 + 7 ⋅ 0.316 + 12 D = 0 D ≈ − 0.226 D \approx -0.226 D ≈ − 0.226
So we have A ≈ 0.083 A \approx 0.083 A ≈ 0.083 , B ≈ − 0.271 B \approx -0.271 B ≈ − 0.271 , C ≈ − 0.316 C \approx -0.316 C ≈ − 0.316 , D ≈ − 0.226 D \approx -0.226 D ≈ − 0.226 and finally we get the general solution of differential equation (1):
y = c 1 e 1.63 x cos x + c 2 e 1.63 x sin x + c 3 e − 3.27 x + 0.083 x 3 − 0.271 x 2 − 0.316 x − 0.226. y = c_1 e^{1.63x} \cos x + c_2 e^{1.63x} \sin x + c_3 e^{-3.27x} + 0.083x^3 - 0.271x^2 - 0.316x - 0.226. y = c 1 e 1.63 x cos x + c 2 e 1.63 x sin x + c 3 e − 3.27 x + 0.083 x 3 − 0.271 x 2 − 0.316 x − 0.226.
Answer:
y = c 1 e 1.63 x cos x + c 2 e 1.63 x sin x + c 3 e − 3.27 x + 0.083 x 3 − 0.271 x 2 − 0.316 x − 0.226. y = c_1 e^{1.63x} \cos x + c_2 e^{1.63x} \sin x + c_3 e^{-3.27x} + 0.083x^3 - 0.271x^2 - 0.316x - 0.226. y = c 1 e 1.63 x cos x + c 2 e 1.63 x sin x + c 3 e − 3.27 x + 0.083 x 3 − 0.271 x 2 − 0.316 x − 0.226.
Answer provided by https://www.AssignmentExpert.com