Internal problem ID [12268]
Internal file name [OUTPUT/10921_Thursday_September_28_2023_01_08_53_AM_70360323/index.tex
]
Book: APPLIED DIFFERENTIAL EQUATIONS The Primary Course by Vladimir A.
Dobrushkin. CRC Press 2015
Section: Chapter 4, Second and Higher Order Linear Differential Equations. Problems page
221
Problem number: Problem 19(a).
ODE order: 2.
ODE degree: 1.
The type(s) of ODE detected by this program : "second_order_integrable_as_is"
Maple gives the following as the ode type
[[_2nd_order, _exact, _nonlinear], [_2nd_order, _reducible, _mu_x_y1], [_2nd_order, _reducible, _mu_y_y1], [_2nd_order, _reducible, _mu_xy]]
Unable to solve or complete the solution.
\[ \boxed {x y^{\prime \prime }+\left (6 x y^{2}+1\right ) y^{\prime }+2 y^{3}=-1} \]
Integrating both sides of the ODE w.r.t \(x\) gives \begin {align*} \int \left (x y^{\prime \prime }+\left (6 x y^{2}+1\right ) y^{\prime }+2 y^{3}\right )d x &= \int \left (-1\right )d x\\ 2 x y^{3}+y^{\prime } x = -x + c_{1} \end {align*}
Which is now solved for \(y\). This is Abel first kind ODE, it has the form \[ y^{\prime }= f_0(x)+f_1(x) y +f_2(x)y^{2}+f_3(x)y^{3} \] Comparing the above to given ODE which is \begin {align*} y^{\prime }&=-2 y^{3}-\frac {x -c_{1}}{x}\tag {1} \end {align*}
Therefore \begin {align*} f_0(x) &= -1+\frac {c_{1}}{x}\\ f_1(x) &= 0\\ f_2(x) &= 0\\ f_3(x) &= -2 \end {align*}
Since \(f_2(x)=0\) then we check the Abel invariant to see if it depends on \(x\) or not. The Abel invariant is given by \begin {align*} -\frac {f_{1}^{3}}{f_{0}^{2} f_{3}} \end {align*}
Which when evaluating gives \begin {align*} \frac {c_{1}^{3}}{54 x^{6} \left (-1+\frac {c_{1}}{x}\right )^{5}} \end {align*}
Since the Abel invariant depends on \(x\) then unable to solve this ode at this time.
Writing the ode as \[ x y^{\prime \prime }+\left (6 x y^{2}+1\right ) y^{\prime }+2 y^{3} = -1 \] Integrating both sides of the ODE w.r.t \(x\) gives \begin {align*} \int \left (x y^{\prime \prime }+\left (6 x y^{2}+1\right ) y^{\prime }+2 y^{3}\right )d x &= \int \left (-1\right )d x\\ y-\frac {y^{\prime \prime \prime } x^{3}}{6}+\frac {y^{\prime \prime } x^{2}}{2}+2 x y^{3} = -x +c_{1} \end {align*}
Which is now solved for \(y\).
Maple trace
`Methods for second order ODEs: --- Trying classification methods --- trying 2nd order Liouville trying 2nd order WeierstrassP trying 2nd order JacobiSN differential order: 2; trying a linearization to 3rd order trying 2nd order ODE linearizable_by_differentiation trying 2nd order, 2 integrating factors of the form mu(x,y) trying differential order: 2; missing variables -> trying 2nd order, dynamical_symmetries, fully reducible to Abel through one integrating factor of the form G(x,y)/(1+H(x,y)*y)^2 --- trying a change of variables {x -> y(x), y(x) -> x} and re-entering methods for dynamical symmetries --- -> trying 2nd order, dynamical_symmetries, fully reducible to Abel through one integrating factor of the form G(x,y)/(1+H(x,y)*y) trying 2nd order, integrating factors of the form mu(x,y)/(y)^n, only the singular cases trying symmetries linear in x and y(x) trying differential order: 2; exact nonlinear -> Calling odsolve with the ODE`, 2*_b(_a)^3*_a+(diff(_b(_a), _a))*_a+_a+c__1 = 0, _b(_a)` *** Sublevel 2 *** Methods for first order ODEs: --- Trying classification methods --- trying a quadrature trying 1st order linear trying Bernoulli trying separable trying inverse linear trying homogeneous types: trying Chini differential order: 1; looking for linear symmetries trying exact trying Abel Looking for potential symmetries Looking for potential symmetries Looking for potential symmetries trying inverse_Riccati trying an equivalence to an Abel ODE differential order: 1; trying a linearization to 2nd order --- trying a change of variables {x -> y(x), y(x) -> x} differential order: 1; trying a linearization to 2nd order trying 1st order ODE linearizable_by_differentiation --- Trying Lie symmetry methods, 1st order --- `, `-> Computing symmetries using: way = 3 `, `-> Computing symmetries using: way = 4 `, `-> Computing symmetries using: way = 2 trying symmetry patterns for 1st order ODEs -> trying a symmetry pattern of the form [F(x)*G(y), 0] -> trying a symmetry pattern of the form [0, F(x)*G(y)] -> trying symmetry patterns of the forms [F(x),G(y)] and [G(y),F(x)] -> trying a symmetry pattern of the form [F(x),G(x)] -> trying a symmetry pattern of the form [F(y),G(y)] -> trying a symmetry pattern of the form [F(x)+G(y), 0] -> trying a symmetry pattern of the form [0, F(x)+G(y)] -> trying a symmetry pattern of the form [F(x),G(x)*y+H(x)] -> trying a symmetry pattern of conformal type trying 2nd order, integrating factor of the form mu(x,y) -> Calling odsolve with the ODE`, 2*_b(_a)^3*_a+(diff(_b(_a), _a))*_a+_a+c__1 = 0, _b(_a)` *** Sublevel 2 *** Methods for first order ODEs: --- Trying classification methods --- trying a quadrature trying 1st order linear trying Bernoulli trying separable trying inverse linear trying homogeneous types: trying Chini differential order: 1; looking for linear symmetries trying exact trying Abel Looking for potential symmetries Looking for potential symmetries Looking for potential symmetries trying inverse_Riccati trying an equivalence to an Abel ODE differential order: 1; trying a linearization to 2nd order --- trying a change of variables {x -> y(x), y(x) -> x} differential order: 1; trying a linearization to 2nd order trying 1st order ODE linearizable_by_differentiation --- Trying Lie symmetry methods, 1st order --- `, `-> Computing symmetries using: way = 3 `, `-> Computing symmetries using: way = 4 `, `-> Computing symmetries using: way = 2 trying symmetry patterns for 1st order ODEs -> trying a symmetry pattern of the form [F(x)*G(y), 0] -> trying a symmetry pattern of the form [0, F(x)*G(y)] -> trying symmetry patterns of the forms [F(x),G(y)] and [G(y),F(x)] -> trying a symmetry pattern of the form [F(x),G(x)] -> trying a symmetry pattern of the form [F(y),G(y)] -> trying a symmetry pattern of the form [F(x)+G(y), 0] -> trying a symmetry pattern of the form [0, F(x)+G(y)] -> trying a symmetry pattern of the form [F(x),G(x)*y+H(x)] -> trying a symmetry pattern of conformal type differential order: 2; looking for linear symmetries -> trying 2nd order, the S-function method *** Sublevel 2 *** Methods for first order ODEs: --- Trying classification methods --- trying a quadrature <- quadrature successful *** Sublevel 2 *** Methods for first order ODEs: --- Trying classification methods --- trying a quadrature trying 1st order linear <- 1st order linear successful -> Calling odsolve with the ODE`, diff(y(x), x) = (-2*x*y(x)^3+c__1-x)/x, y(x)` *** Sublevel 2 *** Methods for first order ODEs: --- Trying classification methods --- trying a quadrature trying 1st order linear trying Bernoulli trying separable trying inverse linear trying homogeneous types: trying Chini differential order: 1; looking for linear symmetries trying exact trying Abel Looking for potential symmetries Looking for potential symmetries Looking for potential symmetries trying inverse_Riccati trying an equivalence to an Abel ODE differential order: 1; trying a linearization to 2nd order --- trying a change of variables {x -> y(x), y(x) -> x} differential order: 1; trying a linearization to 2nd order trying 1st order ODE linearizable_by_differentiation --- Trying Lie symmetry methods, 1st order --- `, `-> Computing symmetries using: way = 3 `, `-> Computing symmetries using: way = 4 `, `-> Computing symmetries using: way = 2 trying symmetry patterns for 1st order ODEs -> trying a symmetry pattern of the form [F(x)*G(y), 0] -> trying a symmetry pattern of the form [0, F(x)*G(y)] -> trying symmetry patterns of the forms [F(x),G(y)] and [G(y),F(x)] -> trying a symmetry pattern of the form [F(x),G(x)] -> trying a symmetry pattern of the form [F(y),G(y)] -> trying a symmetry pattern of the form [F(x)+G(y), 0] -> trying a symmetry pattern of the form [0, F(x)+G(y)] -> trying a symmetry pattern of the form [F(x),G(x)*y+H(x)] -> trying a symmetry pattern of conformal type -> trying 2nd order, No Point Symmetries Class V --- trying a change of variables {x -> y(x), y(x) -> x} and re-entering methods for dynamical symmetries --- -> trying 2nd order, No Point Symmetries Class V -> trying 2nd order, No Point Symmetries Class V --- trying a change of variables {x -> y(x), y(x) -> x} and re-entering methods for dynamical symmetries --- -> trying 2nd order, No Point Symmetries Class V -> trying 2nd order, No Point Symmetries Class V --- trying a change of variables {x -> y(x), y(x) -> x} and re-entering methods for dynamical symmetries --- -> trying 2nd order, No Point Symmetries Class V trying 2nd order, integrating factor of the form mu(x,y)/(y)^n, only the general case -> trying 2nd order, dynamical_symmetries, only a reduction of order through one integrating factor of the form G(x,y)/(1+H(x,y)*y)^ --- trying a change of variables {x -> y(x), y(x) -> x} and re-entering methods for dynamical symmetries --- -> trying 2nd order, dynamical_symmetries, only a reduction of order through one integrating factor of the form G(x,y)/(1+H(x,y)* --- Trying Lie symmetry methods, 2nd order --- `, `-> Computing symmetries using: way = 3 `, `-> Computing symmetries using: way = 5 `, `-> Computing symmetries using: way = formal`
✗ Solution by Maple
dsolve(x*diff(y(x),x$2)+(6*x*y(x)^2+1)*diff(y(x),x)+2*y(x)^3+1=0,y(x), singsol=all)
\[ \text {No solution found} \]
✗ Solution by Mathematica
Time used: 0.0 (sec). Leaf size: 0
DSolve[x*y''[x]+(6*x*y[x]^2+1)*y'[x]+2*y[x]^3+1==0,y[x],x,IncludeSingularSolutions -> True]
Not solved