Hello there yesterday in my exam I had the following problem:
Knowing that $y_1=x$ is a solution, solve: $$(x^2 - 1)y''+2xy'-2y=0$$ My try was to use Liouville-Ostrogradski formula (see: https://www.encyclopediaofmath.org/index.php/Liouville-Ostrogradski_formula) So bassically what I have is (I will denote with C the second Wronskian since is just a constant): $$y_1y_2'-y_1'y_2=Ce^{-\int{\frac{2x}{x^2-1}}dx}$$ $$xy_2'-y_2=\frac{C}{x^2-1}\rightarrow \frac{y_2'}{x}-\frac{y_2}{x^2}=\frac{C}{x^2(x^2-1)}$$ $$\left(\frac{y_2}{x}\right)'=\frac{C}{(x^2-1)x^2}=\frac{C}{(x^2-1)}-\frac{C}{x^2}$$ By integrating : $$\frac{y_2}{x}=\left(\frac{C}{2}\ln\left|\frac{x-1}{x+1}\right|+\frac{C}{x}+C_2\right)$$$$y_2=c_1x\ln\left|\frac{x-1}{x+1}\right|+c_2x +c_3$$ Was my solution correct and complete? I am wondering if I must prove that this differential equation has only those two solutions, but I have no ideea how.
Concerning completeness, for a linear ODE like this, given any two solutions $f_1$ and $f_2$, any linear combinations of them, $Af_1 + Bf_2$ is also going to be a solution. Therefore the set of all solutions of the ODE form a vector space.
It is a little harder to show, but if the ODE is of order $n$, then the space of solutions is $n$ dimensional. That is 2-dimensional in your case.
So to find all solutions, you just need to find 2 linearly independent solutions. All the other solutions will be a linear combination of these two.
You are given that one solution is $y_1 = x$. From your work (with the constant correction that has been discussed by others), a second solution is $$y_2 = x\ln\left|\frac{x-1}{x+1}\right| + 2$$ Therefore the general solution is $$y = Ax + B\left(x\ln\left|\frac{x-1}{x+1}\right| + 2\right)$$ for arbitrary constants $A,B$.