Internal problem ID [9871]
Internal file name [OUTPUT/8818_Monday_June_06_2022_05_33_17_AM_19505267/index.tex]
Book: Differential Gleichungen, E. Kamke, 3rd ed. Chelsea Pub. NY, 1948
Section: Chapter 4, linear fourth order
Problem number: 1549.
ODE order: 4.
ODE degree: 1.
The type(s) of ODE detected by this program : "higher_order_missing_y"
Maple gives the following as the ode type
[[_high_order, _missing_y]]
\[ \boxed {x y^{\prime \prime \prime \prime }+5 y^{\prime \prime \prime }=24} \] Since \(y\) is missing from the ode then we can use the substitution \(y^{\prime } = v \left (x \right )\) to reduce the order by one. The ODE becomes \begin {align*} x v^{\prime \prime \prime }\left (x \right )+5 v^{\prime \prime }\left (x \right ) = 0 \end {align*}
Since \(v \left (x \right )\) is missing from the ode then we can use the substitution \(v^{\prime }\left (x \right ) = w \left (x \right )\) to reduce the order by one. The ODE becomes \begin {align*} x w^{\prime \prime }\left (x \right )+5 w^{\prime }\left (x \right ) = 0 \end {align*}
Integrating both sides of the ODE w.r.t \(x\) gives \begin {align*} \int \left (x w^{\prime \prime }\left (x \right )+5 w^{\prime }\left (x \right )\right )d x &= 0 \\ w^{\prime }\left (x \right ) x +4 w \left (x \right ) = c_{1} \end {align*}
Which is now solved for \(w \left (x \right )\). In canonical form the ODE is \begin {align*} w' &= F(x,w)\\ &= f( x) g(w)\\ &= \frac {-4 w +c_{1}}{x} \end {align*}
Where \(f(x)=\frac {1}{x}\) and \(g(w)=-4 w +c_{1}\). Integrating both sides gives \begin{align*} \frac {1}{-4 w +c_{1}} \,dw &= \frac {1}{x} \,d x \\ \int { \frac {1}{-4 w +c_{1}} \,dw} &= \int {\frac {1}{x} \,d x} \\ -\frac {\ln \left (w -\frac {c_{1}}{4}\right )}{4}&=\ln \left (x \right )+c_{2} \\ \end{align*} Raising both side to exponential gives \begin {align*} \frac {\sqrt {2}}{\left (4 w -c_{1} \right )^{\frac {1}{4}}} &= {\mathrm e}^{\ln \left (x \right )+c_{2}} \end {align*}
Which simplifies to \begin {align*} \frac {\sqrt {2}}{\left (4 w -c_{1} \right )^{\frac {1}{4}}} &= c_{3} x \end {align*}
But since \(v^{\prime }\left (x \right ) = w \left (x \right )\) then we now need to solve the ode \(v^{\prime }\left (x \right ) = \frac {\left (4+c_{1} c_{3}^{4} {\mathrm e}^{4 c_{2}} x^{4}\right ) {\mathrm e}^{-4 c_{2}}}{4 c_{3}^{4} x^{4}}\). Integrating both sides gives \begin {align*} v \left (x \right ) &= \int { \frac {\left (4+c_{1} c_{3}^{4} {\mathrm e}^{4 c_{2}} x^{4}\right ) {\mathrm e}^{-4 c_{2}}}{4 c_{3}^{4} x^{4}}\,\mathop {\mathrm {d}x}}\\ &= \frac {c_{1} x}{4}-\frac {{\mathrm e}^{-4 c_{2}}}{3 x^{3} c_{3}^{4}}+c_{4} \end {align*}
But since \(y^{\prime } = v \left (x \right )\) then we now need to solve the ode \(y^{\prime } = \frac {c_{1} x}{4}-\frac {{\mathrm e}^{-4 c_{2}}}{3 x^{3} c_{3}^{4}}+c_{4}\). Integrating both sides gives \begin {align*} y &= \int { -\frac {-3 c_{1} c_{3}^{4} x^{4}-12 c_{4} x^{3} c_{3}^{4}+4 \,{\mathrm e}^{-4 c_{2}}}{12 x^{3} c_{3}^{4}}\,\mathop {\mathrm {d}x}}\\ &= \frac {c_{1} x^{2}}{8}+c_{4} x +\frac {{\mathrm e}^{-4 c_{2}}}{6 x^{2} c_{3}^{4}}+c_{5} \end {align*}
This is higher order nonhomogeneous ODE. Let the solution be \[ y = y_h + y_p \] Where \(y_h\) is the solution to the homogeneous ODE And \(y_p\) is a particular solution to the nonhomogeneous ODE. \(y_h\) is the solution to \[ x y^{\prime \prime \prime \prime }+5 y^{\prime \prime \prime } = 0 \] Let the particular solution be \[ y_p = U_1 y_1+U_2 y_2+U_3 y_3+U_4 y_4 \] Where \(y_i\) are the basis solutions found above for the homogeneous solution \(y_h\) and \(U_i(x)\) are functions to be determined as follows \[ U_i = (-1)^{n-i} \int { \frac {F(x) W_i(x) }{a W(x)} \, dx} \] Where \(W(x)\) is the Wronskian and \(W_i(x)\) is the Wronskian that results after deleting the last row and the \(i\)-th column of the determinant and \(n\) is the order of the ODE or equivalently, the number of basis solutions, and \(a\) is the coefficient of the leading derivative in the ODE, and \(F(x)\) is the RHS of the ODE. Therefore, the first step is to find the Wronskian \(W \left (x \right )\). This is given by \begin {equation*} W(x) = \begin {vmatrix} y_1&y_2&y_3&y_4\\ y_1'&y_2'&y_3'&y_4'\\ y_1''&y_2''&y_3''&y_4''\\ y_1'''&y_2'''&y_3'''&y_4'''\\ \end {vmatrix} \end {equation*} Substituting the fundamental set of solutions \(y_i\) found above in the Wronskian gives \begin {align*} W &= \left [\begin {array}{cccc} 1 & x & \frac {1}{x^{2}} & x^{2} \\ 0 & 1 & -\frac {2}{x^{3}} & 2 x \\ 0 & 0 & \frac {6}{x^{4}} & 2 \\ 0 & 0 & -\frac {24}{x^{5}} & 0 \end {array}\right ] \\ |W| &= \frac {48}{x^{5}} \end {align*}
The determinant simplifies to \begin {align*} |W| &= \frac {48}{x^{5}} \end {align*}
Now we determine \(W_i\) for each \(U_i\). \begin {align*} W_1(x) &= \det \,\left [\begin {array}{ccc} x & \frac {1}{x^{2}} & x^{2} \\ 1 & -\frac {2}{x^{3}} & 2 x \\ 0 & \frac {6}{x^{4}} & 2 \end {array}\right ] \\ &= -\frac {12}{x^{2}} \end {align*}
\begin {align*} W_2(x) &= \det \,\left [\begin {array}{ccc} 1 & \frac {1}{x^{2}} & x^{2} \\ 0 & -\frac {2}{x^{3}} & 2 x \\ 0 & \frac {6}{x^{4}} & 2 \end {array}\right ] \\ &= -\frac {16}{x^{3}} \end {align*}
\begin {align*} W_3(x) &= \det \,\left [\begin {array}{ccc} 1 & x & x^{2} \\ 0 & 1 & 2 x \\ 0 & 0 & 2 \end {array}\right ] \\ &= 2 \end {align*}
\begin {align*} W_4(x) &= \det \,\left [\begin {array}{ccc} 1 & x & \frac {1}{x^{2}} \\ 0 & 1 & -\frac {2}{x^{3}} \\ 0 & 0 & \frac {6}{x^{4}} \end {array}\right ] \\ &= \frac {6}{x^{4}} \end {align*}
Now we are ready to evaluate each \(U_i(x)\). \begin {align*} U_1 &= (-1)^{4-1} \int { \frac {F(x) W_1(x) }{a W(x)} \, dx}\\ &= (-1)^{3} \int { \frac { \left (24\right ) \left (-\frac {12}{x^{2}}\right )}{\left (x\right ) \left (\frac {48}{x^{5}}\right )} \, dx} \\ &= - \int { \frac {-\frac {288}{x^{2}}}{\frac {48}{x^{4}}} \, dx}\\ &= - \int {\left (-6 x^{2}\right ) \, dx}\\ &= 2 x^{3} \end {align*}
\begin {align*} U_2 &= (-1)^{4-2} \int { \frac {F(x) W_2(x) }{a W(x)} \, dx}\\ &= (-1)^{2} \int { \frac { \left (24\right ) \left (-\frac {16}{x^{3}}\right )}{\left (x\right ) \left (\frac {48}{x^{5}}\right )} \, dx} \\ &= \int { \frac {-\frac {384}{x^{3}}}{\frac {48}{x^{4}}} \, dx}\\ &= \int {\left (-8 x\right ) \, dx}\\ &= -4 x^{2} \end {align*}
\begin {align*} U_3 &= (-1)^{4-3} \int { \frac {F(x) W_3(x) }{a W(x)} \, dx}\\ &= (-1)^{1} \int { \frac { \left (24\right ) \left (2\right )}{\left (x\right ) \left (\frac {48}{x^{5}}\right )} \, dx} \\ &= - \int { \frac {48}{\frac {48}{x^{4}}} \, dx}\\ &= - \int {\left (x^{4}\right ) \, dx}\\ &= -\frac {x^{5}}{5} \end {align*}
\begin {align*} U_4 &= (-1)^{4-4} \int { \frac {F(x) W_4(x) }{a W(x)} \, dx}\\ &= (-1)^{0} \int { \frac { \left (24\right ) \left (\frac {6}{x^{4}}\right )}{\left (x\right ) \left (\frac {48}{x^{5}}\right )} \, dx} \\ &= \int { \frac {\frac {144}{x^{4}}}{\frac {48}{x^{4}}} \, dx}\\ &= \int {\left (3\right ) \, dx}\\ &= 3 x \end {align*}
Now that all the \(U_i\) functions have been determined, the particular solution is found from \[ y_p = U_1 y_1+U_2 y_2+U_3 y_3+U_4 y_4 \] Hence \begin {equation*} \begin {split} y_p &= \left (2 x^{3}\right ) \\ &+\left (-4 x^{2}\right ) \left (x\right ) \\ &+\left (-\frac {x^{5}}{5}\right ) \left (\frac {1}{x^{2}}\right ) \\ &+\left (3 x\right ) \left (x^{2}\right ) \end {split} \end {equation*} Therefore the particular solution is \[ y_p = \frac {4 x^{3}}{5} \] Therefore the general solution is \begin{align*} y &= y_h + y_p \\ &= \left (y &= \frac {c_{1} x^{2}}{8}+c_{4} x +\frac {{\mathrm e}^{-4 c_{2}}}{6 x^{2} c_{3}^{4}}+c_{5}\right ) + \left (\frac {4 x^{3}}{5}\right ) \\ \end{align*}
The solution(s) found are the following \begin{align*} \tag{1} y &= \frac {c_{1} x^{2}}{8}+c_{4} x +\frac {{\mathrm e}^{-4 c_{2}}}{6 x^{2} c_{3}^{4}}+c_{5} +\frac {4 x^{3}}{5} \\ \end{align*}
Verification of solutions
\[ y = \frac {c_{1} x^{2}}{8}+c_{4} x +\frac {{\mathrm e}^{-4 c_{2}}}{6 x^{2} c_{3}^{4}}+c_{5} +\frac {4 x^{3}}{5} \] Verified OK.
\[ \begin {array}{lll} & {} & \textrm {Let's solve}\hspace {3pt} \\ {} & {} & x \left (\frac {d}{d x}\frac {d}{d x}y^{\prime \prime }\right )+5 \frac {d}{d x}y^{\prime \prime }=24 \\ \bullet & {} & \textrm {Highest derivative means the order of the ODE is}\hspace {3pt} 4 \\ {} & {} & \frac {d}{d x}\frac {d}{d x}y^{\prime \prime } \end {array} \]
Maple trace
`Methods for high order ODEs: --- Trying classification methods --- trying a quadrature trying high order exact linear fully integrable -> Calling odsolve with the ODE`, diff(_b(_a), _a) = -(5*_b(_a)-24)/_a, _b(_a)` *** Sublevel 2 *** Methods for first order ODEs: --- Trying classification methods --- trying a quadrature trying 1st order linear <- 1st order linear successful <- high order exact linear fully integrable successful`
✓ Solution by Maple
Time used: 0.0 (sec). Leaf size: 26
dsolve(x*diff(diff(diff(diff(y(x),x),x),x),x)+5*diff(diff(diff(y(x),x),x),x)-24=0,y(x), singsol=all)
\[ y \left (x \right ) = \frac {4 x^{3}}{5}-\frac {c_{1}}{24 x^{2}}+\frac {c_{2} x^{2}}{2}+c_{3} x +c_{4} \]
✓ Solution by Mathematica
Time used: 0.051 (sec). Leaf size: 34
DSolve[-24 + 5*Derivative[3][y][x] + x*Derivative[4][y][x] == 0,y[x],x,IncludeSingularSolutions -> True]
\[ y(x)\to \frac {4 x^3}{5}+c_4 x^2-\frac {c_1}{24 x^2}+c_3 x+c_2 \]