Problem Solve x^{2}y^{\prime }+y^{2}=xyy^{\prime }
Solution
Rewriting the ODE as \begin{equation} y^{\prime }\left ( x^{2}-xy\right ) +y^{2}=0\tag{1} \end{equation} Dividing by x^{2}\neq 0 gives \frac{dy}{dx}\left ( 1-\frac{y}{x}\right ) +\frac{y^{2}}{x^{2}}=0 We see this is homogeneous of order 1. This can be confirmed by writing the above as\begin{align*} dy\left ( x^{2}-xy\right ) +y^{2}dx & =0\\ dyx^{2}-xydy+y^{2}dx & =0 \end{align*}
We want to find if a weight m can be found, so that the substitution y=vx^{m} makes the above ODE separable. To find m, we assign weight m to both y and dy, and a weight of 1 to both x and dx, and then try to find if there is an m which makes each term sums to the same total weight. (in other words, we want each term units to be the same).
The term \left ( dy\right ) \left ( x^{2}\right ) has total weight of m+2 (it is the exponents that we add). And the term \left ( x\right ) \left ( y\right ) \left ( dy\right ) has total weight 1+2m and the last term \left ( y^{2}\right ) \left ( dx\right ) has weight 2m+1. Therefore we have this result for the weight of each term (there are 3 terms above). \left \{ m+2,1+2m,1+2m\right \} We see that if m=1 then each term will have the same total weight of 3 giving \left \{ 3,3,3\right \} . So this is homogenous ODE of order m=1. Now that we know the weight, we use the substitution \fbox{$y=vx$} Hence y^{\prime }=v^{\prime }x+v. Substituting these back into (1) gives a new ODE in v which is separable. If it is not separable, it means we made a mistake somewhere.\begin{align*} \left ( v^{\prime }x+v\right ) \left ( x^{2}-x^{2}v\right ) +v^{2}x^{2} & =0\\ v^{\prime }x^{3}-v^{\prime }vx^{3}+vx^{2}-x^{2}v^{2}+v^{2}x^{2} & =0\\ v^{\prime }x^{3}-v^{\prime }vx^{3}+vx^{2} & =0 \end{align*}
Dividing by x^{3} for x\neq 0 gives\begin{align*} v^{\prime }-v^{\prime }v+\frac{v}{x} & =0\\ v^{\prime }\left ( 1-v\right ) & =-\frac{v}{x}\\ \frac{dv}{dx}\frac{\left ( 1-v\right ) }{v} & =-\frac{1}{x}\\ dv\frac{\left ( 1-v\right ) }{v} & =-\frac{1}{x}dx\\ dv\frac{\left ( v-1\right ) }{v} & =\frac{1}{x}dx \end{align*}
Integrating both sides gives\begin{align*} \int v-\frac{1}{v}dv & =\int \frac{1}{x}dx\\ v-\ln v & =\ln x+C \end{align*}
Taking exponential of both sides gives\begin{align*} e^{v-\ln v} & =Cx\\ \frac{e^{v}}{v} & =Cx \end{align*}
But v=\frac{y}{x}. Therefore the above becomes\begin{align*} \frac{e^{\frac{y}{x}}}{\frac{y}{x}} & =Cx\\ \frac{e^{\frac{y}{x}}}{y} & =C \end{align*}
Hence the solution is\begin{equation} y=C_{1}e^{\frac{y}{x}}\qquad x\neq 0 \tag{2} \end{equation} Where C_{1} is the constant of integration. y can not be solved for directly in the above. But we can solve for x in terms of y if needed as follows\begin{align} \ln y & =\ln C_{1}+\frac{y}{x}\nonumber \\ \ln y-C_{2} & =\frac{y}{x}\nonumber \\ x & =\frac{y}{\ln y-C_{2}} \tag{3} \end{align}
Problem Solve y^{\prime }=\frac{a^{2}}{\left ( x+y\right ) ^{2}}
Solution
Let v\left ( x\right ) =x+y\left ( x\right ) . Hence v^{\prime }=1+y^{\prime } or y^{\prime }=v^{\prime }-1. Substituting this back into the ODE gives\begin{align*} v^{\prime }-1 & =\frac{a^{2}}{v^{2}}\\ \frac{dv}{dx} & =\frac{a^{2}}{v^{2}}+1 \end{align*}
This is separable.\begin{align*} \frac{dv}{\frac{a^{2}}{v^{2}}+1} & =dx\\ \frac{v^{2}}{a^{2}+v^{2}}dv & =dx \end{align*}
By long division \frac{v^{2}}{a^{2}+v^{2}}=1-\frac{a^{2}}{a^{2}+v^{2}}. The above becomes \left ( 1-\frac{a^{2}}{a^{2}+v^{2}}\right ) dv=dx Integrating both sides gives\begin{align*} \int \left ( 1-\frac{a^{2}}{a^{2}+v^{2}}\right ) dv & =\int dx\\ \int dv-a^{2}\int \frac{1}{a^{2}+v^{2}}dv & =\int dx \end{align*}
But \int \frac{1}{a^{2}+v^{2}}dv=\frac{1}{a^{2}}\int \frac{1}{1+\left ( \frac{v}{a}\right ) ^{2}}dv=\frac{1}{a^{2}}\left ( a\arctan \left ( \frac{v}{a}\right ) \right ) =\frac{1}{a}\arctan \left ( \frac{v}{a}\right ) , hence the above becomes\begin{align*} v-a^{2}\left ( \frac{1}{a}\arctan \left ( \frac{v}{a}\right ) \right ) & =x+C\\ v-a\arctan \left ( \frac{v}{a}\right ) & =x+C\\ a\arctan \left ( \frac{v}{a}\right ) & =v-x-C\\ \arctan \left ( \frac{v}{a}\right ) & =\frac{v-x}{a}+C_{1} \end{align*}
Where C_{1}=\frac{-C}{a}, a new constant. Taking the \tan of both sides gives \frac{v}{a}=\tan \left ( \frac{v-x}{a}+C_{1}\right ) But v=x+y, and the above becomes\begin{align*} \frac{x+y}{a} & =\tan \left ( \frac{\left ( x+y\right ) -x}{a}+C_{1}\right ) \\ \frac{x+y}{a} & =\tan \left ( \frac{y}{a}+C_{1}\right ) \end{align*}
Therefore the final solution is y=a\tan \left ( \frac{y}{a}+C_{1}\right ) -x\qquad a\neq 0 Where C_{1} is arbitrary constant.
Problem Solve y^{\prime \prime }+\left ( y^{\prime }\right ) ^{2}+1=0
Solution
Since y is missing from the ODE, we can convert this to a first order using y^{\prime }=p\left ( x\right ) . Therefore y^{\prime \prime }=\frac{dp}{dx} and the ODE becomes\begin{align*} \frac{dp}{dx}+p^{2}+1 & =0\\ \frac{dp}{dx} & =-\left ( 1+p^{2}\right ) \\ \frac{dp}{1+p^{2}} & =-dx \end{align*}
Integrating both sides gives\begin{align*} \int \frac{dp}{1+p^{2}} & =-\int dx\\ \arctan \left ( p\right ) & =-x+C_{1}\\ p & =\tan \left ( -x+C_{1}\right ) \end{align*}
But p=y^{\prime }. Hence we need now to solve \frac{dy}{dx}=\tan \left ( -x+C_{1}\right ) . Integrating both sides gives\begin{align*} y & =\int \tan \left ( -x+C_{1}\right ) dx\\ & =\int \frac{\sin \left ( -x+C_{1}\right ) }{\cos \left ( -x+C_{1}\right ) }dx\\ & =\int \frac{-\sin \left ( x-C_{1}\right ) }{\cos \left ( x-C_{1}\right ) }dx\\ & =\int \frac{\frac{d}{dx}\left ( \cos \left ( x-C_{1}\right ) \right ) }{\cos \left ( x-C_{1}\right ) }dx \end{align*}
But \int \frac{V^{\prime }}{V}dx=\ln \left ( V\right ) , hence the above becomes y=\ln \left ( \cos \left ( x-C_{1}\right ) \right ) +C_{2} Replacing -C_{1} by new constant C_{3}, the final solution becomes y=\ln \left ( \cos \left ( x+C_{3}\right ) \right ) +C_{2} Where C_{2},C_{3} are constants of integration.
Problem Solve xy^{\prime }+y+x^{4}y^{4}e^{x}=0
Solution
Dividing by x\neq 0 and rewriting gives\begin{equation} y^{\prime }+\frac{1}{x}y=\left ( -x^{3}e^{x}\right ) y^{4}\tag{1} \end{equation} A Bernoulli ODE has the form y^{\prime }+a\left ( x\right ) y=b\left ( x\right ) y^{n} where n\neq 1. Comparing the above to Bernoulli ODE form, show it is Bernoulli ODE where a\left ( x\right ) =\frac{1}{x},b\left ( x\right ) =-x^{3}e^{x}. Dividing (1) by y^{4} gives \frac{1}{y^{4}}y^{\prime }+\frac{1}{x}y^{-3}=-x^{3}e^{x} Letting v=y^{-3} or \frac{dv}{dx}=-3y^{-4}\frac{dy}{dx}. Hence \frac{dy}{dx}=-\frac{dv}{dx}\frac{y^{4}}{3}. Substituting this in the above gives\begin{align*} \frac{1}{y^{4}}\left ( -\frac{dv}{dx}\frac{y^{4}}{3}\right ) +\frac{1}{x}v & =-x^{3}e^{x}\\ -\frac{1}{3}\frac{dv}{dx}+\frac{1}{x}v & =-x^{3}e^{x}\\ \frac{dv}{dx}-\frac{3}{x}v & =3x^{3}e^{x} \end{align*}
This is now linear in v. The integrating factor \mu =e^{\int -\frac{3}{x}dx}=e^{-3\ln x}=\frac{1}{x^{3}}. Multiplying both sides of the above by this integrating factor making the left side complete differential\begin{align*} \frac{d}{dx}\left ( \frac{1}{x^{3}}v\right ) & =\frac{1}{x^{3}}3x^{3}e^{x}\\ \frac{d}{dx}\left ( \frac{1}{x^{3}}v\right ) & =3e^{x} \end{align*}
Integrating gives\begin{align*} \frac{1}{x^{3}}v & =3e^{x}+C\\ v & =3x^{3}e^{x}+Cx^{3}\\ & =x^{3}\left ( 3e^{x}+C\right ) \end{align*}
But v=y^{-3}, hence the above becomes\begin{align*} \frac{1}{y^{3}} & =x^{3}\left ( 3e^{x}+C\right ) \\ y^{3} & =\frac{1}{x^{3}\left ( 3e^{x}+C\right ) } \end{align*}
This shows that there are 3 solutions since the above is a cubic equation. But we can leave the solution in implicit form\begin{align*} y & =\sqrt [\frac{1}{3}]{\frac{1}{x^{3}\left ( 3e^{x}+C\right ) }}\\ & =\frac{1}{x}\sqrt [\frac{1}{3}]{\frac{1}{3e^{x}+C}} \end{align*}
Problem Find both a general solution and a singular solution of
x^{2}\left ( y^{\prime }\right ) ^{2}-2\left ( xy-4\right ) y^{\prime }+y^{2}=0 Solution
Rewriting the ODE as y^{2}-2xyy^{\prime }+8y^{\prime }+x^{2}\left ( y^{\prime }\right ) ^{2}=0 Let y^{\prime }=p and the above becomes y^{2}+y\left ( -2xp\right ) +\left ( 8p+x^{2}p^{2}\right ) =0 This is quadratic in y. Solving for y=\frac{-b}{2a}\pm \frac{1}{2a}\sqrt{b^{2}-4ac}\begin{align*} y & =xp\pm \frac{1}{2}\sqrt{4x^{2}p^{2}-4\left ( 8p+x^{2}p^{2}\right ) }\\ & =xp\pm \sqrt{x^{2}p^{2}-8p-x^{2}p^{2}}\\ & =xp\pm 2\sqrt{-2p} \end{align*}
case one \begin{align} y & =xp+2\sqrt{-2p}\nonumber \\ & =xp+f\left ( p\right ) \tag{1} \end{align}
This can be written as y=G\left ( x,p\right ) Where G\left ( x,p\right ) =xp+f\left ( p\right ) . This form of ODE is called the Clairaut ODE. Taking derivative w.r.t. x gives y^{\prime }=\frac{\partial G}{\partial x}+\frac{\partial G}{\partial p}\frac{dp}{dx} But y^{\prime }=p and the above becomes p=\frac{\partial G}{\partial x}+\frac{\partial G}{\partial p}\frac{dp}{dx} But \frac{\partial G}{\partial x}=p, hence the above reduces to\begin{equation} 0=\frac{\partial G}{\partial p}\frac{dp}{dx}\tag{2} \end{equation} Then either \frac{\partial G}{\partial p}=0 or \frac{dp}{dx}=0.
When \frac{dp}{dx}=0 or y^{\prime \prime }=0 therefore the solution is\begin{equation} y=C_{1}x+C_{2} \tag{3} \end{equation} But we are solving a first order ODE. So we expect it to have one constant of integration only. By comparing (3) with equation (1) which is y=xp+f\left ( p\right ) shows that C_{2}=f\left ( C_{1}\right ) =2\sqrt{-2C_{1}} Then the solution will now contain one constant of integration C_{1}. Hence the first solution is y=C_{1}x+2\sqrt{-2C_{1}} The second possibility comes from \frac{\partial G}{\partial p}=0. This gives \begin{align*} x+f^{\prime }\left ( p\right ) & =0\\ x+2\frac{d}{dp}\left ( -2p\right ) ^{\frac{1}{2}} & =0\\ x+2\frac{1}{2}\left ( -2p\right ) ^{\frac{-1}{2}}\left ( -2\right ) & =0\\ x-\frac{2}{\sqrt{-2p}} & =0\\ x\sqrt{-2p} & =2\\ -2px^{2} & =4\\ p & =-\frac{2}{x^{2}} \end{align*}
Now that we found p, we substitute it back into (1) given by y=xp+2\sqrt{-2p}. Hence the second solution is found directly as follows\begin{align*} y & =xp+2\sqrt{-2p}\\ & =-\frac{2}{x}+2\sqrt{-2\left ( -\frac{2}{x^{2}}\right ) }\\ & =-\frac{2}{x}+2\sqrt{\frac{4}{x^{2}}}\\ & =-\frac{2}{x}+\frac{4}{x}\\ & =\frac{2}{x} \end{align*}
Summary of case one From above we obtained the following two solutions\begin{align*} y_{1} & =C_{1}x+2\sqrt{-2C_{1}}\\ y_{2} & =\frac{2}{x} \end{align*}
Where y_{2}\left ( x\right ) is the singular solution since it can’t be obtained from the first solution with the constants of integrations by changing them to any value.
We now do the same steps for the case of y=xp-2\sqrt{-2p}. This follows the same steps as above as the only difference is the sign and hence the steps will not be repeated. It gives the solution y_{3}=C_{1}x-2\sqrt{-2C_{1}} With the same singular solution. Therefore there are three solutions to this ODE and these are summarized below\begin{align*} y_{1} & =C_{1}x+2\sqrt{-2C_{1}}\\ y_{2} & =\frac{2}{x}\\ y_{3} & =C_{1}x-2\sqrt{-2C_{1}} \end{align*}
With y_{2}\left ( x\right ) being the singular solution. Singular solutions do not have constant of integration in them and can not be obtained from the general solution by any substitution for constants of integration. The general solution contain constant of integrations in them.
Problem
Find the general real solution to the following equation where A\left ( x\right ) is a known function A\left ( x\right ) y^{\prime \prime }+A^{\prime }\left ( x\right ) y^{\prime }+\frac{y}{A\left ( x\right ) }=0 Solution
Let us first assume A\left ( x\right ) is constant not zero. The above reduces to y^{\prime \prime }+\frac{y}{A^{2}}=0 This is harmonic oscillator It has the form of y^{\prime \prime }+\omega ^{2}y=0 with \omega =\frac{1}{A} being the natural frequency. The solution to this is easily found to be \begin{align} y\left ( x\right ) & =C_{1}\cos \left ( \omega x\right ) +C_{2}\sin \left ( \omega x\right ) \nonumber \\ & =C_{1}\cos \left ( \frac{x}{A}\right ) +C_{2}\sin \left ( \frac{x}{A}\right ) \tag{1} \end{align}
Since A is not constant, then we can try a similar solution1 but use f\left ( x\right ) for the arguments of the trigonometric functions\begin{equation} y\left ( x\right ) =C_{1}\cos \left ( f\left ( x\right ) \right ) +C_{2}\sin \left ( f\left ( x\right ) \right ) \tag{2} \end{equation} where f\left ( x\right ) is function of x to be determined. Hence. From now on, we will write f instead of f\left ( x\right ) to simplify notation.\begin{align*} y^{\prime } & =-C_{1}f^{\prime }\sin \left ( f\right ) +C_{2}f^{\prime }\cos \left ( f\right ) \\ y^{\prime \prime } & =-C_{1}f^{\prime \prime }\sin \left ( f\right ) -C_{1}\left ( f^{\prime }\right ) ^{2}\cos \left ( f\right ) +C_{2}f^{\prime \prime }\cos \left ( f\right ) -C_{2}\left ( f^{\prime }\right ) ^{2}\sin \left ( f\right ) \end{align*}
Substituting these back into the original ODE gives\begin{multline*} A^{2}\left ( -C_{1}f^{\prime \prime }\sin \left ( f\right ) -C_{1}\left ( f^{\prime }\right ) ^{2}\cos \left ( f\right ) +C_{2}f^{\prime \prime }\cos \left ( f\right ) -C_{2}\left ( f^{\prime }\right ) ^{2}\sin \left ( f\right ) \right ) +\\ AA^{\prime }\left ( -C_{1}f^{\prime }\sin \left ( f\right ) +C_{2}f^{\prime }\cos \left ( f\right ) \right ) +C_{1}\cos \left ( f\right ) +C_{2}\sin \left ( f\right ) =0 \end{multline*} Collecting terms gives \cos \left ( f\right ) \left ( -C_{1}A^{2}\left ( f^{\prime }\right ) ^{2}+C_{2}A^{2}f^{\prime \prime }+C_{2}AA^{\prime }f^{\prime }+C_{1}\right ) +\sin \left ( f\right ) \left ( -C_{1}A^{2}f^{\prime \prime }-C_{2}A^{2}\left ( f^{\prime }\right ) ^{2}-C_{1}AA^{\prime }f^{\prime }+C_{2}\right ) =0 Since this is zero for all \sin and \cos then \begin{align*} -C_{1}A^{2}\left ( f^{\prime }\right ) ^{2}+C_{2}A^{2}f^{\prime \prime }+C_{2}AA^{\prime }f^{\prime }+C_{1} & =0\\ -C_{1}A^{2}f^{\prime \prime }-C_{2}A^{2}\left ( f^{\prime }\right ) ^{2}-C_{1}AA^{\prime }f^{\prime }+C_{2} & =0 \end{align*}
Multiplying the first equation by C_{2} and the second by C_{1} gives\begin{align*} -C_{2}C_{1}A^{2}\left ( f^{\prime }\right ) ^{2}+C_{2}^{2}A^{2}f^{\prime \prime }+C_{2}^{2}AA^{\prime }f^{\prime }+C_{1}C_{2} & =0\\ -C_{1}^{2}A^{2}f^{\prime \prime }-C_{1}C_{2}A^{2}\left ( f^{\prime }\right ) ^{2}-C_{1}^{2}AA^{\prime }f^{\prime }+C_{1}C_{2} & =0 \end{align*}
Subtracting the second equation from the first gives\begin{align*} \left ( -C_{2}C_{1}A^{2}\left ( f^{\prime }\right ) ^{2}+C_{2}^{2}A^{2}f^{\prime \prime }+C_{2}^{2}AA^{\prime }f^{\prime }+C_{1}C_{2}\right ) -\left ( -C_{1}^{2}A^{2}f^{\prime \prime }-C_{1}C_{2}A^{2}\left ( f^{\prime }\right ) ^{2}-C_{1}^{2}AA^{\prime }f^{\prime }+C_{1}C_{2}\right ) & =0\\ -C_{2}C_{1}A^{2}\left ( f^{\prime }\right ) ^{2}+C_{2}^{2}A^{2}f^{\prime \prime }+C_{2}^{2}AA^{\prime }f^{\prime }+C_{1}C_{2}+C_{1}^{2}A^{2}f^{\prime \prime }+C_{1}C_{2}A^{2}\left ( f^{\prime }\right ) ^{2}+C_{1}^{2}AA^{\prime }f^{\prime }-C_{1}C_{2} & =0\\ C_{2}^{2}A^{2}f^{\prime \prime }+C_{2}^{2}AA^{\prime }f^{\prime }+C_{1}^{2}A^{2}f^{\prime \prime }+C_{1}^{2}AA^{\prime }f^{\prime } & =0\\ f^{\prime \prime }\left ( C_{2}^{2}A^{2}+C_{1}^{2}A^{2}\right ) +f^{\prime }\left ( C_{2}^{2}AA^{\prime }+C_{1}^{2}AA^{\prime }\right ) & =0 \end{align*}
Let us call C_{2}^{2}A^{2}+C_{1}^{2}A^{2}=\mu and C_{2}^{2}AA^{\prime }+C_{1}^{2}AA^{\prime }=\beta for the moment. The above becomes \mu f^{\prime \prime }+\beta f^{\prime }=0 Since f is missing, then we can solve the above by assuming f^{\prime }=v. The above becomes v^{\prime }+\frac{\beta }{\mu }v=0. This is linear in v. The integrating factor is I=e^{\int \frac{\beta }{\mu }dx}. Hence the ode becomes \begin{align*} \frac{d}{dx}\left ( e^{\int \frac{\beta }{\mu }dx}v\right ) & =0\\ v & =C_{3}e^{-\int \frac{\beta }{\mu }dx} \end{align*}
Since the proposed solution in (2) contains integration of constants already, we can choose C_{3}=1 without affecting the final solution. Hence f^{\prime }\left ( x\right ) =e^{-\int \frac{\beta }{\mu }dx} Therefore\begin{align} f\left ( x\right ) & =\int e^{-\int \frac{\beta }{\mu }dx}dx+C_{4}\nonumber \\ & =\int \left ( C_{3}e^{-\int \frac{C_{2}^{2}AA^{\prime }+C_{1}^{2}AA^{\prime }}{C_{2}^{2}A^{2}+C_{1}^{2}A^{2}}dx}dx\right ) dx+C_{4} \tag{3} \end{align}
Again, since the proposed solution in (2) contains integration of constants already, we can choose C_{4}=0. The above becomes\begin{align*} f\left ( x\right ) & =\int e^{-\int \frac{\beta }{\mu }dx}dx\\ & =\int e^{-\int \frac{C_{2}^{2}AA^{\prime }+C_{1}^{2}AA^{\prime }}{C_{2}^{2}A^{2}+C_{1}^{2}A^{2}}dx}dx \end{align*}
The expression \frac{C_{2}^{2}AA^{\prime }+C_{1}^{2}AA^{\prime }}{C_{2}^{2}A^{2}+C_{1}^{2}A^{2}} can be simplified as follows \frac{C_{2}^{2}AA^{\prime }+C_{1}^{2}AA^{\prime }}{C_{2}^{2}A^{2}+C_{1}^{2}A^{2}}=\frac{AA^{\prime }\left ( C_{2}^{2}+C_{1}^{2}\right ) }{A^{2}\left ( C_{2}^{2}+C_{1}^{2}\right ) }=\frac{AA^{\prime }}{A^{2}} Hence (3) becomes\begin{align*} f\left ( x\right ) & =\int e^{-\int \frac{AA^{\prime }}{A^{2}}dx}dx\\ & =\int e^{-\int \frac{A^{\prime }}{A}dx}dx\\ & =\int e^{-\ln A}dx\\ & =\int \frac{1}{A\left ( x\right ) }dx \end{align*}
Therefore the solution from (2) is\begin{align} y\left ( x\right ) & =C_{1}\cos \left ( f\left ( x\right ) \right ) +C_{2}\sin \left ( f\left ( x\right ) \right ) \nonumber \\ & =C_{1}\cos \left ( \int \frac{1}{A\left ( x\right ) }dx\right ) +C_{2}\sin \left ( \int \frac{1}{A\left ( x\right ) }dx\right ) \tag{4} \end{align}
Let us now try to verify this solution by substituting it back into the ODE. From (4), where we now write A instead of A\left ( x\right ) everywhere to simplify the notation\begin{align*} y^{\prime }\left ( x\right ) & =-C_{1}\sin \left ( \int \frac{1}{A}dx\right ) \left ( \int \frac{1}{A}dx\right ) ^{\prime }+C_{2}\cos \left ( \int \frac{1}{A}dx\right ) \left ( \int \frac{1}{A}dx\right ) ^{\prime }\\ & =-C_{1}\sin \left ( \int \frac{1}{A}dx\right ) \frac{1}{A}+C_{2}\cos \left ( \int \frac{1}{A}dx\right ) \frac{1}{A} \end{align*}
And y^{\prime \prime }\left ( x\right ) becomes\begin{multline*} y^{\prime \prime }\left ( x\right ) =-C_{1}\left ( \cos \left ( \int \frac{1}{A}dx\right ) \left ( \int \frac{1}{A}dx\right ) ^{\prime }\frac{1}{A}+\sin \left ( \int \frac{1}{A}dx\right ) \left ( \frac{-A^{\prime }}{A^{2}}\right ) \right ) \\ +C_{2}\left ( -\sin \left ( \int \frac{1}{A}dx\right ) \left ( \int \frac{1}{A}dx\right ) ^{\prime }\frac{1}{A}+\cos \left ( \int \frac{1}{A}dx\right ) \left ( \frac{-A^{\prime }}{A^{2}}\right ) \right ) \end{multline*} or\begin{align*} y^{\prime \prime }\left ( x\right ) & =-C_{1}\left ( \cos \left ( \int \frac{1}{A}dx\right ) \frac{1}{A^{2}}-\sin \left ( \int \frac{1}{A}dx\right ) \frac{A^{\prime }}{A^{2}}\right ) +C_{2}\left ( -\sin \left ( \int \frac{1}{A}dx\right ) \frac{1}{A^{2}}-\cos \left ( \int \frac{1}{A}dx\right ) \frac{A^{\prime }}{A^{2}}\right ) \\ & =-C_{1}\cos \left ( \int \frac{1}{A}dx\right ) \frac{1}{A^{2}}+C_{1}\sin \left ( \int \frac{1}{A}dx\right ) \frac{A^{\prime }}{A^{2}}-C_{2}\sin \left ( \int \frac{1}{A}dx\right ) \frac{1}{A^{2}}-C_{2}\cos \left ( \int \frac{1}{A}dx\right ) \frac{A^{\prime }}{A^{2}} \end{align*}
Substituting the above expression for y,y^{\prime },y^{\prime \prime } into the original ODE A^{2}y^{\prime \prime }+AA^{\prime }y^{\prime }+y=0 gives\begin{multline*} A^{2}\left ( -C_{1}\cos \left ( \int \frac{1}{A}dx\right ) \frac{1}{A^{2}}+C_{1}\sin \left ( \int \frac{1}{A}dx\right ) \frac{A^{\prime }}{A^{2}}-C_{2}\sin \left ( \int \frac{1}{A}dx\right ) \frac{1}{A^{2}}-C_{2}\cos \left ( \int \frac{1}{A}dx\right ) \frac{A^{\prime }}{A^{2}}\right ) \\ +AA^{\prime }\left ( -C_{1}\sin \left ( \int \frac{1}{A}dx\right ) \frac{1}{A}+C_{2}\cos \left ( \int \frac{1}{A}dx\right ) \frac{1}{A}\right ) +C_{1}\cos \left ( \int \frac{1}{A\left ( x\right ) }dx\right ) +C_{2}\sin \left ( \int \frac{1}{A\left ( x\right ) }dx\right ) =0 \end{multline*} Simplifying gives\begin{multline*} -C_{1}\cos \left ( \int \frac{1}{A}dx\right ) +C_{1}\sin \left ( \int \frac{1}{A}dx\right ) A^{\prime }-C_{2}\sin \left ( \int \frac{1}{A}dx\right ) -C_{2}\cos \left ( \int \frac{1}{A}dx\right ) A^{\prime }\\ -C_{1}A^{\prime }\sin \left ( \int \frac{1}{A}dx\right ) +C_{2}A^{\prime }\cos \left ( \int \frac{1}{A}dx\right ) +C_{1}\cos \left ( \int \frac{1}{A}dx\right ) +C_{2}\sin \left ( \int \frac{1}{A\left ( x\right ) }dx\right ) =0 \end{multline*} Canceling C_{1}terms gives -C_{2}\sin \left ( \int \frac{1}{A}dx\right ) -C_{2}\cos \left ( \int \frac{1}{A}dx\right ) A^{\prime }+C_{2}A^{\prime }\cos \left ( \int \frac{1}{A}dx\right ) +C_{2}\sin \left ( \int \frac{1}{A}dx\right ) =0 Which simplifies to -C_{2}\cos \left ( \int \frac{1}{A}dx\right ) A^{\prime }+C_{2}A^{\prime }\cos \left ( \int \frac{1}{A}dx\right ) =0 Or 0=0 Solution (4) has been verified.
Problem Find the general real solution to the equation xy^{\prime \prime }+\frac{3}{x}y=1+x^{3} Solution
We start by writing the ODE as\begin{equation} x^{2}y^{\prime \prime }+3y=x+x^{4}\tag{1} \end{equation} The solution is given by y=y_{h}+y_{p} where y_{h} is solution to homogeneous ODE x^{2}y_{h}^{\prime \prime }+3y_{h}=0 and y_{p} is a particular solution to x^{2}y_{p}^{\prime \prime }+3y_{p}=x+x^{4}. We start by solving the homogeneous x^{2}y^{\prime \prime }+3y=0 This is Euler type ODE. Using the standard substitution y=Ax^{r}, then y^{\prime }=Arx^{r-1},y^{\prime \prime }=Ar\left ( r-1\right ) x^{r-2} and the above becomes\begin{align*} x^{2}Ar\left ( r-1\right ) x^{r-2}+3Ax^{r} & =0\\ Ar\left ( r-1\right ) x^{r}+3Ax^{r} & =0 \end{align*}
Since x^{r}\neq 0 and A\neq 0 then the above simplifies to\begin{align*} r\left ( r-1\right ) +3 & =0\\ r^{2}-r+3 & =0 \end{align*}
Hence \begin{align*} r & =\frac{-b}{2a}\pm \frac{1}{2a}\sqrt{b^{2}-4ac}\\ & =\frac{1}{2}\pm \frac{1}{2}\sqrt{1-12}\\ & =\frac{1}{2}\pm \frac{1}{2}i\sqrt{11} \end{align*}
Hence the solution is \begin{align*} y & =C_{1}x^{\frac{1}{2}+\frac{i}{2}\sqrt{11}}+C_{2}x^{\frac{1}{2}-\frac{i}{2}\sqrt{11}}\\ & =C_{1}x^{\frac{1}{2}}x^{\frac{i}{2}\sqrt{11}}+C_{2}x^{\frac{1}{2}}x^{\frac{-i}{2}\sqrt{11}}\\ & =C_{1}\sqrt{x}e^{\ln x^{\frac{i}{2}\sqrt{11}}}+C_{2}\sqrt{x}e^{\ln x^{\frac{-i}{2}\sqrt{11}}}\\ & =C_{1}\sqrt{x}e^{\frac{i}{2}\sqrt{11}\ln x}+C_{2}\sqrt{x}e^{\frac{-i}{2}\sqrt{11}\ln x} \end{align*}
Using Euler formula the above can now be written in terms of \sin and \cos \begin{align} y & =\sqrt{x}\left ( C_{1}e^{\frac{i}{2}\sqrt{11}\ln x}+C_{2}e^{\frac{-i}{2}\sqrt{11}\ln x}\right ) \nonumber \\ y_{h} & =\sqrt{x}\left ( C_{3}\cos \left ( \frac{1}{2}\sqrt{11}\ln x\right ) +C_{4}\sin \left ( \frac{1}{2}\sqrt{11}\ln x\right ) \right ) \tag{2} \end{align}
Now we find the particular solution using the method of undetermined coefficients. Since the RHS is polynomial x+x^{4} then we guess y_{p}=A+Bx+Cx^{2}+Dx^{3}+Ex^{4} Then y^{\prime }=B+2Cx+3Dx^{2}+4Ex^{3} and y^{\prime \prime }=2C+6Dx+12Ex^{2}. Substituting these back in (1) \begin{align*} x^{2}\left ( 2C+6Dx+12Ex^{2}\right ) +3\left ( A+Bx+Cx^{2}+Dx^{3}+Ex^{4}\right ) & =x+x^{4}\\ 2Cx^{2}+6Dx^{3}+12Ex^{4}+3A+3Bx+3Cx^{2}+3Dx^{3}+3Ex^{4} & =x+x^{4}\\ 3A+x\left ( 3B\right ) +x^{2}\left ( 2C+3C\right ) +x^{3}\left ( 6D+3D\right ) +\left ( 3E+12E\right ) x^{4} & =x+x^{4}\\ 3A+x\left ( 3B\right ) +x^{2}\left ( 5C\right ) +x^{3}\left ( 9D\right ) +15Ex^{4} & =x+x^{4} \end{align*}
By comparing coefficients the following equations are generated\begin{align*} A & =0\\ 3B & =1\\ 5C & =0\\ 9D & =0\\ 15E & =1 \end{align*}
Hence A=0,B=\frac{1}{3},C=0,D=0,E=\frac{1}{15}. Therefore y_{p}=\frac{1}{3}x+\frac{1}{15}x^{4} Hence the final solution is\begin{align*} y & =y_{h}+y_{p}\\ & =\sqrt{x}\left ( C_{3}\cos \left ( \frac{1}{2}\sqrt{11}\ln x\right ) +C_{4}\sin \left ( \frac{1}{2}\sqrt{11}\ln x\right ) \right ) +\frac{1}{3}x+\frac{1}{15}x^{4} \end{align*}
Problem For what values of k does the equation\begin{equation} y^{\prime \prime }-\left ( \frac{1}{4}+\frac{k}{x}\right ) y=0 \tag{1} \end{equation} defined for 0<x<\infty have a solution vanishing at x=0 and at x=\infty ?
Solution
Let us look what happens at x\rightarrow \infty , then the term \frac{1}{4}\gg \frac{k}{x} and the ODE simplifies to y^{\prime \prime }-\frac{1}{4}y=0 Which has the solutions y=\left \{ e^{\frac{1}{2}x},e^{\frac{-1}{2}x}\right \} . We reject the first one since it does not vanish at x\rightarrow \infty , and use y=e^{\frac{-1}{2}x}. Now we assume the solution to (1) is of the form \begin{equation} y=P\left ( x\right ) e^{\frac{-1}{2}x} \tag{2} \end{equation} And we now try to find P\left ( x\right ) . Substituting this solution back into (1), given that
\begin{align*} y^{\prime } & =P^{\prime }e^{\frac{-x}{2}}-\frac{1}{2}Pe^{\frac{-x}{2}}\\ y^{\prime \prime } & =P^{\prime \prime }e^{\frac{-x}{2}}-\frac{1}{2}P^{\prime }e^{\frac{-x}{2}}-\frac{1}{2}P^{\prime }e^{\frac{-x}{2}}+\frac{1}{4}Pe^{\frac{-x}{2}}\\ & =P^{\prime \prime }e^{\frac{-x}{2}}-P^{\prime }e^{\frac{-x}{2}}+\frac{1}{4}Pe^{\frac{-x}{2}} \end{align*}
Substituting the above into (1) and canceling common term e^{\frac{-x}{2}} gives\begin{align} \left ( P^{\prime \prime }-P^{\prime }+\frac{1}{4}P\right ) -\left ( \frac{1}{4}+\frac{k}{x}\right ) P & =0\nonumber \\ P^{\prime \prime }-P^{\prime }-\frac{k}{x}P & =0\nonumber \\ xP^{\prime \prime }-xP^{\prime }-kP & =0 \tag{3} \end{align}
To solve this for P\left ( x\right ) , we use Frobenius series . Assuming\begin{align*} P\left ( x\right ) & =\sum _{n=0}^{\infty }c_{n}x^{n+r}\\ P^{\prime }\left ( x\right ) & =\sum _{n=0}^{\infty }\left ( n+r\right ) c_{n}x^{n+r-1}\\ P^{\prime \prime }\left ( x\right ) & =\sum _{n=0}^{\infty }\left ( n+r\right ) \left ( n+r-1\right ) c_{n}x^{n+r-2} \end{align*}
Hence (3) becomes\begin{align*} x\sum _{n=0}^{\infty }\left ( n+r\right ) \left ( n+r-1\right ) c_{n}x^{n+r-2}-x\sum _{n=0}^{\infty }\left ( n+r\right ) c_{n}x^{n+r-1}-k\sum _{n=0}^{\infty }c_{n}x^{n+r} & =0\\ \sum _{n=0}^{\infty }\left ( n+r\right ) \left ( n+r-1\right ) c_{n}x^{n+r-1}-\sum _{n=0}^{\infty }\left ( n+r\right ) c_{n}x^{n+r}-k\sum _{n=0}^{\infty }c_{n}x^{n+r} & =0\\ \sum _{n=0}^{\infty }\left ( n+r\right ) \left ( n+r-1\right ) c_{n}x^{n+r-1}-\sum _{n=1}^{\infty }\left ( n+r-1\right ) c_{n-1}x^{n+r-1}-k\sum _{n=1}^{\infty }c_{n-1}x^{n+r-1} & =0 \end{align*}
For n=0, and assuming c_{0}\neq 0 then\begin{align*} \left ( n+r\right ) \left ( n+r-1\right ) c_{n} & =0\\ \left ( r\right ) \left ( r-1\right ) c_{0} & =0\\ r\left ( r-1\right ) & =0 \end{align*}
Hence r=1 or r=0.
Case r=1\begin{align*} P & =\sum _{n=0}^{\infty }c_{n}x^{n+1}\\ & =\sum _{n=1}^{\infty }c_{n-1}x^{n} \end{align*}
Hence\begin{align*} P^{\prime } & =\sum _{n=1}^{\infty }nc_{n-1}x^{n-1}\\ P^{\prime \prime } & =\sum _{n=1}^{\infty }\left ( n\right ) \left ( n-1\right ) c_{n-1}x^{n-2} \end{align*}
And now (3) becomes\begin{align*} x\sum _{n=1}^{\infty }\left ( n\right ) \left ( n-1\right ) c_{n-1}x^{n-2}-x\sum _{n=1}^{\infty }nc_{n-1}x^{n-1}-k\sum _{n=1}^{\infty }c_{n-1}x^{n} & =0\\ \sum _{n=1}^{\infty }\left ( n\right ) \left ( n-1\right ) c_{n-1}x^{n-1}-\sum _{n=1}^{\infty }nc_{n-1}x^{n}-k\sum _{n=1}^{\infty }c_{n-1}x^{n} & =0\\ \sum _{n=2}^{\infty }\left ( n\right ) \left ( n-1\right ) c_{n-1}x^{n-1}-\sum _{n=1}^{\infty }nc_{n-1}x^{n}-k\sum _{n=1}^{\infty }c_{n-1}x^{n} & =0\\ \sum _{n=1}^{\infty }\left ( n+1\right ) \left ( n\right ) c_{n}x^{n}-\sum _{n=1}^{\infty }nc_{n-1}x^{n}-k\sum _{n=1}^{\infty }c_{n-1}x^{n} & =0 \end{align*}
Hence for n\geq 1 we obtain\begin{align*} \left ( n+1\right ) \left ( n\right ) c_{n}-nc_{n-1}-kc_{n-1} & =0\\ c_{n} & =\frac{\left ( n+k\right ) c_{n-1}}{n\left ( n+1\right ) } \end{align*}
For n=1 c_{1}=\frac{\left ( k+1\right ) c_{0}}{2} For n=2 c_{2}=\frac{\left ( k+2\right ) c_{1}}{2\left ( 3\right ) }=\frac{\left ( k+2\right ) }{2\left ( 3\right ) }\frac{\left ( k+1\right ) }{2}c_{0}=\frac{\left ( k+1\right ) \left ( k+2\right ) }{\left ( 2\right ) \left ( 2\right ) \left ( 3\right ) }c_{0} For n=3 c_{3}=\frac{\left ( k+3\right ) c_{2}}{3\left ( 4\right ) }=\frac{\left ( k+3\right ) }{\left ( 3\right ) \left ( 4\right ) }\frac{\left ( k+1\right ) \left ( k+2\right ) }{\left ( 2\right ) \left ( 2\right ) \left ( 3\right ) }c_{0}=\frac{\left ( k+1\right ) \left ( k+2\right ) \left ( k+3\right ) }{\left ( 2\right ) \left ( 2\right ) \left ( 3\right ) \left ( 3\right ) \left ( 4\right ) }c_{0} For n=4 c_{4}=\frac{\left ( k+4\right ) c_{n-1}}{\left ( 4\right ) \left ( 5\right ) }=\frac{\left ( k+4\right ) c_{3}}{\left ( 4\right ) \left ( 5\right ) }=\frac{\left ( k+1\right ) \left ( k+2\right ) \left ( k+3\right ) \left ( k+4\right ) }{\left ( 2\right ) \left ( 2\right ) \left ( 3\right ) \left ( 3\right ) \left ( 4\right ) \left ( 4\right ) \left ( 5\right ) }c_{0} And so on. Hence\begin{align} P\left ( x\right ) & =\sum _{n=1}^{\infty }c_{n-1}x^{n}\nonumber \\ & =c_{0}x+c_{1}x^{2}+c_{2}x^{3}+c_{3}x^{4}+\cdots \nonumber \\ & =c_{0}\left ( x+\frac{\left ( k+1\right ) }{2}x^{2}+\frac{\left ( k+1\right ) \left ( k+2\right ) }{\left ( 2\right ) \left ( 2\right ) \left ( 3\right ) }x^{3}+\frac{\left ( k+1\right ) \left ( k+2\right ) \left ( k+3\right ) }{\left ( 2\right ) \left ( 2\right ) \left ( 3\right ) \left ( 3\right ) \left ( 4\right ) }x^{4}+\frac{\left ( k+1\right ) \left ( k+2\right ) \left ( k+3\right ) \left ( k+4\right ) }{\left ( 2\right ) \left ( 2\right ) \left ( 3\right ) \left ( 3\right ) \left ( 4\right ) \left ( 4\right ) \left ( 5\right ) }x^{5}+\cdots \right ) \nonumber \\ & =c_{0}\left ( x+\left ( k+1\right ) \frac{x^{2}}{2!}+\frac{\left ( k+1\right ) \left ( k+2\right ) }{2!}\frac{x^{3}}{3!}+\frac{\left ( k+1\right ) \left ( k+2\right ) \left ( k+3\right ) }{3!}\frac{x^{4}}{4!}+\frac{\left ( k+1\right ) \left ( k+2\right ) \left ( k+3\right ) \left ( k+4\right ) }{4!}\frac{x^{5}}{5!}+\cdots \right ) \tag{4} \end{align}
But e^{x}=1+x+\frac{x^{2}}{2!}+\frac{x^{3}}{3!}+\cdots . Or e^{x}-1=x+\frac{x^{2}}{2!}+\frac{x^{3}}{3!}+\cdots . So there is an exponential term inside (4). Hence to make (4) vanish at x\rightarrow \infty , then k needs to be a negative integer. Taking k=-1 makes all terms with k in them vanish, leaving P\left ( x\right ) =c_{0}x So now the solution from (2) becomes y\left ( x\right ) =c_{0}xe^{-\frac{x}{2}} Which goes to zero as x\rightarrow \infty since an exponential decays to zero faster that x going to infinity.
We now need to check if negative k integer value (specifically k=-1 which we picked from above) will also make the solution vanish as x\rightarrow 0. When x\rightarrow 0 the ODE becomes\begin{equation} y^{\prime \prime }-\frac{k}{x}y=0 \tag{5} \end{equation} Since \frac{k}{x}\gg \frac{1}{4} close to x=0. Since k is negative integer -1 then the above becomes y^{\prime \prime }+\frac{k}{x}y=0 To see this will go to zero as x\rightarrow 0, Intuitively since \frac{k}{x} is now positive and very large, then this is like a harmonic oscillator with very large stiffness. (Spring mass system). When the stiffness becomes very large, the solution goes to zero (the natural frequency goes to infinity, since \omega =\sqrt{\frac{k}{x}} which means the period goes to zero since \omega =2\pi T) which implies no motion. So this shows that negative integer value of k found from first part makes the solution vanish at both x\rightarrow \infty and at x\rightarrow 0. Actually for x\rightarrow 0 we just needed k to be negative in order to change the sign. But for x\rightarrow \infty we found we needed k to be a negative integer which we choose -1. So this will work for x=0 and x=\infty .
I first tried to solve the give ODE directly using series method. I left this here as an appendix, not to be graded but as a reference.
x is singular point. But it is a regular singular point since \lim _{x\rightarrow 0}x^{2}\frac{k}{x}=x and hence the limit exist. Therefore assuming solution is Frobenius series y=x^{r}\sum _{n=0}^{\infty }c_{n}x^{n}=\sum _{n=0}^{\infty }c_{n}x^{n+r} Therefore y^{\prime }=\sum _{n=0}^{\infty }\left ( n+r\right ) c_{n}x^{n+r-1} and y^{\prime \prime }=\sum _{n=0}^{\infty }\left ( n+r\right ) \left ( n+r-1\right ) c_{n}x^{n+r-2}, then (1) becomes\begin{align*} \sum _{n=0}^{\infty }\left ( n+r\right ) \left ( n+r-1\right ) c_{n}x^{n+r-2}-\left ( \frac{1}{4}+\frac{k}{x}\right ) \sum _{n=0}^{\infty }c_{n}x^{n+r} & =0\\ \sum _{n=0}^{\infty }\left ( n+r\right ) \left ( n+r-1\right ) c_{n}x^{n+r-2}-\frac{k}{x}\sum _{n=0}^{\infty }c_{n}x^{n+r}-\frac{1}{4}\sum _{n=0}^{\infty }c_{n}x^{n+r} & =0\\ \sum _{n=0}^{\infty }\left ( n+r\right ) \left ( n+r-1\right ) c_{n}x^{n+r-2}-k\sum _{n=0}^{\infty }c_{n}x^{n+r-1}-\frac{1}{4}\sum _{n=0}^{\infty }c_{n}x^{n+r} & =0 \end{align*}
But k\sum _{n=0}^{\infty }c_{n}x^{n+r-1}=k\sum _{n=1}^{\infty }c_{n-1}x^{n+r-2} and \sum _{n=0}^{\infty }c_{n}x^{n+r}=\sum _{n=2}^{\infty }c_{n-2}x^{n+r-2} and the above becomes\begin{equation} \sum _{n=0}^{\infty }\left ( n+r\right ) \left ( n+r-1\right ) c_{n}x^{n+r-2}-k\sum _{n=1}^{\infty }c_{n-1}x^{n+r-2}-\frac{1}{4}\sum _{n=2}^{\infty }c_{n-2}x^{n+r-2}=0 \tag{2} \end{equation} The first step is to obtain the indicial equation. As the nature of the roots will tell us how to proceed. The indicial equation is obtained from n=0 in (2) with the assumption that c_{0}\neq 0. This leads to\begin{align*} \left ( n+r\right ) \left ( n+r-1\right ) c_{n} & =0\\ r\left ( r-1\right ) c_{0} & =0 \end{align*}
c_{0} is always taken as non-zero. This leads to r\left ( r-1\right ) =0 With solutions r_{1}=1 or r_{2}=0. (We take r_{1} as the larger root first, since Frobenius series solution can only guarantee solution for the larger root, when the roots differ by an integer as this is the case).
Since r_{1}-r_{2} is an integer, then this tells us we can obtain a first solution y_{1}\left ( x\right ) associated with r_{1}=1 from the Frobenius series\begin{equation} y_{1}\left ( x\right ) =\sum _{n=0}^{\infty }c_{n}x^{n+1} \tag{3} \end{equation} But to find the second solution y_{2}\left ( x\right ) associated with r_{2}=0 we can try either reduction of order method or use \begin{equation} y_{2}\left ( x\right ) =Ay_{1}\left ( x\right ) \ln \left ( x\right ) +\sum _{n=0}^{\infty }d_{n}x^{n} \tag{4} \end{equation} Where A is some constant, which can be zero, and d_{n} are the coefficients for the second series. We have to do the above when the roots of the indicial equation differ by integer. Otherwise, the second solution would have been found using Frobenius series y_{2}\left ( x\right ) \sum _{n=0}^{\infty }c_{n}x^{n+r_{2}} like with the first solution.
OK, Now we will first find y_{1}\left ( x\right ) from (3)
case r_{1}=1
Using (3)\begin{align*} y^{\prime } & =\sum _{n=0}^{\infty }\left ( n+1\right ) c_{n}x^{n}\\ y^{\prime \prime } & =\sum _{n=0}^{\infty }n\left ( n+1\right ) c_{n}x^{n-1}\\ & =\sum _{n=1}^{\infty }n\left ( n+1\right ) c_{n}x^{n-1} \end{align*}
Substituting the above into (1) gives\begin{align*} \sum _{n=1}^{\infty }n\left ( n+1\right ) c_{n}x^{n-1}-\left ( \frac{1}{4}+\frac{k}{x}\right ) \sum _{n=0}^{\infty }c_{n}x^{n+1} & =0\\ \sum _{n=1}^{\infty }n\left ( n+1\right ) c_{n}x^{n-1}-\frac{1}{4}\sum _{n=0}^{\infty }c_{n}x^{n+1}-\frac{k}{x}\sum _{n=0}^{\infty }c_{n}x^{n+1} & =0\\ \sum _{n=1}^{\infty }n\left ( n+1\right ) c_{n}x^{n-1}-\frac{1}{4}\sum _{n=0}^{\infty }c_{n}x^{n+1}-k\sum _{n=0}^{\infty }c_{n}x^{n} & =0\\ \sum _{n=0}^{\infty }\left ( n+1\right ) \left ( n+2\right ) c_{n+1}x^{n}-\frac{1}{4}\sum _{n=1}^{\infty }c_{n-1}x^{n}-k\sum _{n=0}^{\infty }c_{n}x^{n} & =0 \end{align*}
For n=0\begin{align*} \left ( 1\right ) \left ( 2\right ) c_{1}-kc_{0} & =0\\ c_{1} & =\frac{k}{2}c_{0} \end{align*}
For n>0 we obtain the recursion equation\begin{align*} \left ( n+1\right ) \left ( n+2\right ) c_{n+1}-\frac{1}{4}c_{n-1}-kc_{n} & =0\\ c_{n+1} & =\frac{\frac{1}{4}c_{n-1}+kc_{n}}{\left ( n+1\right ) \left ( n+2\right ) } \end{align*}
For n=1 c_{2}=\frac{\frac{1}{4}c_{0}+kc_{1}}{\left ( 2\right ) \left ( 3\right ) }=\frac{\frac{1}{4}c_{0}+k\left ( \frac{k}{2}c_{0}\right ) }{6}=\frac{\frac{1}{4}c_{0}+\frac{k^{2}}{2}c_{0}}{6}=c_{0}\frac{\frac{1}{4}+\frac{k^{2}}{2}}{6}=c_{0}\frac{1+2k^{2}}{24} For n=2\begin{align*} c_{3} & =\frac{\frac{1}{4}c_{1}+kc_{2}}{\left ( 3\right ) \left ( 4\right ) }\\ & =\frac{\frac{1}{4}\frac{k}{2}c_{0}+kc_{0}\frac{1+2k^{2}}{24}}{12}\\ & =\frac{\frac{k}{8}c_{0}+kc_{0}\frac{1+2k^{2}}{24}}{12}\\ & =c_{0}\frac{\frac{3k+k+2k^{3}}{24}}{\left ( 12\right ) }\\ & =c_{0}\frac{4k+2k^{3}}{288} \end{align*}
And so on. Hence \begin{align*} y_{1}\left ( x\right ) & =c_{0}x+c_{1}x^{2}+c_{2}x^{3}+c_{3}x^{4}+\cdots \\ & =c_{0}x+\frac{k}{2}c_{0}x^{2}+c_{0}\frac{1+2k^{2}}{24}x^{3}+c_{0}\frac{4k+2k^{3}}{288}x^{4}+\cdots \\ & =c_{0}x\left ( 1+\frac{k}{2}x+\frac{1+2k^{2}}{24}x^{2}+\frac{4k+2k^{3}}{288}x^{3}+\cdots \right ) \\ & =c_{0}x\left ( 1+\frac{k}{2}x+\left ( \frac{1}{24}+\frac{1}{12}k^{2}\right ) x^{2}+\frac{k}{288}\left ( 4+2k^{2}\right ) x^{3}+\cdots \right ) \end{align*}
I could not find closed form function for the above.
Now that we found y_{1}\left ( x\right ) , then y_{2}\left ( x\right ) is, from (4), repeated here\begin{equation} y_{2}\left ( x\right ) =Ay_{1}\left ( x\right ) \ln \left ( x\right ) +\sum _{n=0}^{\infty }d_{n}x^{n} \tag{4} \end{equation} Since we want the solution to vanish at x=0 then we set A=0 and y_{2}\left ( x\right ) simplifies to\begin{equation} y_{2}\left ( x\right ) =\sum _{n=0}^{\infty }d_{n}x^{n} \tag{4} \end{equation} Where d_{0}\neq 0. Hence y^{\prime }\left ( x\right ) =\sum _{n=0}^{\infty }nd_{n}x^{n-1} and y^{\prime \prime }=\sum _{n=0}^{\infty }n\left ( n-1\right ) d_{n}x^{n-2}. Rewriting the ODE as xy^{\prime \prime }-\left ( \frac{x}{4}+k\right ) y=0 and now substituting the derivatives into this gives\begin{align*} x\sum _{n=0}^{\infty }n\left ( n-1\right ) d_{n}x^{n-2}-\left ( \frac{x}{4}+k\right ) \sum _{n=0}^{\infty }d_{n}x^{n} & =0\\ \sum _{n=0}^{\infty }n\left ( n-1\right ) d_{n}x^{n-1}-\frac{x}{4}\sum _{n=0}^{\infty }d_{n}x^{n}-k\sum _{n=0}^{\infty }d_{n}x^{n} & =0\\ \sum _{n=2}^{\infty }n\left ( n-1\right ) d_{n}x^{n-1}-\frac{1}{4}\sum _{n=0}^{\infty }d_{n}x^{n+1}-k\sum _{n=0}^{\infty }d_{n}x^{n} & =0\\ \sum _{n=1}^{\infty }\left ( n+1\right ) \left ( n\right ) d_{n+1}x^{n}-\frac{1}{4}\sum _{n=1}^{\infty }d_{n-1}x^{n}-k\sum _{n=0}^{\infty }d_{n}x^{n} & =0 \end{align*}
For n=0 we obtain kd_{0}=0 which implies d_{0}=0 since k\neq 0.
For n>0\begin{align*} \left ( n+1\right ) \left ( n\right ) d_{n+1}-\frac{1}{4}d_{n-1}-kd_{n} & =0\\ d_{n+1} & =\frac{\frac{1}{4}d_{n-1}+kd_{n}}{\left ( n\right ) \left ( n+1\right ) } \end{align*}
For n=1 d_{2}=\frac{\frac{1}{4}d_{0}+kd_{1}}{2}=\frac{k}{2}d_{1} For n=2 d_{3}=\frac{\frac{1}{4}d_{1}+kd_{2}}{\left ( 2\right ) \left ( 3\right ) }=\frac{\frac{1}{4}d_{1}+k\left ( \frac{k}{2}d_{1}\right ) }{6}=\frac{d_{1}+2k^{2}d_{1}}{32}=d_{1}\frac{1+2k^{2}}{32} For n=3 d_{4}=\frac{\frac{1}{4}d_{2}+kd_{3}}{\left ( 3\right ) \left ( 4\right ) }=\frac{d_{2}+4kd_{3}}{48}=\frac{\frac{k}{2}d_{1}+4k\left ( d_{1}\frac{1+2k^{2}}{32}\right ) }{48}=d_{1}\frac{\frac{1}{8}k\left ( 2k^{2}+5\right ) }{48}=d_{1}\frac{\left ( 2k^{3}+5k\right ) }{384} And so on. Hence the second solution is\begin{align*} y_{2}\left ( x\right ) & =\sum _{n=0}^{\infty }d_{n}x^{n}\\ & =d_{0}+d_{1}x+d_{2}x^{2}+d_{3}x^{3}+d_{4}x^{4}+\cdots \\ & =d_{1}x+\frac{k}{2}d_{1}x^{2}+d_{1}\frac{1+2k^{2}}{32}x^{3}+d_{1}\frac{\left ( 2k^{3}+5k\right ) }{384}x^{4}+\cdots \\ & =d_{1}x\left ( 1+\frac{k}{2}x+d_{1}\frac{1+2k^{2}}{32}x^{2}+d_{1}\frac{\left ( 2k^{3}+5k\right ) }{384}x^{3}+\right ) \end{align*}
I am not sure if the above solution for y_{2}\left ( x\right ) is correct. I need to check this again later.