60.5.25 problem 1562

Internal problem ID [11522]
Book : Differential Gleichungen, E. Kamke, 3rd ed. Chelsea Pub. NY, 1948
Section : Chapter 4, linear fourth order
Problem number : 1562
Date solved : Sunday, March 30, 2025 at 08:24:08 PM
CAS classification : [[_high_order, _with_linear_symmetries]]

\begin{align*} x^{4} y^{\prime \prime \prime \prime }+4 x^{3} y^{\prime \prime \prime }-\left (4 n^{2}-1\right ) x^{2} y^{\prime \prime }+\left (4 n^{2}-1\right ) x y^{\prime }-4 y x^{4}&=0 \end{align*}

✓ Maple. Time used: 0.022 (sec). Leaf size: 65

ode:=x^4*diff(diff(diff(diff(y(x),x),x),x),x)+4*x^3*diff(diff(diff(y(x),x),x),x)-(4*n^2-1)*x^2*diff(diff(y(x),x),x)+(4*n^2-1)*x*diff(y(x),x)-4*y(x)*x^4 = 0; 
dsolve(ode,y(x), singsol=all);
 

\[ y = \left (\operatorname {BesselY}\left (n , \left (\frac {1}{2}-\frac {i}{2}\right ) \sqrt {2}\, x \right ) c_3 +c_1 \operatorname {BesselJ}\left (n , \left (\frac {1}{2}-\frac {i}{2}\right ) \sqrt {2}\, x \right )\right ) \operatorname {BesselJ}\left (n , \left (\frac {1}{2}+\frac {i}{2}\right ) \sqrt {2}\, x \right )+\operatorname {BesselY}\left (n , \left (\frac {1}{2}+\frac {i}{2}\right ) \sqrt {2}\, x \right ) \left (\operatorname {BesselY}\left (n , \left (\frac {1}{2}-\frac {i}{2}\right ) \sqrt {2}\, x \right ) c_4 +c_2 \operatorname {BesselJ}\left (n , \left (\frac {1}{2}-\frac {i}{2}\right ) \sqrt {2}\, x \right )\right ) \]

✓ Mathematica. Time used: 0.509 (sec). Leaf size: 140

ode=-4*x^4*y[x] + (-1 + 4*n^2)*x*D[y[x],x] - (-1 + 4*n^2)*x^2*D[y[x],{x,2}] + 4*x^3*Derivative[3][y][x] + x^4*Derivative[4][y][x] == 0; 
ic={}; 
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
 

\[ y(x)\to c_1 \, _0F_3\left (;\frac {1}{2},1-\frac {n}{2},\frac {n}{2}+1;\frac {x^4}{64}\right )+\frac {1}{8} i c_2 x^2 \, _0F_3\left (;\frac {3}{2},\frac {3}{2}-\frac {n}{2},\frac {n}{2}+\frac {3}{2};\frac {x^4}{64}\right )+c_3 \left (\frac {i}{2}\right )^{-n} \operatorname {Gamma}(1-n)^2 \left (\text {ber}_{-n}(x){}^2+\text {bei}_{-n}(x){}^2\right )+c_4 \left (\frac {i}{2}\right )^n \operatorname {Gamma}(n+1)^2 \left (\text {ber}_n(x){}^2+\text {bei}_n(x){}^2\right ) \]

✗ Sympy

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

NotImplementedError : The given ODE -x*(4*n**2*Derivative(y(x), (x, 2)) + 4*x**2*y(x) - x**2*Derivative(y(x), (x, 4)) - 4*x*Derivative(y(x), (x, 3)) - Derivative(y(x), (x, 2)))/(4*n**2 - 1) + Derivative(y(x), x) cannot be solved by the factorable group method