3.226 problem 1230

Internal problem ID [9559]
Internal file name [OUTPUT/8499_Monday_June_06_2022_03_14_44_AM_17303029/index.tex]

Book: Differential Gleichungen, E. Kamke, 3rd ed. Chelsea Pub. NY, 1948
Section: Chapter 2, linear second order
Problem number: 1230.
ODE order: 2.
ODE degree: 1.

The type(s) of ODE detected by this program : "exact linear second order ode", "second_order_integrable_as_is"

Maple gives the following as the ode type

[[_2nd_order, _exact, _linear, _homogeneous]]

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

3.226.1 Solving as second order integrable as is ode

Integrating both sides of the ODE w.r.t \(x\) gives \begin {align*} \int \left (\left (x^{2}+1\right ) y^{\prime \prime }+a x y^{\prime }+\left (a -2\right ) y\right )d x &= 0 \\ \left (a x -2 x \right ) y+\left (x^{2}+1\right ) y^{\prime } = c_{1} \end {align*}

Which is now solved for \(y\).

Entering Linear first order ODE solver. In canonical form a linear first order is \begin {align*} y^{\prime } + p(x)y &= q(x) \end {align*}

Where here \begin {align*} p(x) &=\frac {a x -2 x}{x^{2}+1}\\ q(x) &=\frac {c_{1}}{x^{2}+1} \end {align*}

Hence the ode is \begin {align*} y^{\prime }+\frac {\left (a x -2 x \right ) y}{x^{2}+1} = \frac {c_{1}}{x^{2}+1} \end {align*}

The integrating factor \(\mu \) is \begin{align*} \mu &= {\mathrm e}^{\int \frac {a x -2 x}{x^{2}+1}d x} \\ &= {\mathrm e}^{\frac {\left (a -2\right ) \ln \left (x^{2}+1\right )}{2}} \\ \end{align*} Which simplifies to \[ \mu = \left (x^{2}+1\right )^{\frac {a}{2}-1} \] The ode becomes \begin {align*} \frac {\mathop {\mathrm {d}}}{ \mathop {\mathrm {d}x}}\left ( \mu y\right ) &= \left (\mu \right ) \left (\frac {c_{1}}{x^{2}+1}\right ) \\ \frac {\mathop {\mathrm {d}}}{ \mathop {\mathrm {d}x}} \left (\left (x^{2}+1\right )^{\frac {a}{2}-1} y\right ) &= \left (\left (x^{2}+1\right )^{\frac {a}{2}-1}\right ) \left (\frac {c_{1}}{x^{2}+1}\right )\\ \mathrm {d} \left (\left (x^{2}+1\right )^{\frac {a}{2}-1} y\right ) &= \left (c_{1} \left (x^{2}+1\right )^{-2+\frac {a}{2}}\right )\, \mathrm {d} x \end {align*}

Integrating gives \begin {align*} \left (x^{2}+1\right )^{\frac {a}{2}-1} y &= \int {c_{1} \left (x^{2}+1\right )^{-2+\frac {a}{2}}\,\mathrm {d} x}\\ \left (x^{2}+1\right )^{\frac {a}{2}-1} y &= c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right ) + c_{2} \end {align*}

Dividing both sides by the integrating factor \(\mu =\left (x^{2}+1\right )^{\frac {a}{2}-1}\) results in \begin {align*} y &= \left (x^{2}+1\right )^{-\frac {a}{2}+1} c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \left (x^{2}+1\right )^{-\frac {a}{2}+1} \end {align*}

which simplifies to \begin {align*} y &= \left (c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \right ) \left (x^{2}+1\right )^{-\frac {a}{2}+1} \end {align*}

Summary

The solution(s) found are the following \begin{align*} \tag{1} y &= \left (c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \right ) \left (x^{2}+1\right )^{-\frac {a}{2}+1} \\ \end{align*}

Verification of solutions

\[ y = \left (c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \right ) \left (x^{2}+1\right )^{-\frac {a}{2}+1} \] Verified OK.

3.226.2 Solving as type second_order_integrable_as_is (not using ABC version)

Writing the ode as \[ \left (x^{2}+1\right ) y^{\prime \prime }+a x y^{\prime }+\left (a -2\right ) y = 0 \] Integrating both sides of the ODE w.r.t \(x\) gives \begin {align*} \int \left (\left (x^{2}+1\right ) y^{\prime \prime }+a x y^{\prime }+\left (a -2\right ) y\right )d x &= 0 \\ y^{\prime }+a x y+y^{\prime } x^{2}-2 y x = c_{1} \end {align*}

