32.27 problem 236

Internal problem ID [11060]
Internal file name [OUTPUT/10317_Wednesday_January_24_2024_10_07_09_PM_18207104/index.tex]

Book: Handbook of exact solutions for ordinary differential equations. By Polyanin and Zaitsev. Second edition
Section: Chapter 2, Second-Order Differential Equations. section 2.1.2-7 Equation of form \((a_4 x^4+a_3 x^3+a_2 x^2 x+a_1 x+a_0) y''+f(x)y'+g(x)y=0\)
Problem number: 236.
ODE order: 2.
ODE degree: 1.

The type(s) of ODE detected by this program : "unknown"

Maple gives the following as the ode type

[[_2nd_order, _with_linear_symmetries]]

Unable to solve or complete the solution.

\[ \boxed {\left (x^{2}+1\right )^{2} y^{\prime \prime }+2 x \left (x^{2}+1\right ) y^{\prime }+\left (\left (x^{2}+1\right ) \left (a^{2} x^{2}-\lambda \right )+m^{2}\right ) y=0} \]

Maple trace

`Methods for second order ODEs: 
--- Trying classification methods --- 
trying a quadrature 
checking if the LODE has constant coefficients 
checking if the LODE is of Euler type 
trying a symmetry of the form [xi=0, eta=F(x)] 
checking if the LODE is missing y 
-> Trying a Liouvillian solution using Kovacics algorithm 
<- No Liouvillian solutions exists 
-> Trying a solution in terms of special functions: 
   -> Bessel 
   -> elliptic 
   -> Legendre 
   -> Whittaker 
      -> hyper3: Equivalence to 1F1 under a power @ Moebius 
   -> hypergeometric 
      -> heuristic approach 
      -> hyper3: Equivalence to 2F1, 1F1 or 0F1 under a power @ Moebius 
   -> Mathieu 
      -> Equivalence to the rational form of Mathieu ODE under a power @ Moebius 
trying a solution in terms of MeijerG functions 
-> Heun: Equivalence to the GHE or one of its 4 confluent cases under a power @ Moebius 
<- Heun successful: received ODE is equivalent to the  HeunC  ODE, case  a <> 0, e <> 0, c = 0 `
 

Solution by Maple

Time used: 0.5 (sec). Leaf size: 68

dsolve((x^2+1)^2*diff(y(x),x$2)+2*x*(x^2+1)*diff(y(x),x)+( (x^2+1)*(a^2*x^2-lambda)+m^2)*y(x)=0,y(x), singsol=all)
 

\[ y \left (x \right ) = \left (\operatorname {HeunC}\left (0, \frac {1}{2}, m , -\frac {a^{2}}{4}, \frac {1}{4}+\frac {m^{2}}{4}-\frac {\lambda }{4}, -x^{2}\right ) c_{2} x +\operatorname {HeunC}\left (0, -\frac {1}{2}, m , -\frac {a^{2}}{4}, \frac {1}{4}+\frac {m^{2}}{4}-\frac {\lambda }{4}, -x^{2}\right ) c_{1} \right ) \left (x^{2}+1\right )^{\frac {m}{2}} \]

Solution by Mathematica

Time used: 0.605 (sec). Leaf size: 124

DSolve[(x^2+1)^2*y''[x]+2*x*(x^2+1)*y'[x]+( (x^2+1)*(a^2*x^2-\[Lambda])+m^2)*y[x]==0,y[x],x,IncludeSingularSolutions -> True]
 

\[ y(x)\to \left (x^2+1\right )^{\frac {\sqrt {m^2}}{2}} \left (c_2 x \text {HeunC}\left [\frac {1}{4} \left (\lambda -m^2-3 \sqrt {m^2}-2\right ),-\frac {a^2}{4},\frac {3}{2},\sqrt {m^2}+1,0,-x^2\right ]+c_1 \text {HeunC}\left [\frac {1}{4} \left (\lambda -m^2-\sqrt {m^2}\right ),-\frac {a^2}{4},\frac {1}{2},\sqrt {m^2}+1,0,-x^2\right ]\right ) \]