1.67 problem 67
Internal
problem
ID
[7759]
Book
:
Own
collection
of
miscellaneous
problems
Section
:
section
1.0
Problem
number
:
67
Date
solved
:
Monday, October 21, 2024 at 04:02:52 PM
CAS
classification
:
[[_2nd_order, _missing_x], [_2nd_order, _reducible, _mu_x_y1]]
Solve
\begin{align*} 3 y y^{\prime \prime }+y&=5 \end{align*}
1.67.1 Solved as second order missing x ode
Time used: 498.742 (sec)
This is missing independent variable second order ode. Solved by reduction of order by using
substitution which makes the dependent variable \(y\) an independent variable. Using
\begin{align*} y' &= p \end{align*}
Then
\begin{align*} y'' &= \frac {dp}{dx}\\ &= \frac {dp}{dy}\frac {dy}{dx}\\ &= p \frac {dp}{dy} \end{align*}
Hence the ode becomes
\begin{align*} 3 y p \left (y \right ) \left (\frac {d}{d y}p \left (y \right )\right )+y = 5 \end{align*}
Which is now solved as first order ode for \(p(y)\).
The ode \(p^{\prime } = -\frac {y -5}{3 p y}\) is separable as it can be written as
\begin{align*} p^{\prime }&= -\frac {y -5}{3 p y}\\ &= f(y) g(p) \end{align*}
Where
\begin{align*} f(y) &= -\frac {y -5}{3 y}\\ g(p) &= \frac {1}{p} \end{align*}
Integrating gives
\begin{align*} \int { \frac {1}{g(p)} \,dp} &= \int { f(y) \,dy}\\ \int { p\,dp} &= \int { -\frac {y -5}{3 y} \,dy}\\ \frac {p^{2}}{2}&=-\frac {y}{3}+\ln \left (y^{{5}/{3}}\right )+c_1 \end{align*}
Solving for \(p\) from the above solution(s) gives (after possible removing of solutions that do not
verify)
\begin{align*} p&=-\frac {\sqrt {18 \ln \left (y^{{5}/{3}}\right )+18 c_1 -6 y}}{3}\\ p&=\frac {\sqrt {18 \ln \left (y^{{5}/{3}}\right )+18 c_1 -6 y}}{3} \end{align*}
For solution (1) found earlier, since \(p=y^{\prime }\) then we now have a new first order ode to solve which is
\begin{align*} y^{\prime } = -\frac {\sqrt {18 \ln \left (y^{{5}/{3}}\right )+18 c_1 -6 y}}{3} \end{align*}
Since initial conditions \(\left (x_0,y_0\right ) \) are given, then the result can be written as Since unable to evaluate
the integral, and no initial conditions are given, then the result becomes
\[ \int _{}^{y}-\frac {3}{\sqrt {18 \ln \left (\tau ^{{5}/{3}}\right )+18 c_1 -6 \tau }}d \tau = x +c_2 \]
Singular solutions
are found by solving
\begin{align*} -\frac {\sqrt {18 \ln \left (y^{{5}/{3}}\right )+18 c_1 -6 y}}{3}&= 0 \end{align*}
for \(y\). This is because we had to divide by this in the above step. This gives the following
singular solution(s), which also have to satisfy the given ODE.
\begin{align*} y = \left ({\mathrm e}^{-\frac {5 \operatorname {LambertW}\left (-\frac {243^{{4}/{5}} {\mathrm e}^{-\frac {3 c_1}{5}}}{405}\right )}{3}-c_1}\right )^{{3}/{5}} \end{align*}
For solution (2) found earlier, since \(p=y^{\prime }\) then we now have a new first order ode to solve which is
\begin{align*} y^{\prime } = \frac {\sqrt {18 \ln \left (y^{{5}/{3}}\right )+18 c_1 -6 y}}{3} \end{align*}
Since initial conditions \(\left (x_0,y_0\right ) \) are given, then the result can be written as Since unable to evaluate
the integral, and no initial conditions are given, then the result becomes
\[ \int _{}^{y}\frac {3}{\sqrt {18 \ln \left (\tau ^{{5}/{3}}\right )+18 c_1 -6 \tau }}d \tau = x +c_3 \]
Singular solutions
are found by solving
\begin{align*} \frac {\sqrt {18 \ln \left (y^{{5}/{3}}\right )+18 c_1 -6 y}}{3}&= 0 \end{align*}
for \(y\). This is because we had to divide by this in the above step. This gives the following
singular solution(s), which also have to satisfy the given ODE.
\begin{align*} y = \left ({\mathrm e}^{-\frac {5 \operatorname {LambertW}\left (-\frac {243^{{4}/{5}} {\mathrm e}^{-\frac {3 c_1}{5}}}{405}\right )}{3}-c_1}\right )^{{3}/{5}} \end{align*}
Will add steps showing solving for IC soon.
The solution
\[
\int _{}^{y}-\frac {3}{\sqrt {18 \ln \left (\tau ^{{5}/{3}}\right )+18 c_1 -6 \tau }}d \tau = x +c_2
\]
was found not to satisfy the ode or the IC. Hence it is removed. The solution
\[
y = \left ({\mathrm e}^{-\frac {5 \operatorname {LambertW}\left (-\frac {243^{{4}/{5}} {\mathrm e}^{-\frac {3 c_1}{5}}}{405}\right )}{3}-c_1}\right )^{{3}/{5}}
\]
was found not to satisfy the ode or the IC. Hence it is removed.
1.67.2 Maple step by step solution
1.67.3 Maple dsolve solution
Solving time : 597.045
(sec)
Leaf size : maple_leaf_size
\[ \text {No solution found} \]
1.67.4 Mathematica DSolve solution
Solving time : 0.324
(sec)
Leaf size : 41
DSolve[{3*y[x]*D[y[x],{x,2}]+y[x]==5,{}},
y[x],x,IncludeSingularSolutions->True]
\[
\text {Solve}\left [\int _1^{y(x)}\frac {1}{\sqrt {c_1+\frac {2}{3} (5 \log (K[1])-K[1])}}dK[1]{}^2=(x+c_2){}^2,y(x)\right ]
\]