Which is now solved for \(y\).

Entering Linear first order ODE solver. In canonical form a linear first order is \begin {align*} y^{\prime } + p(x)y &= q(x) \end {align*}

Where here \begin {align*} p(x) &=\frac {a x -2 x}{x^{2}+1}\\ q(x) &=\frac {c_{1}}{x^{2}+1} \end {align*}

Hence the ode is \begin {align*} y^{\prime }+\frac {\left (a x -2 x \right ) y}{x^{2}+1} = \frac {c_{1}}{x^{2}+1} \end {align*}

The integrating factor \(\mu \) is \begin{align*} \mu &= {\mathrm e}^{\int \frac {a x -2 x}{x^{2}+1}d x} \\ &= {\mathrm e}^{\frac {\left (a -2\right ) \ln \left (x^{2}+1\right )}{2}} \\ \end{align*} Which simplifies to \[ \mu = \left (x^{2}+1\right )^{\frac {a}{2}-1} \] The ode becomes \begin {align*} \frac {\mathop {\mathrm {d}}}{ \mathop {\mathrm {d}x}}\left ( \mu y\right ) &= \left (\mu \right ) \left (\frac {c_{1}}{x^{2}+1}\right ) \\ \frac {\mathop {\mathrm {d}}}{ \mathop {\mathrm {d}x}} \left (\left (x^{2}+1\right )^{\frac {a}{2}-1} y\right ) &= \left (\left (x^{2}+1\right )^{\frac {a}{2}-1}\right ) \left (\frac {c_{1}}{x^{2}+1}\right )\\ \mathrm {d} \left (\left (x^{2}+1\right )^{\frac {a}{2}-1} y\right ) &= \left (c_{1} \left (x^{2}+1\right )^{-2+\frac {a}{2}}\right )\, \mathrm {d} x \end {align*}

Integrating gives \begin {align*} \left (x^{2}+1\right )^{\frac {a}{2}-1} y &= \int {c_{1} \left (x^{2}+1\right )^{-2+\frac {a}{2}}\,\mathrm {d} x}\\ \left (x^{2}+1\right )^{\frac {a}{2}-1} y &= c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right ) + c_{2} \end {align*}

Dividing both sides by the integrating factor \(\mu =\left (x^{2}+1\right )^{\frac {a}{2}-1}\) results in \begin {align*} y &= \left (x^{2}+1\right )^{-\frac {a}{2}+1} c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \left (x^{2}+1\right )^{-\frac {a}{2}+1} \end {align*}

which simplifies to \begin {align*} y &= \left (c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \right ) \left (x^{2}+1\right )^{-\frac {a}{2}+1} \end {align*}

Summary

The solution(s) found are the following \begin{align*} \tag{1} y &= \left (c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \right ) \left (x^{2}+1\right )^{-\frac {a}{2}+1} \\ \end{align*}

Verification of solutions

\[ y = \left (c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \right ) \left (x^{2}+1\right )^{-\frac {a}{2}+1} \] Verified OK.

3.226.3 Solving as exact linear second order ode ode

An ode of the form \begin {align*} p \left (x \right ) y^{\prime \prime }+q \left (x \right ) y^{\prime }+r \left (x \right ) y&=s \left (x \right ) \end {align*}

is exact if \begin {align*} p''(x) - q'(x) + r(x) &= 0 \tag {1} \end {align*}

For the given ode we have \begin {align*} p(x) &= x^{2}+1\\ q(x) &= a x\\ r(x) &= a -2\\ s(x) &= 0 \end {align*}

Hence \begin {align*} p''(x) &= 2\\ q'(x) &= a \end {align*}

Therefore (1) becomes \begin {align*} 2- \left (a\right ) + \left (a -2\right )&=0 \end {align*}

Hence the ode is exact. Since we now know the ode is exact, it can be written as \begin {align*} \left (p \left (x \right ) y^{\prime }+\left (q \left (x \right )-p^{\prime }\left (x \right )\right ) y\right )' &= s(x) \end {align*}

Integrating gives \begin {align*} p \left (x \right ) y^{\prime }+\left (q \left (x \right )-p^{\prime }\left (x \right )\right ) y&=\int {s \left (x \right )\, dx} \end {align*}

Substituting the above values for \(p,q,r,s\) gives \begin {align*} \left (a x -2 x \right ) y+\left (x^{2}+1\right ) y^{\prime }&=c_{1} \end {align*}

