2.1.64 Problem 64

Internal problem ID [8776]
Book : Own collection of miscellaneous problems
Section : section 1.0
Problem number : 64
Date solved : Sunday, March 30, 2025 at 01:32:58 PM
CAS classification : [[_Emden, _Fowler], [_2nd_order, _with_linear_symmetries]]

✓ Maple. Time used: 0.052 (sec). Leaf size: 106

ode:=y(x)^2*diff(diff(y(x),x),x) = x; 
dsolve(ode,y(x), singsol=all);
 

$$y=\mathrm{RootOf}(\ln \left(x\right)+2^{1/3}{\int }_{}^{\mathit{\text{\_Z}}}\frac{1}{2^{1/3}\mathit{\text{\_f}}+2\mathrm{RootOf}(\mathrm{AiryBi}\left(\frac{2{\mathit{\text{\_Z}}}^2\mathit{\text{\_f}}+2^{2/3}}{2\mathit{\text{\_f}}}\right)c_1\mathit{\text{\_Z}}+\mathit{\text{\_Z}}\mathrm{AiryAi}\left(\frac{2{\mathit{\text{\_Z}}}^2\mathit{\text{\_f}}+2^{2/3}}{2\mathit{\text{\_f}}}\right)+\mathrm{AiryBi}(1,\frac{2{\mathit{\text{\_Z}}}^2\mathit{\text{\_f}}+2^{2/3}}{2\mathit{\text{\_f}}})c_1+\mathrm{AiryAi}(1,\frac{2{\mathit{\text{\_Z}}}^2\mathit{\text{\_f}}+2^{2/3}}{2\mathit{\text{\_f}}}))}d\mathit{\text{\_f}}-c_2)x$$

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 2nd order, integrating factors of the form mu(x,y)/(y)^n, only the sing\ 
ular cases 
trying symmetries linear in x and y(x) 
trying differential order: 2; exact nonlinear 
trying 2nd order, integrating factor of the form mu(y) 
trying 2nd order, integrating factor of the form mu(x,y) 
trying 2nd order, integrating factor of the form mu(x,y)/(y)^n, only the gener\ 
al case 
trying 2nd order, integrating factor of the form mu(y,y) 
-> Calling odsolve with the ODE, -(_y1^3-4)/_y1^3*y(x)+2/_y1^3*(x*diff(y(x),x)+ 
_y1), y(x) 
   *** Sublevel 2 *** 
   Methods for first order ODEs: 
   --- Trying classification methods --- 
   trying a quadrature 
   trying 1st order linear 
   <- 1st order linear successful 
trying differential order: 2; mu polynomial in y 
trying 2nd order, integrating factor of the form mu(x,y) 
differential order: 2; looking for linear symmetries 
differential order: 2; found: 1 linear symmetries. Trying reduction of order 
   2nd order, trying reduction of order with given symmetries: 
[x, y] 
   -> Calling odsolve with the ODE, diff(y(x),x) = y(x)/x, y(x) 
      *** Sublevel 2 *** 
      Methods for first order ODEs: 
      --- Trying classification methods --- 
      trying a quadrature 
      trying 1st order linear 
      <- 1st order linear successful 
   -> Computing symmetries using: way = 3 
   -> Computing symmetries using: way = exp_sym 
Try integration with the canonical coordinates of the symmetry [x, y] 
-> Calling odsolve with the ODE, diff(_b(_a),_a) = -_b(_a)^2*(-_a^2+_b(_a))/_a^ 
2, _b(_a), explicit 
   *** 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 
   <- Abel successful 
<- differential order: 2; canonical coordinates successful
 

✗ Mathematica

ode=y[x]^2*D[y[x],{x,2}]==x; 
ic={}; 
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
 

Not solved

✗ Sympy

from sympy import * 
x = symbols("x") 
y = Function("y") 
ode = Eq(-x + y(x)**2*Derivative(y(x), (x, 2)),0) 
ics = {} 
dsolve(ode,func=y(x),ics=ics)
 

NotImplementedError : solve: Cannot solve -x + y(x)**2*Derivative(y(x), (x, 2))