Internal problem ID [9411]
Internal file name [OUTPUT/8349_Monday_June_06_2022_02_48_43_AM_22060907/index.tex
]
Book: Differential Gleichungen, E. Kamke, 3rd ed. Chelsea Pub. NY, 1948
Section: Chapter 2, linear second order
Problem number: 1080.
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 {y^{\prime \prime }-\left (\frac {f^{\prime }\left (x \right )}{f \left (x \right )}+2 a \right ) y^{\prime }+\left (\frac {f^{\prime }\left (x \right ) a}{f \left (x \right )}+a^{2}-b^{2} f \left (x \right )^{2}\right ) y=0} \]
Maple trace
`Methods for second order ODEs: --- Trying classification methods --- trying a symmetry of the form [xi=0, eta=F(x)] checking if the LODE is missing y -> Heun: Equivalence to the GHE or one of its 4 confluent cases under a power @ Moebius -> trying a solution of the form r0(x) * Y + r1(x) * Y where Y = exp(int(r(x), dx)) * 2F1([a1, a2], [b1], f) -> Trying changes of variables to rationalize or make the ODE simpler trying a symmetry of the form [xi=0, eta=F(x)] checking if the LODE is missing y -> Heun: Equivalence to the GHE or one of its 4 confluent cases under a power @ Moebius -> trying a solution of the form r0(x) * Y + r1(x) * Y where Y = exp(int(r(x), dx)) * 2F1([a1, a2], [b1], f) trying a symmetry of the form [xi=0, eta=F(x)] trying 2nd order exact linear trying symmetries linear in x and y(x) trying to convert to a linear ODE with constant coefficients <- unable to find a useful change of variables trying a symmetry of the form [xi=0, eta=F(x)] trying 2nd order exact linear trying symmetries linear in x and y(x) trying to convert to a linear ODE with constant coefficients trying 2nd order, integrating factor of the form mu(x,y) trying a symmetry of the form [xi=0, eta=F(x)] checking if the LODE is missing y -> Heun: Equivalence to the GHE or one of its 4 confluent cases under a power @ Moebius -> trying a solution of the form r0(x) * Y + r1(x) * Y where Y = exp(int(r(x), dx)) * 2F1([a1, a2], [b1], f) -> Trying changes of variables to rationalize or make the ODE simpler trying a symmetry of the form [xi=0, eta=F(x)] checking if the LODE is missing y -> Heun: Equivalence to the GHE or one of its 4 confluent cases under a power @ Moebius -> trying a solution of the form r0(x) * Y + r1(x) * Y where Y = exp(int(r(x), dx)) * 2F1([a1, a2], [b1], f) trying a symmetry of the form [xi=0, eta=F(x)] trying 2nd order exact linear trying symmetries linear in x and y(x) trying to convert to a linear ODE with constant coefficients <- unable to find a useful change of variables trying a symmetry of the form [xi=0, eta=F(x)] trying to convert to an ODE of Bessel type -> trying reduction of order to Riccati trying Riccati sub-methods: trying Riccati_symmetries -> 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)] <- symmetry pattern of the form [F(x),G(x)*y+H(x)] successful <- Riccati with symmetry pattern of the form [F(x),G(x)*y+H(x)] successful <- reduction of order to Riccati successful <- reduction of order to Riccati successful`
✓ Solution by Maple
Time used: 0.093 (sec). Leaf size: 65
dsolve(diff(diff(y(x),x),x)-(diff(f(x),x)/f(x)+2*a)*diff(y(x),x)+(a*diff(f(x),x)/f(x)+a^2-b^2*f(x)^2)*y(x)=0,y(x), singsol=all)
\[ y \left (x \right ) = {\mathrm e}^{\int \frac {{\mathrm e}^{2 b \left (\int f \left (x \right )d x \right )} \left (\left (-f \left (x \right ) b +a \right ) {\mathrm e}^{-2 b \left (\int f \left (x \right )d x -c_{1} \right )}-f \left (x \right ) b -a \right )}{-{\mathrm e}^{2 b \left (\int f \left (x \right )d x \right )}+{\mathrm e}^{2 c_{1} b}}d x} c_{2} \]
✓ Solution by Mathematica
Time used: 0.112 (sec). Leaf size: 47
DSolve[y[x]*(a^2 - b^2*f[x]^2 + (a*Derivative[1][f][x])/f[x]) - (2*a + Derivative[1][f][x]/f[x])*y'[x] + y''[x] == 0,y[x],x,IncludeSingularSolutions -> True]
\[ y(x)\to e^{a x} \left (c_1 \exp \left (b \int _1^xf(K[1])dK[1]\right )+c_2 \exp \left (-b \int _1^xf(K[2])dK[2]\right )\right ) \]