We now have a first order ode to solve which is \begin {align*} \left (a x -2 x \right ) y+\left (x^{2}+1\right ) y^{\prime } = c_{1} \end {align*}

Entering Linear first order ODE solver. In canonical form a linear first order is \begin {align*} y^{\prime } + p(x)y &= q(x) \end {align*}

Where here \begin {align*} p(x) &=\frac {a x -2 x}{x^{2}+1}\\ q(x) &=\frac {c_{1}}{x^{2}+1} \end {align*}

Hence the ode is \begin {align*} y^{\prime }+\frac {\left (a x -2 x \right ) y}{x^{2}+1} = \frac {c_{1}}{x^{2}+1} \end {align*}

The integrating factor \(\mu \) is \begin{align*} \mu &= {\mathrm e}^{\int \frac {a x -2 x}{x^{2}+1}d x} \\ &= {\mathrm e}^{\frac {\left (a -2\right ) \ln \left (x^{2}+1\right )}{2}} \\ \end{align*} Which simplifies to \[ \mu = \left (x^{2}+1\right )^{\frac {a}{2}-1} \] The ode becomes \begin {align*} \frac {\mathop {\mathrm {d}}}{ \mathop {\mathrm {d}x}}\left ( \mu y\right ) &= \left (\mu \right ) \left (\frac {c_{1}}{x^{2}+1}\right ) \\ \frac {\mathop {\mathrm {d}}}{ \mathop {\mathrm {d}x}} \left (\left (x^{2}+1\right )^{\frac {a}{2}-1} y\right ) &= \left (\left (x^{2}+1\right )^{\frac {a}{2}-1}\right ) \left (\frac {c_{1}}{x^{2}+1}\right )\\ \mathrm {d} \left (\left (x^{2}+1\right )^{\frac {a}{2}-1} y\right ) &= \left (c_{1} \left (x^{2}+1\right )^{-2+\frac {a}{2}}\right )\, \mathrm {d} x \end {align*}

Integrating gives \begin {align*} \left (x^{2}+1\right )^{\frac {a}{2}-1} y &= \int {c_{1} \left (x^{2}+1\right )^{-2+\frac {a}{2}}\,\mathrm {d} x}\\ \left (x^{2}+1\right )^{\frac {a}{2}-1} y &= c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right ) + c_{2} \end {align*}

Dividing both sides by the integrating factor \(\mu =\left (x^{2}+1\right )^{\frac {a}{2}-1}\) results in \begin {align*} y &= \left (x^{2}+1\right )^{-\frac {a}{2}+1} c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \left (x^{2}+1\right )^{-\frac {a}{2}+1} \end {align*}

which simplifies to \begin {align*} y &= \left (c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \right ) \left (x^{2}+1\right )^{-\frac {a}{2}+1} \end {align*}

Summary

The solution(s) found are the following \begin{align*} \tag{1} y &= \left (c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \right ) \left (x^{2}+1\right )^{-\frac {a}{2}+1} \\ \end{align*}

Verification of solutions

\[ y = \left (c_{1} x \operatorname {hypergeom}\left (\left [\frac {1}{2}, 2-\frac {a}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right )+c_{2} \right ) \left (x^{2}+1\right )^{-\frac {a}{2}+1} \] Verified OK.

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)] 
   One independent solution has integrals. Trying a hypergeometric solution free of integrals... 
   -> hyper3: Equivalence to 2F1, 1F1 or 0F1 under a power @ Moebius 
   <- hyper3 successful: received ODE is equivalent to the 2F1 ODE 
   -> Trying to convert hypergeometric functions to elementary form... 
   <- elementary form for at least one hypergeometric solution is achieved - returning with no uncomputed integrals 
<- linear_1 successful`
 

✓ Solution by Maple

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

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

\[ y \left (x \right ) = c_{1} \left (x^{2}+1\right )^{1-\frac {a}{2}}+c_{2} x \operatorname {hypergeom}\left (\left [1, \frac {a}{2}-\frac {1}{2}\right ], \left [\frac {3}{2}\right ], -x^{2}\right ) \]

✓ Solution by Mathematica

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

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

\[ y(x)\to \left (x^2+1\right )^{\frac {1}{2}-\frac {a}{4}} \left (c_1 P_{\frac {a-4}{2}}^{\frac {a-2}{2}}(i x)+c_2 Q_{\frac {a-4}{2}}^{\frac {a-2}{2}}(i x)\right ) \]