27.28 problem 38

Internal problem ID [10872]

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-2 Equation of form \(y''+f(x)y'+g(x)y=0\)
Problem number: 38.
ODE order: 2.
ODE degree: 1.

CAS Maple gives this as type [[_2nd_order, _with_linear_symmetries]]

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

Solution by Maple

Time used: 0.078 (sec). Leaf size: 254

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

\[ y = c_{1} {\mathrm e}^{-\frac {\left (\frac {\left (a^{5} b^{5}+\left (b^{2} a^{2}\right )^{\frac {5}{2}}\right ) \left (a x +\frac {3}{2}\right ) x b}{\left (b^{2} a^{2}\right )^{\frac {5}{2}}}+6 a \right ) x}{6}} \operatorname {HeunT}\left (\frac {b^{3} 3^{\frac {2}{3}} a^{2}}{2 \left (b^{2} a^{2}\right )^{\frac {4}{3}}}, -\frac {6 \sqrt {b^{2}}\, \sqrt {a^{2}}}{b a}, -\frac {b^{2} 3^{\frac {1}{3}}}{4 \left (b^{2} a^{2}\right )^{\frac {2}{3}}}, \frac {3^{\frac {2}{3}} a \,b^{2} \left (2 a x +1\right )}{6 \left (b^{2} a^{2}\right )^{\frac {5}{6}}}\right )+c_{2} {\mathrm e}^{\frac {\left (\frac {\left (a x +\frac {3}{2}\right ) \left (a^{5} b^{5}-\left (b^{2} a^{2}\right )^{\frac {5}{2}}\right ) x b}{\left (b^{2} a^{2}\right )^{\frac {5}{2}}}-6 a \right ) x}{6}} \operatorname {HeunT}\left (\frac {b^{3} 3^{\frac {2}{3}} a^{2}}{2 \left (b^{2} a^{2}\right )^{\frac {4}{3}}}, \frac {6 \sqrt {b^{2}}\, \sqrt {a^{2}}}{b a}, -\frac {b^{2} 3^{\frac {1}{3}}}{4 \left (b^{2} a^{2}\right )^{\frac {2}{3}}}, -\frac {\left (a x +\frac {1}{2}\right ) a 3^{\frac {2}{3}} b^{2}}{3 \left (b^{2} a^{2}\right )^{\frac {5}{6}}}\right ) \]

Solution by Mathematica

Time used: 2.136 (sec). Leaf size: 57

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

\[ y(x)\to e^{-a x} (a x+1) \left (c_2 \int _1^x\frac {e^{-\frac {1}{6} b K[1]^2 (2 a K[1]+3)}}{(a K[1]+1)^2}dK[1]+c_1\right ) \]